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Some Applications of the Wave Front Shearing 
Interferometer 


James B. Saunders 


(May 27, 


eter. 


1965) 


This paper gives the results of several applications of the wave front shearing prism interferom- 
The instrument is very compact and easy to apply. 


It is applied to the testing of chromatic 


aberration of simple and compound lenses; and to the testing of wave forms that characterize the 


monochromatic aberrations (spherical, coma, and astigmatism). 


ferent type lenses. 
telescope objectives. 


Key Words: 


chromatic aberration of lenses. 


Previous applications of lens testing interferometers 
have been limited because of the difficulties associ- 
ated with their use. The Twyman lens testing inter- 
ferometer, for instance, requires several precision 
elements that have apertures equal to or exceeding 
that of the lens to be tested. The curvature of the 
convex reference mirror should match the back focus 
of the lens to a reasonable approximation, thus neces- 
sitating the use of a different reference mirror for 
lenses of different focal lengths. The relatively large 
separation of the elements invites vibrations and is 
vulnerable to temperature and atmospheric disturb- 
ances. Several modifications of the wave front shear- 
ing prism interferometer (abbreviated to WFSPI), that 


AA “i 


FIGURE 1. 
Interference is « 
leave the source in different directions 


»btained by the combination of one component 


Results are shown for several dif- 


This interferometer is equally applicable to the testing of small lenses and large 


Interferometer, testing of lenses, lens aberrations, aberrations, prism interferometer, 


are small and relatively free from some of the above 
mentioned objections, have been described by 
Saunders [1].' 

All the elements of the WFSPI are combined into 
one small compound prism. Since it operates at a 
point of convergence of the light beam (fig. 1) it can 
be very small. A small prism can be made so that 
imperfections in it are quite negligible. A prism 2.5 
em in length will accommodate any lens that provides 
a working distance (image point to back surface of lens) 
in excess of 1.7 cm if its aperture does not exceed 
this distance. 


: Pane sin bei kets indicate the literature references at the end of this paper 


A, & BD 


Optics of the interferometer. 
wr 4 


of ray 4 with one component, BB’), of ray B that 
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The maximum separation of any two interfering 
beams is the lateral shear, S, which is usually a small 
fraction of the diameter of the lens being tested. 
Thus, the effects of temperature and atmospheric 
disturbances are easily suppressed. 

Interference fringes may be obtained with white 
light. However, a narrow band interference filter is 
preferred for most work. 

The WFSPI divides a wave front into two compo- 
nents, shears one component relative to the other, 
and recombines the two beams of light so as to pro- 
duce interference in the overlapping area. Each ray 
that enters the prism is divided into four components, 
two emerging from one face and two from another, as 
shown in figure lb. None of these four components 
are recombined and, therefore, do not interfere with 
each other. Interference is obtained by the combi- 
nation of one component of each of two rays that leave 
the source in different directions (rays A and B) and 
are subsequently recombined by the prism, as shown 
in figure le. An observer’s eye, receiving either pair 
of the recombined rays, sees interference fringes in 
whatever plane it is focused on. When the eye is 
focused on the lens, it sees fringes in the plane of the 
lens (fig. 1d). 

A reference surface, such as is required in the 
Twyman interferometer, is not used since this inter- 
ferometer yields results by comparing a wave front 
with an image of itself. The shape and width of the 
fringes are measures of the deviation of the original 
wave front from a sphere. A simple mathematical 
operation [2] has been described that permits a com- 
putation of the wave front from observations of a 
single fringe pattern. 

The fringe configuration produced by a lens that is 
afflicted with chromatic aberration will change with 
change in the wavelength of the light source. Figure 
2 shows a plot of image distance (focal length of the 
lens if the source is at infinity) versus wavelength for 
a single element lens that has a focal length of approxi- 
mately 42 cm and an aperture of 9cm. The slit of a 
prism monochromator was used, with a zirconium arc, 
as the light source. A small filament[3], mounted on a 
milling machine table and applied as described by 
Platzeck [4] and Gaviola, was also used to locate the 
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WAVELENGTH 


FIGURE 2. 


onvex lens 


Image distance versus wavelength. 
Plano-« fecal length 42 cm, aperture 9 em 


The inserted photographs are 
the fringe patterns for the indicated wavelengths 


focal positions. The photographs inserted in figure 
2 represent the fringe configurations for the three 
indicated wavelengths. 

It can be shown by ray tracings that the image dis- 
tance, q, corresponding to any wavelength, A,, is 
given by the formula 


qo =d+an;'—A,N,d/dbn,, 


where d is the distance from the lens to the prism (fig. 
la); a is the distance from the entrance face of the 
prism to the point where the principal ray suffers 
total internal reflection; n, is the refractive index of 
the prism; N, and 6 are, respectively, the number of 
fringes and distance between the two reference points, 
P! and P; (fig. 3). The sign of N, is chosen so that it 
increases when d is increased. The angle ¢ is the 
relative deviation, inside the prism, of any two com- 
ponents of an original ray that emerges from the same 
face of the prism (such as rays A; and A; of fig. 1b). 
The difference in focal distance, Aqi2, for any two 
wavelengths, A; and Ag, is given by 


Aqie — a(n, 1— ny 1)— AAiNin “. = AoNon my \/db. 


A change of one in the value of N corresponds to 
approximately 0.5 mm in the focal distance for the 
prism used while obtaining the data for figure 2. This 
prism produces a shear of 0.006 radians. 

When the chromatic aberration is small, a prism of 
high sensitivity (large shear angle) should be used. 
Figure 4 is a plot of image distance versus wavelength 
for an air spaced doublet (focal length=463 mm, 
aperture =64 mm). The shear is approximately 0.047 
radians. One fringe corresponds to approximately 
0.1 mm shift in the image point. The photographs 
from which the data were obtained are shown in 
figure 5. 

The shape of a wave front, produced by a lens, is 
representative of all monochromatic aberrations of the 
lens in whatever manner it is used when producing 
this wave form. The interferometer of figure la pro- 
duces a set of interference fringes that is characteristic 
of the shape of the wave front. The photograph shown 
in figure 6 depicts fringes produced by an f/2 camera 
lens with the source located 20 m from the lens. The 
curve marked f/2 in figure 7 is a plot of the deviations 
of the corresponding wave front from a best fitting 
reference sphere. It is apparent that if the between- 
the-lens diaphragm had been adjusted to an f/3 system 
the fringe pattern would have been reduced in area 
to that of the overlapping area of the two inserted 
circles in figure 6. When this reduced fringe pattern 
is analyzed by the method described in reference 2, 
the computed deviations from a best fitting sphere are 
found to be represented by curve f/3 in figure 7. 

Since the curve that represents the deviation of the 
wave front from a best fitting sphere also represents 
the image forming characteristics of the lens, it is 
suggested that the quality of any lens may be repre- 
sented by a relatively simple number. The manner of 
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FIGURE 3. 


Reference points P{ and Ps are 


used for measuring 


change in focus. 


Points P) and PP 


separation of reference 


and P 
points P 
beam in the plane of the lens 


lens 














FIGURE Image distance wavelength for a 2-element 


achromatic lens. 


versus 


selecting this number may be best decided later, but 
it appears that the weighted arithmetical mean devi- 
ation of the wave front from a statistically chosen best 
fitting sphere might represent the image forming 
characteristics of the lens. This number would vary 
with the position of the object, relative to the lens and 
its axis; however, so does the quality of the image. 
The weights should be chosen so as to give equal 
weight per unit area over the whole lens. The chosen 
number should also include other factors, such as 


bare the two sheared images of a reference point PP. The 
ind P, equals (S +b 


where S is the lateral shear of the 


Phe reference points may be marked on the surface of the 


FIGURE 5. Photograph of interference fringes for different wave- 


lengths. 


aperture for instance, and might therefore be a com- 
pound number. 

The weighted mean deviations of wave fronts are 
readily obtained from the corresponding fringe patterns 


produced by this interferometer. The weighted arith- 
metical mean deviation for curve f/2 in figure 7 is 2.1 
and for f/3 it is 0.1A. The relative magnitude of these 
two numbers is believed to represent, to a close approx- 
imation, the corresponding relative image forming 
qualities. 
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FIGURE 6. Interference fringes in a f]2 camera lens. 


The inserted circles represent the edge of an iris diaphragm when adjusted to transmit 
an {/3 beam of light 


— — 





( 


-5 ( 5 
DISTANCE FROM CENTER 


FIGURE 7. Deviations of the wave front, produced by an f]/2 lens 
when wide open and at f/3, from best fitting reference spheres. 


The photographs of figure 8 are the fringes produced 
by an aerial camera lens of 8-in. aperture and 50-in. 
focal length that was designed and made at NBS. 
This lens has four air-spaced elements, and was 
designed to be used with another element —a concave 
lens near the focus—as an image flattener. The direc- 
tions of shear are indicated for each photograph. Only 
four of the eight photographs that were made are shown 
here. Figure 9 shows the eight curves, representing 
the deviations from best fitting spheres along the indi- 
cated diameters. The eight equally spaced directions 
are indicated in the drawing. The abscissas of each 
of the four pairs of curves in figure 9 represent the 
same direction along the indicated diameter. The 
difference in wave front shape, for the several diam- 
eters, clearly indicates assymmetry in the wave front. 

An astronomical objective may be tested in the 
laboratory by autocollimation if a sufficiently large 
optical flat is available. If such a flat is not avail- 
able, it may be tested in the manner illustrated in 
figure la, by placing a small source at a considerable 
distance from the lens. Figure 10 shows the photo- 
graphs of the fringes and figure 11 shows the corre- 


180° — 9° 270°—90° 


AERIAL CAMERA LENS FOCAL LENGTH 50" APERTURE 6” 


FIGURE 8. Interference fringes of a lens when sheared in opposite 
directions along two axes that are normal to each other. 


O°- 180° 





~ ~A-45°- 225° 
225°-45°-- “ 


YP | 


DEVIATIONS 











CENTER , cm 


DISTANCE FROM 


FIGURE 9. Curves representing deviations of wave fronts from close 
fitting spheres. 
The direction of shear. relative to the initial position of the lens, was increased 45° be 
tween successive exposures for making the corresponding photographs 
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FIGURE 10. Fringes produced by a lens when used properly, 


and when the axis is reversed, b. 
Aperture equals 12 in.; fo« 


il length equals 15 ft.: air spaced, two element objective 








STANCE ALONG DIAMETER,cm 


FIGURE | Deviations of front versus distance along one 
diameter of an astronomical refractor. 


computed from the two 


wave 


The prese lata photographs of figure LO 
sponding deviation curves for a 12-in. aperture, 15-ft 
focal length refractor, when it is used as intended and 
when the axis is reversed. 


Figure 12 shows interference fringes produced by 
the 32-in. aperture, 40-ft focal length reflecting tele- 
scope at the University of Virginia at. Charlottesville, 


Va. The 12-sec exposure used for this photograph 
practically elimates the atmospheric effect by record- 
ing the average position of the fringes. The photo- 
graph is too recent to have permitted inclusion of the 
data analysis in this report. 


FIGURE 12. First photograph of interference fringes produced by a 


star for testing an astronomicai reflector. 


Aperture 32 in.: focal length equal 40 ft. Cassigrainian type telescope 


Test made with this interferometer show that it is 
as easy to assemble as a Foucault knife edge test. 
It gives quantitative results that are as accurate as 
other interferometer tests. The fringe pattern is prac- 
tically independent of vibrations of the instrument. 
Technical personnel are not required for its use. It 
is relatively easy to make, the cost of material for it is 
negligible, and it is easy to reduce the effects of tem- 
perature gradients and atmospheric turbidity to negli- 
gible effects, and it is easy to reduce the effects of tem- 
perature gradients and atmospheric turbidity to negli- 
gible significance. Because of the simplicity of its 
operation, size, ruggedness and stability, this inter- 
ferometer readily lends itself to combination with 
electronics for automatic testing for lenses. 
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Precision Method for Evaluating Primary Aberrations 
of Lenses With a Twyman Interferometer 


James B. Saunders 


(May 21, 1965) 


A simplified procedure for separate evaluation of the three primary aberrations of lenses with the 


Twyman lens testing interferometer is described. 


of reference points that give best results. 
in the literature. 


Each of the aberration coefficients is found to be a 
function of observations at only four points on the lens. 


Equations are given for the optimum choice 


These equations are applied to data previously reported 
rhe results indicate that this procedure is sufficiently precise to reveal high order 


aberrations in a lens that was assumed to be practically free from such aberrations. 


Key Words: Interferometer, lens testing, Twyman interferometer. 


1. Introduction 


The formula for optical path difference in the 
Twyman lens testing interferometer, used by Pro- 
fessor Kingslake, contains seven parameters. It 
would seem that data from seven reference points 
would be necessary for evaluating any one of the 
parameters. By judiciously choosing the location of 
the reference points each of the primary monochro- 
matic aberrations may be evaluated from data taken 
at only four points. It will be shown that the results 
obtained from four points, with a given lens, are of the 
same order of accuracy as those previously obtained, 
by least squares, from 44 points. This represents a 
very large reduction in computations. 

Each of the primary monochromatic aberrations is 
separately evaluated. The results indicate that the 
sensitivity of the Twyman interferometer is sufficient, 
not only to measure the aberration coefficients ac- 
curately, but also to measure changes in these aber- 
rations with changes in the aperture of the lens. 

The discovery that four properly chosen points are 
adequate for evaluation of each aberration coefficient 
with the Twyman interferometer was the result of an 
anologous discovery that two points are adequate for 
the same measurements with the Kosters prism inter- 
ferometer. 


2. Fundamental Equations 


The terminology previously used by Kingslake ! will 
be used in this paper. The data which he obtained 


with a simple lens? will be used for demonstration of 


the procedures described here. This lens was as- 
sumed to be practically free from higher order aber- 
rations over the aperture used for obtaining results. 

'R. Kingslake, The interferometer patterns due to the primary aberrations, Trans. Opt 


Soc. 27, 94 (1925-6) 
2 R. Kingslake, The analysis of an interferogram, Trans. Opt. Soc 28, | (1926-7) 


The formula for optical path difference is 
A(x? + y*)? + Byix? + y*) + Cir? + 3y*) 
+ Dix?+y?)+Ey+Fx+G=p () 


where A, B, and C represent, respectively, coefficients 
of spherical aberration, coma, and astigmatism. 
These coefhicients are related to the aberration con- 
stants, a;, a2, and a3, by the equations, 


Longitudinal (primary) 
spherical aberration 


4A=a,= sf 





__ coma 


FBIh => sth 


_ distance between focal lines 


2f-C/h? = as oe 





where s is half the aperture of the lens, A is the image 
distance from the optic axis and f is the distance from 
the image to the back surface of the lens. The xy- 
rectangular coordinate system, used in eq (1), is chosen 
so that the xyplane is in the principal plane of the 
lens—with the y-axis in the meridional plane. The 
origin is on the optic axis of the lens. The coefficients, 
D, E, F, and G relate to adjustments of the interferom- 
eter. The order of interference is G, at the origin 
(on the axis of the lens), and is p at any point (x, y). 

Figure 1 shows a photograph of the interferogram 
used by Kingslake in footnote 2. This interferogram 
contains the effects of any and all aberrations that 
are associated with this lens when it is used in the 
manner described. If this interferogram should con- 
tain higher order aberrations, eq (1) could not repre- 
sent the fringes. The parameters, A, B, and C would 
be found to vary with the aperture of the lens; whereas, 
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FicurRE 1. Interference fringes of a simple ophthalmic lens as seen 


in a Twyman lens testing interferometer. 
ircle represents the 


The inscribed « “useful aperture of the lens.” 


This area was pre 
viously assumed to be practically free from higher order aberrations 


about the margin, represent regions of positive and negative orders of interference 


if the lens is free from higher order aberrations these 
parameters would be independent of lens aperture 
and, therefore, also independent of the location of the 
chosen reference points. 

The family of curves in figure 2 represents the inter- 
ferogram of a lens that is absolutely free from higher 


FIGURE 2. 


The (+) and (—) signs, 


order aberrations. These curves (fringes) are com- 
puted from eq (1) and represent the area inside the 
circle in figure 1. The values for the parameters are 
those obtained by Kingslake (see footnote 2). Each 
curve corresponds to an integral order of interference 
(integral value for p). 


3. Choice of Reference Points 


Data from a family of at least four reference points 
are required to obtain a unique value for each of the 
aberration coefficients, 4, B, and C in eq (1). The 
relative distribution of these points (fig. 3) can be so 
chosen that when the observed data (values for x, y, 
and p) at four points are substituted for their counter- 
parts in eq (1), the resultant four equations will permit 
the elimination of all parameters except the one de- 
sired. There are many different arrangements of 
these four reference points that will permit solutions 
of the aberration parameters. Preliminary study 
showed that, in each case, one such arrangement of 
points could be found that yields optimum accuracy. 
Therefore, only one arrangement of points will be 
given here for evaluating each of the three aberration 
parameters. These arrangements are which 
give optimum accuracy. 


those 


Graphical representation of computed fringes inside 


the circle of figure 1. 


These curves were 


that contains only primary aberrations 


computed from the optical path difference formula (eq (1) for a lens 
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4. Astigmatism 


If we draw a circle of radius r, with its center at 
the origin, in figure 3, it will intersect the coordinate 
axes at (r, 0) (0, r),(—r, 0), and (0,—r. If we choose 
these four points in the interferogram of figure 1, we 
may evaluate the parameter C from data observed at 
these points. 

The observed data are the orders of interference, 
p, and the coordinates, x and y, at the four reference 
points. The data, labeled (11, yi, pi), (%2, v2, po), 
(x3, ¥3, Ps), and (x4, ys, ps), may be considered as the 
coordinates in a three dimension rectangular coordi- 
nate system. We get optimum accuracy when the 


following relationships exist between the x and y 


coordinates: 
Ls = 
Cf ta 


On substituting these values in eq (1), successively 
for the four points, we obtain the following four 
equations: 


rA + P°C+PrD +rF +G=p, 


A+ PB+3PC4+PD+4+rE +G=p2 


r4A + PC+PD —rF +G=ps3 


rA—PB+3rC+P°D—rE +G=ps 
If we add the first equation in 3 to the third and the 
second to the fourth, we obtain, 


2Art + 2Cr + 2Dr + 2G = pit ps 
: (4) 
2Ar*+ 6Cr + 2Dr? + 2G = pot ps 


By subtracting the first member equation of 4 from 
the second member we obtain 


4Cr? = po+ pa— pi— Ps. (5) 


Values for x and y, along the two axes, are given by 
Kingslake for integral values of p. These values are 
represented graphically in figure 4. The values for 
the p’s in table | are interpolated from figure 4. When 
the p’s, that correspond to any chosen value for r, 
are substituted in eq (5) we obtain corresponding values 
for C. These values are given in column 6, table 1. 
The probable error in C is given in column 7 as AC, 
and is based on the assumption that the probable 
error in a single observation is the same as that ob- 
tained by Kingslake. The probable error, AP, of a 
single observation was found to be 0.013A. Values 
for C are plotted in figure 5 as ordinates for corre- 
sponding values of r as abscissas. The magnitude of 
the probable error, AC, is represented by the length 
of arrows attached to the corresponding points. Ob- 
viously, the most precise value for C is obtained by 
using large values for r. 
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FIGURE 3. The distribution of reference points for most precise 
evaluation of primary aberrations. 

The four larger circles, on the axes, are for astigmatism; the four smaller circles, on the 

x-axis, are for spherical aberration; and the four solid disks, on the y-axis, are for coma 
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FIGURE 4. Cross sections of the wave front along the x and y-axes. 
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TABLE 1. Value for astigmatism, C, and probable error, AC, 
wavelength) for different lens apertures (in millimeters) 


| 
Orders of interferenc Coefhcient of 


astigmatism 
AC 


0.650 
163 
072 
O41 
026 
O18 
013 
ol 
008 
OO7 


Since the fringes in figure 2 are precisely repre- 
sented by eq (1), the values of the corresponding 


parameters should be independent of the choice of 


reference points. The value for C (labeled Cy in fig. 5) 
obtained from any family of four reference points in 
figure 2 will be the same as that found by Kingslake 
because this value was used for computing these 
curves. The corresponding probable error in C;, 
(i.e., AC;) given by Kingslake is also indicated in figure 
5 to illustrate the relative probable error for results 
obtained from four favorably chosen reference points 
and from 44 reference points treated statistically by 
the Method of Least Squares. Since the coefficient, 
C, varies with r (i.e., from one family of reference 
points to another) eq (1) obviously does not apply to 
the fringe pattern of figure 1. The lens that produced 
this set of fringes is obviously afflicted with higher 
order aberrations. 








8 
RADIUS VECT 


OR.r (mm) 
FIGURE 5. Values of the coefficient of astigmatism, C 
aperture, r, of the lens. 


The vertical arrows represent computed probable errors in ¢ 
line represents the value for C obtained from figure 2 anc 


,» versus semi- 


The horizontal dashed 
1 given in footnote 2 


5. Coma 


Coma cannot be evaluated from a family of reference 
points that are all equally distant from the origin, as 
is the case for astigmatism. However, if we choose a 
family of four points on the y-axis (fig. 3) made up of 
two pairs such that the two points of each pair are 
equally distant from the origin but of opposite sign, 
then we can obtain an equation for B as a function 
of the coordinates of these four points. If the ordi- 
nates of the inner pair of points are one half the value 
of those for the outer points and the outer points are 
near the margin of the lens the resultant value for B 
will be the most precise value obtainable from any 
family of four points. The marginal point should be 
far enough from the margin to permit an accurate 
evaluation of the corresponding order of interference. 

The observed x, y, p-coordinates for the four points 
are (X2, Y2, P2), (%5, ¥5, Ps), (X6, Ye, Pe), and (x4, ¥4, P4), 
where 


Xo = Xs —=X%e—X4—0 | 


2ys vaiaae 26 a r| 


We may substitute the coordinates for each of these 
four points into eq (1) and obtain four equations, from 
which we may obtain the following formula for B: 


3P°-B = 2(p2— ps+ 2p6— 2ps). 
If we compute values for B that correspond to dif- 


ferent values of r (i.e., for different apertures of the 
lens) we obtain the results shown in figure 6. The 


A, WAVELENGTH 








Figure 6. Values of the coefficients of coma, B (on the left) and 
spherical aberration, A (on the right) versus the semiaperture, 
r, of the lens. 

The vertical arrows represent computed probable errors. The horizontal, dashed lines 

represent values obtained from figure 2 and given in footnote 2 
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probable error, AB, of B (based on the same assump- 
tions used for obtaining AC, above) is indicated by the 
vertical arrows attached to the corresponding points. 
The value of B, computed from figure 2, (indicated 
as Bx in fig. 5), is constant for all families of reference 
points. This value for Bx is the same as that given 
by Kingslake for the same lens. 


6. Spherical Aberration 


Spherical aberration cannot be evaluated from a 
family of reference points that are equally distant 
from the origin, as is the case for astigmatism; nor 
can it be evaluated from a family of four points on the 
y-axis, as is the case for coma. However, if we choose 
a family of four points on the x-axis, made up of two 
pairs such that the two points of each pair are equally 
distant from some convenient point, x (see fig. 3) 
but in opposite directions from xo, then we can obtain 
an equation for A as a function of the coordinates of 
these four points alone. Also, if the distance from 
xo to each of the inner pair of points (x7 and xx) are 
one-half the value of those for the outer points (x; 
and xy), and the outermost point, x;, is near the mar- 
gin of the lens, the resultant value for A will be the 
most precise value obtainable from any such family 
of points. 

The x, y, p-coordinates of the four points are (x, 
¥1, Pr), (X7, ¥7, Pz), (Xs, Ys, Ps), and (xy, ¥9, py), where 


X; = (8/5)x7 = 4x9 = — 8xg = — 2xg= Fr, 


vi} 0 


This family of points has the same relationship to 
xo as the family of points in eq (5) has to the origin 
(optic axis). 

On successively substituting the coordinates of 
these four points into eq (1) we obtain four equations 
from which we may obtain the following formula for A: 


81Ar*= 128[p: — py + 2(ps — pz). 


Values for A, computed for different values of r (i.e., 
for different apertures) are shown in figure 6, together 
with the corresponding probable errors, AA. The 
value for A, computed from figure 2 (indicated as Ax 
in fig. 6), is constant for all values of r. The value 
Ax is that given by Kingslake for the same lens. 


7. Conclusions 


The method for obtaining the coefficients of primary, 
monochromatic aberrations, described here, requires 
a very small amount of calculations. A 10-in. slide 
rule is adequate for most work. 

The coefficients of all three of the aberrations are 
found to vary with the aperture of the lens. This in- 
dicated that the method is sufficiently precise to eval- 
uate even higher orders of aberrations provided the 
necessary formula for optical path difference is 
available. 


(Paper 69C4—205) 
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Comparators for Voltage Transformer Calibrations at NBS 
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(July 28, 1965) 


An inductive comparator is described for measuring ratio-correction-factor and phase-angle error 
of voltage transformers by comparison with reference transformers. The simple circuit utilizes a 
resistance-capacitance network and an inductive voltage divider. The relative magnitude of the ratio- 
correction-factors is obtained from the dial readings of the inductive voltage divider, and the phase 
angle between the voltage phasors is equal to +wCR. The accuracy of these measurements at 60 
and 400 hertz is within 2 parts-per-million for ratio and within 10 microradians for phase angle at the 
limits of the phase-angle range of + 38 milliradians. Only a few ranges of reference transformers are 
required since fractional ratios from 1/1 to 1/4 may be measured as readily as ratios near unity, without 
diminution in accuracy. A general discussion of a resistive comparator, which has been in con- 
tinuous operation at NBS since 1939, is included. 


Key Words: Voltage divider, inductive voltage divider, resistive voltage divider, voltage trans- 


former calibration, ratio-correction-factor, phase-angle errors, voltage transformer 
comparator, voltage-transformer-test-set, comparator, test-set, calibration. 


1. Introduction 


The method most widely used by public utilities 
and standardizing laboratories to determine the ratio 
and phase angle of voltage transformers is by com- 
parison with calibrated reference transformers. How- 
ever, it is a relative rather than an absolute method, 
in that the accuracy of the results obtained depends 
upon the knowledge of accuracy of the reference 
transformers [1].! 

At the NBS laboratories, a group of reference 
transformers is available by means of which test trans- 
formers having any ratio from 1:1 to 2000:1 can be 
calibrated. The accurate initial 
the voltage ratio and phase angle for these reference 
transformers has been discussed elsewhere and will 
not be repeated here [2, 3]. 

The purposes of this paper are (1) to describe the 
resistive comparator which has been in continuous 
operation since 1939; (2) to describe a new inductive 
comparator; and (3) to analyze various factors affecting 
the accuracy with which measurements can be made. 


2. Resistive Comparator 


This comparator was designed and constructed at 
NBS by F. K. Harris and has been in continuous 
operation since 1939, except for a very brief period 
during which new switches were installed to replace 
the worn out original ones [4, 5]. 

The schematic circuit of this comparator is shown 
in figure 1.2 The in-phase voltage balance condition 
is achieved by adjusting the tap of the resistive di- 
vider, and the quadrature voltage, by varying the mu- 
tual inductor, M. 


' Figures in brackets indicate the literature references at the end of this paper. 
2 Simplified by not showing the shielding arrangements around each segment of R, 
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FIGURE 1. 


2.1. Description of Components 


The R, and R,, circuits are supplied from the sec- 
ondary of transformer S. Each can be varied from 
2000 to 4000 2 in increments of 100 0. The exact 
values are determined by the ratios of the test and 
reference transformers. The 100-0 and 1000-2 unit 
resistors of R, are of woven wire and mounted on 
ceramic forms. Thus the residual, L;, of Rs is re- 
duced to a minimum [6]. The R» circuit, nominally 
equal to Rs, include the primary circuit of M which 
has a d-c resistance of 200 2. The capacitor, C, 
shunting the lower portion of the resistance, serves 
to compensate for the self-inductance introduced by 
M [5]. Multiple taps are brought out from R» to sup- 
ply the shield potential of corresponding sections of 
R, and of the detector circuit [7]. The shielding is 
shown in figure 4. 

As shown in figure 2, the portion of the Ry, circuit 
between 1800 and 2000 2 carries a continuously ad- 
justable tap, R,, which is connected through the de- 
tector circuit and M to the secondary of the test 
transformer. 
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FIGURE 2. 


The mutual inductor used in this circuit consists 
of a 10 X 20 mH per section unit (+ 200 mH maximum) 
and a continuously adjustable inductometer (+25 
mH maximum) [8, 9]. 

Both the R,;/R, ratio of the resistive divider and 
the mutual inductance, M, are periodically redeter- 
mined to insure that maximum accuracy is maintained 
for transformer testing. The residual reactances, 
L, and L,, were also accurately determined at the 
time the comparator was built. 

The design is such that the resistors must dissipate 
considerable power. Therefore, the entire resistive 
divider, R; and R», is immersed in oil to minimize 
the effect of temperature coefficient on the R;/Rz 
ratio. 

The detector circuit consists of a doubly shielded 
impedance-matching transformer (1000-4000 turns) 
and a wave-analyzer used as a tuned null detector. 


2.2. Comparator Desk 


Figure 3 shows the general layout of the desk. The 
resistive divider is located to the right with the R, 
adjustment dials in the back and the R, dials in front. 
The step mutual inductor and the inductometer are 
located in the left portion of the desk. The detector 
is in the rear right. Since this comparator is adapt- 
able for calibrating voltage-transformer-test-sets, the 
panel shown at the lower left side of the desk is for 
arranging the connections for calibrating either trans- 
formers or test-sets. Figure 4 shows the internal 
wiring of the resistive divider and the individual 
shields around units of Rs. 


2.3. Theoretical Relations 


The equations for obtaining the ratio-correction- 
factor and phase-angle error are developed as follows. 
By definition, the voltage ratio of a transformer is the 
ratio of the primary terminal voltage to the secondary 


258 


FIGURE 3. 


terminal voltage so that 


| V pri | =N’= Nf 

|} Sec! ; 
where N’ is the actual voltage ratio, N the nominal 
ratio, and f the ratio-correction-factor. If @ represents 
the angle between the primary and reversed secondary 
voltage phasors, then 


’ Vy piss 
y,= Tol (cos 0,+] sin 6s) 


ore 


‘ 
y, = "41 (cos 6, +; sin 6;) 
Ns 


By inspection of figure 1, the following mesh current 
equations are obtained 

Ve=1N2Z54hZ 

Ve=1Z,+14{Z2+ Zs) + I3joM 

V,=InjoM + 13Zm 


3 Subscripts s and x denote that the symbols refer to the reference or test transformers, 
respectively 





FIGURE 4. 


1,=0 when a null is indicated on the detector, then 


# IZ; IsZ m 


Since ‘ea Rw by the addition of €. 
nominally, the balance equation becomes 


and Rn=R, 


(cos Oz +7 sin Or) (cos 0,+ 7 sin 0.) 





rf. Nets 
{Re |! +40 (-E +R) || 


fs cos 6, 
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. NR, 
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e+ i= D+ 8-6 
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| a © 1 , 
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R, R, | R 


As previously indicated, R, and R» circuits con- 
stitute a burden on the reference transformer. There- 
fore, f; and @, are determined with these burdens. 

The design of the divider is such that the balance 
conditions require 

N; 


N’ 


Re > l 
R, ol + 


this requirement is met by using reference transform- 
ers whose voltage ratios have been adjusted to values 
2 percent less than the integral values usually en- 
countered as nominal ratios in instrument 
transformers. 

The minimum resolution of this setup is such that 
the precision for ratio measurements is 2 ppm and for 
phase angle, 15 wrad at 25 and 60 Hz. However, sus- 
ceptibility of the mutual inductor to the gradients of 
stray magnetic fields, large residuals and uncertainty 
in ratio calibration of the resistive divider limit the 
accuracy of these measurements to 10 ppm and 60 
rad, respectively. These errors increase rapidly 
with frequency. Therefore, this comparator is not 
suitable for accurate calibrations beyond 60 Hz. 


3. Inductive Comparator 


This comparator avoids a number of the objection- 
able features of the resistive comparator. The shield- 
ing is simple and is easily accomplished: the method 
is capable of extension to cover the audio-frequency 
range; the apparatus is less bulky and requires very 
little power for operation; the components are in- 
herently more accurate and stable and readily avail- 
able commercially; the ratio difference is indicated 
directly on the decade dials; and the set (except the 
detector) is self-contained in one portable case. 

The schematic diagram of the circuit is shown in 
figure 5, which indicates the positions of resistor, R, 
and detector when the voltage phasor of the test trans- 
former, X, leads that of the reference transformer, S. 
For lagging angles, the positions of R and detector are 


*The (f,—1) and @, for the NBS reference transfowners are very small 
the approximations made are justifiable 


Therefore. 
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FIGURE 5. 





simply interchanged. The in-phase voltage balance 
is obtained by adjusting the tap on the inductive volt- 
age divider, and the quadrature voltage balance is 
obtained by adjusting either the resistor or the ca- 
pacitor. 


3.1. Theoretical Relations 


As previously defined, the voltage ratio is 


V, Q ms A 


“a. +j sin 6s), 


V.= ee 


y 
(cos 6,+j sin 6,)=— 
Nifr m 


m being the actual ratio of the inductive voltage di- 
vider. Corrections to the ratios and the phase angles 
of the inductive voltage divider used are insignificant 
and can be neglected. Then from figure 5, 


V._1—jwCR ene 
oe ee ee 


The balance equation of the circuit becomes 


N,fr(cos 0,+ j sin 0s) 
mN,f(cos 0,+] sin 0,) 





= 1—w?C*R?—jwCR 


and on separation and by neglecting second order 
terms 


Naft 
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radians, and 
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interchanged, 
6, ~ 6,—wCR radians, and 
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3.2. Loading Effects 


Since the input impedance of the inductive voltage 
divider at 60 Hz is in excess of 0.8 MQ (nearly 1.00 


power factor), the effect of loading of the f, and 6, 
of the test transformer can be expressed as 


to Transformer Secondary 
(External Impedance) 


(“ee Total Impedance Referred ) 





= (Correction to fy and 6,). 


A transformer with an equivalent resistance of 1 2 
(which is easily measurable) will result in a correction 
of approximately 1 ppm [5]. 

The effect of loading by the quadrature balancing 
network on either the inductive voltage divider or 
the reference transformer is usually insignificant. 
At the limits of phase-angle range of + 37.8 mrad (130 
min), the correction to ratio and phase angle of the 
inductive voltage divider is ~ 1 ppm and <4 ywrad, 
respectively; or, to the reference transformer, < < 1 
ppm and < <1 wrad. 


3.3. Description of Components 


All the components in this circuit are readily avail- 
able commercially or accessible in most testing lab- 
oratories. In figure 5, R is a 0-10 kQ, 10-turn, 
continuously adjustable resistor with small residual 
reactance. C is a 101 nF decade capacitor with a 
dissipation factor of < 0.0003. Both R and C were 
calibrated before and after assembly. The deviations 
from nominal values were found to be within 52 and 
5 pF, respectively. 

The inductive voltage divider is a six-decade unit 
designed for optimum performance at 60 Hz.> For 
tests at higher frequencies (up to 1 kHz), it can easily 
be replaced by a unit of appropriate design, or by 
application of appropriate loading corrections to the 
ratio-correction-factor and phase angle. However, 
due to the simplicity of the circuitry and infrequent 
requests for calibrations at higher frequencies, a tem- 


* Since maximum voltage on this divider is limited to 150 V at 60 Hz, it must be connected 
across the test transformer rather than the reference transformer whose secondary voltage 
may exceed this value in some instances 
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porary setup is generally used at NBS. This consists 
of a high frequency divider, commercially available 
air capacitors, and an a-c decade resistor. 

The input impedance of the inductive voltage di- 
vider was measured in a capacitive voltage divider 
setup: see figure 7 for circuit diagram [10]. The volt- 
age-ratio corrections and phase-angle errors were 
obtained ‘by a tracking method [11]. However, for Ficune 6. 
a well designed inductive voltage divider, these cor- 
rections are usually < 0.2 ppm of input in ratio lin- 
earity deviation and <8 wrad of input in phase 
deviation. 

The detector circuit consists of a wave-analyzer used 
as a tuned null indicator and a specially constructed 
shielded transformer, shown in figure 6. This trans- 
former is used to match impedances and to isolate 
the measuring circuit from the detector circuit electro- 
statically. It consists of a high permeability toroidal 
core, a 700-turn inner primary winding, two electrostatic 
shields, a 3500-turn outer secondary winding, and an 
overall Mu-metal case. The use of a toroid minimizes 
the magnetic pickup from stray fields in the laboratory. 

The polarity checking voltmeter is a 0-150 V, 3-in. 
panel-type instrument of 3 percent accuracy. The 
21 kQ) resistor matches its internal resistance. The 
voltmeter indicates the secondary voltage of the test 
transformer when the selector switch is in the 
“CALIBRATION” position. 

Figures 8 and 9 show the exterior and interior of 
the comparator. The cylindrical shield shown in 
the interior view encloses the 10-kQ resistor. 


3.4. Accuracy Determination 


After the completion of the comparator, the cali- 
bration of each component was rechecked. 
The accuracy of the completed comparator was 
verified by the following procedures: (1) connecting 
a single voltage source to both “V;” and “V,"’ ter- 
minals to obtain the zero corrections for C and R and 
the unity correction for dial settings of the inductive 
voltage divider; (2) calibrating several pairs of trans- 
formers of the same ratio to reference transformers of 
appropriate ratios, thus determining the relative ratio- : 
correction-factors and phase-angle errors with several FIGURE 9. 
values of secondary burden for each pair of transform- 
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ers; (3) comparing each pair of these calibrated trans- 
formers to each other with the comparator alternately 
as reference and test units and with various combina- 
tions of previously used secondary burdens; and (4) 
connecting a single voltage source to “V,”’ terminals 
and to the input of a calibrated inductive voltage 
divider, whose output is connected to “V,” terminals 
in order to check the fractional ratio measurements. 


3.5. Transformer Testing Procedure 


A reference transformer is selected so that V, > N,. 
The selector switch is set to “CHECK” position; and 
the ratio, m, is set to equal N,/Ns. The secondaries 
of the test and reference transformers are connected 
to “V,” and “V,” terminals respectively. With the 
transformers energized and the voltmeter switch de- 
pressed, the voltmeter indicates (V,;—V,)/2 if the po- 
larities of V,; and V, are correctly connected to the 
comparator. [The polarity of V, must be reversed if 
voltmeter indicates (V;+V,)/2.] The test proceeds 
with the selector switch turned to “CALIBRATION.” 
V, is now indicated on the voltmeter. 

If VN; > Ns by more than 10 percent, the test can still 
be made, but the reference transformer must be con- 
nected to “V,” terminal and the test transformer to 
“Vo. At the balanced condition, the ratio, m, will 
indicate the reciprocal of the relative ratio-correction- 
factors, and the relative phase-angle switch will show 
the negative of the markings as shown in figures 6 
and 8. 


3.6. Consideration of Errors 


An examination of errors requires a careful appraisal 
of all components of the measuring circuit as well as a 
justification of estimations, even though in this case 
one would expect most of the associated errors to be 
negligible. 

The first to be considered are those possible errors 
concealed through approximation in developing the 
expressions for obtaining the ratio-correction-factor 
and phase-angle relations, where the effects of the 
residual of the 10-kQ resistor and the dissipation factor 
of the decade capacitor are ignored. Approximate 
magnitudes of these errors can be obtained by adding 
the estimated equivalent residual reactance across R 
and estimated equivalent series resistance to C. 
These effects were found to be less than 0.5 ppm for 
ratio and 5 wrad for phase angle at the limits of the 
phase angle range (+ 38 mrad) at 60 Hz. 

A second source of uncertainty is in the calibration 
of the inductive voltage dividers. However, past 
experience indicates that errors should not exceed 
0.2 ppm of input in ratio linearity deviation and 2 wrad 
of input in phase deviation [12]. 

Errors caused by the leads, switch contacts, and 
imperfect shields are estimated to be <0.5 ppm and 
<2 prad. 

The last important uncertainty that should be con- 
sidered is in obtaining the transformer secondary 


load corrections. Such errors are estimated to be 
0.5 ppm and 3 prad for each transformer. 

Finally, one should add 3 times the standard devia- 
tions (0.3 ppm of ratio and 1] purad) of the results ob- 
tained in repeated balances of the comparator to the 
combination of the systematic-error estimates. 
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4. Conclusion 


A high accuracy inductive comparator has been 
developed for measuring the relative voltage ratio 
and phase angle of voltage transformers and testing 
sets at power and audio frequencies. An _ overall 
accuracy of <2 ppm in ratio measurements and < 10 
prad in phase angle at 60 and 400 Hz is achieved, 
with a slight reduction in the accuracy attainable 
at higher audio frequencies. 

The simple circuit does not require any special 
equipment. Only a few ranges of reference trans- 
formers are required, since fractional ratios from 1/1 
to 1/4 may be measured as readily as ratios near unity. 

This inductive comparator could replace testing 
sets now in general use. It is compact and portable, 
has a broader useful frequency range, and exceeds 
the accuracy of many commercial sets presently avail- 
able by one or two orders of magnitude. 

The resistive comparator will continue to be used 
at NBS for calibration of voltage transformers at 25 
Hz at which frequency the operating voltage is usually 
above the maximum operating voltage of the inductive 
comparator. 


The author acknowledges the valuable help and 
advice received from Mr. F. L. Hermach, Chief, 
Electrical Instrument Section, NBS. Mr. R. J. Berry 
contributed — skillfully in constructing the new 
comparator. 
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Methods are presented for measuring both the relative and absolute voltage dependence of admit- 


tance standards. 
which are described in detail. 


on selected capacitors are presented. 


The practical circuits based on the methods are composed of precision bridges 
Measurements with these circuits yield values for the proportional 
changes in admittance which result from specified changes in voltage. 


The results of measurements 


Uncertainties are of the order of one part in 109. 


\ general instability in the voltage dependence of precision air capacitors was observed at the 


higher accuracy levels. 
surtaces. 


Most of the instability is believed to be caused by changes in the electrode 


Key Words: Bridge, capacitor, coaxial choke, difference transformer, electrode surfaces, thin 
films, three-terminal, transformer, voltage coefhcient, voltage dependence. voltage 


ratio. 


1. Introduction 

Several years ago a computable cross capacitor, 
based on a theorem by A. M. Thompson and D. G. 
Lampard [1],! was constructed at NBS [2], The 
value of this capacitor is computable, in esu, within 
an estimated accuracy of 2 parts per million (ppm). 
A new cross capacitor of improved design is now being 
constructed. It is expected that this capacitor will 
be computable, in esu, within a few parts in 108 or 
better. 

The expected increase in accuracy not only requires 
better bridges and capacitors, but also requires that 
many sources of error be reevaluated. One such 
source of error is the voltage dependence of precision 
air capacitors. 

To meet our particular needs, the quantity chosen 
as a measure of voltage dependence is termed y and 
is defined by the equation 


Y2=¥(1+y) 


where Y; is the admittance of a given standard at 
voltage V;, Y2 is the admittance of the same standard 
at voltage V2, and the condition V2 > V, is imposed to 
restrict the meaning of y. A given standard may 
have any number of values of y, each with a speci- 
fied voltage change. Since y represents a propor- 
tional change in admittance (AY/Y), this measure of 
voltage dependence is in the preferred form for cor- 
rections to precision measurements. 

In the past, most measurements of voltage depend- 
ence have been confined to a determination of the 
ratio, (1+ ya)/(1+ yo), where ya and y» apply to differ- 
ent capacitors but have the same specified voltage 
change. One exception to the above is the work of 
N. L. Kusters and O. Petersons [3]. They developed 
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' Figures in brackets indicate the literature references at the end of this paper 


a method for measuring the voltage dependence of 
an individual capacitor, provided the voltage depend- 
ence results only from electrode deflections. 

The accuracy requirements described earlier re- 
sulted in the need for a method which would be inde- 
pendent of the mechanisms which cause voltage 
dependence. Such a method is described in section 2. 
The method yields values of y corresponding to a 
change in voltage from V to 2V, where V is a variable. 

The practical circuits based on the method are 
described in sections 3 and 4. These circuits were 
used to measure values of y for several capacitors 
corresponding to the following voltage changes: 12.5 
to 25, 25 to 50, 50 to 100, and 100 to 200 V, rms, at 
1592 Hz. Measurement uncertainties had to be kept 
small since the above values of y are used to calculate 
other values of y in different voltage ranges such as 
from zero to 200 V. The latter value is needed for 
the determination of the absolute unit of capacitance. 
Small uncertainties also provide the means for quickly 
detecting instabilities in the voltage dependence of 
capacitors. It was found that the voltage dependence 
of precision air capacitors was, in general, unstable 
at the accuracy level described. This subject is dis- 
cussed in section 5. 


2. Theory 


One method for measuring the voltage ratio of a 
transformer is to connect admittance standards to 
the transformer so as to form a bridge as shown in 
figure 1 [4]. In general, if the transformer has a 
nominally m:n ratio, then m+ n nominally equal admit- 
tances are used, and m+n detector balances are 
required corresponding to a cyclic permutation of 
these admittances. In this sequence of balances 
each admittance appears with voltage E; m times, 
and with voltage E. n times. If E; # E2, a ratio meas- 
urement error results since the admittance of a prac- 
tical standard changes when the voltage applied to it 
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FIGURE 1. Circuit for measuring the voltage ratio of a trans- 


former {4}. 


is changed. The difference between the measured 
ratio and the true ratio is related to the average change 
in admittance of the standards. 

If the transformer ratio can be determined from a 


second set of measurements which is independent of 


voltage, the combination of the two sets of measure- 
ments will yield a value for the average change in ad- 
mittance of the m+n standards corresponding to a 
change in voltage from E; to Ez. The change in admit- 
tance of each standard can then be determined by 
measuring the relative changes between various pairs 
of standards. 

The basis of this paper is a method of the type de- 
scribed above in which the second set of measurements 
consists of determining a group of 1:1 ratios which 
can be combined to determine the m:n ratio of a 
transformer. For simplicity, a transformer having a 
2:1 ratio was chosen so that only two 1:1 ratios need 
be determined. 

It is assumed that taps are brought out from the 
secondary winding of a bridge transformer [5] so as 
to obtain three secondary voltages of approximately 
equal magnitude. These voltages are represented 
by E,4, Ex, and Ec in figure 2. The circuit shown in 
figure 2 consists of three bridges, each of which may 
be balanced, in turn, by the small, adjustable voltage 
e. For the case in which the circuit elements are 
closely matched, e/E, =q <1, and only a low accu- 
racy calibration of q is needed for high accuracy 
determinations of ratios such as E,/E4. 

It is assumed that the changes in Y,, Yo, Y3. . . Yz 
which result from changes in e can be neglected. 
Following from this, we may consider the voltages 
applied to Y;, Ys, and Y; to be either E, or 2E, and 
define the values of Y;, Ys, and Y7 to be Y;, Y¢, and Y; 
at yoltage E4, and Y;(1+ ys), Y(1+ ye), and Y¥i(1+ yz) 
at voltage 2E 4. 

The balance equations for bridge 1 are 


SE l—q 
Y, Este, 1+88+q 
Y2 _ Ea—e2_ 1— qe 
Y, Eg+ees 1+ 8s + qe 

















FIGURE 2. Basic three-bridge-circuit. 


where Es = 1+ 5p. 
A 


Eliminating Y; and Y2, we obtain 


52 =— (qi + qz). (3) 
The balance equations for bridge 2 are 
Ys_ Ex tes _ ] + dp + qs 
5 ae 1+8¢ (4) 


Y,_Ext+es_1+ds+qs 
Y3 Ec 1+6c (5) 








where ce 1+ 8c. 


E Eliminating Y3 and Y,, we obtain 
A 


3 + gs + 26 
$0 +e, (6) 


where €; is a small correction term consisting of second 
and higher order terms in g. 


The balance equations for bridge 3 are 
Y;( 1+ Ys) oo E, —“€s = qs 
Y,+ Y; Ext E, T €s 2+ 68+ 8c 
Yel + ye) _ E,—€¢ ; 1 — qe 
Y3+Y; EsptEct+es 2+688+5c 








Yx1 + 7. E, —€7 Pa ] = 
Y;+ Ye Expt+Ect+e; 2+ 62+ 86¢ 





Eliminating Y;, Yg, and Y;, we obtain 





- 82+ 8ce+qs+qe+ 
y=- B 5. q: d6 en 


2 
where 


3 st+yYet yz 
pene 





and €» is a small correction term consisting of second 
and higher order terms in q and y. 
Solving eqs (3), (6), and (10), we obtain 


> 4qi + 4q2 — q3 — qa — 2qs — 2q6—2q7 _ €1— 2€2 
, 4 2 





(12) 


In order to separate real and imaginary parts, let 
yY¥=y'—-)j’ (13) 
= 4’ — a" 
so that 


f = 4g; + 4q2 — q3 — qs — 2q5 — 2q6 — 2q7 _ €1— 2€3 
4 2 





and 


e;—2e; (16) 


4 ° a 





y" _ 4qi+4q3 — g3 — gi — 245 — 2q6 — 2q7 


For the capacitors and accuracies described in this 
paper, y’ is equal to the proportional change in capaci- 
tance (AC/C), and y" is equal to the change in dis- 
sipation factor, corresponding to a specified increase 
in voltage. 

The values of ys, ys, and yz can be separated from 
the value of Y by measuring the differences, ys— ye. 
Ye—Yz. and yz—ys. Measurements of this type are 
described in section 4. Combining the differences, 
we obtain 


¥s=¥+(y¥s — Ye)/3 — (yz — ys)/3 (17) 


¥6=¥ + (ye — y7)/3 — (ys — ye)/3 (18) 


¥7=¥+ (yz — ¥5)/3 — (Ye — y7)/3. (19) 


3. Practical Three-Bridge-Circuit 


The method described in the previous section was 
used to determine the voltage dependence of three 100 
picofarad (pF) air capacitors, designated 100A, 100B, 
and 100C. This section describes the equipment and 
techniques used to determine the values of 


= _ Yi004 + Yi00B + Yi00c 
Yi100 = 3 





(20) 


corresponding to the following voltage changes: 12.5 
to 25, 25 to 50, 50 to 100, and 100 to 200 V, rms, at 
1592 Hz. 

If small capacitors are to be well defined, they must 
be shielded. The addition of shielding to the circuit 
of figure 2 is shown in figures 3a, 3b, and 3c. 
the admittance standards represented in figures 3a, 
3b, and 3c consists of a completely shielded three- 
terminal air capacitor connected to a pair of coaxial 
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FIGURE 3a. Bridge 1 of the practical three-bridge-circuit. 


FIGURE 3b. Bridge 2 of the practical three-bridge-circuit. 


FIGURE 3c. Bridge 3 of the practical three-bridge-circuit. 
cables. The defining terminals of each standard 
are located at one of the detector junctions and at one 
of the coaxial connectors designated a, 6, c i 
The required circuits are obtained by joining con- 
nectors a. 6b, c ...g to connectors, 4, B,C ...G 
in various configurations. When connector a is joined 
to connector A and b to B, the balance equation is 
given by equation (1), section 2. When a is joined to 
B and b to A, the balance equation is given by eq (2). 
Similarly, eqs (4) and (5) are applicable to figure 3b 
and eqs (7), (8), and (9) are applicable to figure 3c. 
For purposes of evaluating the effects of self and 
mutual impedances in the leads (not shown in the 
figures), it will initially be assumed that: (1) net cur- 





rents through the “coaxial chokes” [6], designated Ch 
in the figures, are negligibly small, (2) the voltages at 
connector A are, to sufficient accuracy, independent 
of whether A is joined to a or to 6, and similarly for the 
voltages at connectors B,C, D . . . G, (3) only negli- 
gible effects result from externally induced emfs and 
external loading, e.g., capacitance from the bridge 
shielding to its surroundings. 

Subject to the above assumptions, the two values of 
e corresponding to the two balance conditions for 
bridge 1 yield a measure of the voltage ratio, V4/V, 
where V’, is the voltage between the inner and outer 
terminals of connector A with e=0 and Vz is the 
voltage between the inner and outer terminals of con- 
nector B with e=0. Similarly, a measure of the volt- 
age ratios,V ¢/Vp and V/V» are obtained using bridges 
2 and 3, respectively. The value of Yioo is calculated 
from the measured values of V4/Vg. Vc/Vp, and VelV re. 
It is assumed that the relationship between voltages 
Vi, Ve, Ve V; is known. Ideally, the voltage 
drops in the leads which precede connectors A, B, 
C ...G should be zero, but if the voltage drops are 
small, eq (12) is sufficiently accurate when the rela- 
tionship between the voltages is as follows: Vs=— Vp, 
Ve=Ve=Ve, Ve=VatVe. Since voltage ratios 
rather than voltages are measured, other satisfactory 
relationships between the voltages can be described, 
e.g., V~ and Vp can each be changed if V ¢/Vp remains 
unchanged. 
Yi00 Which results from voltage drops in the leads will 
be small if the impedances in the leads which join 
various terminals of connectors A, B, C ...G are 
small or properly matched. These features are present 
in the junction box shown within broken lines in figure 


4. Parallel lines represent copper strips separated by 
thin strips of insulation. 

The error in the measurement of Yioo resulting from 
lead impedances within the junction box and resulting 
from contact impedances between connectors A, B, 


C...G and connectors a, b,c. 
approximately | x 10~°. 


. g was found to be 
Most of this error was cor- 


Junction oe 3 


FIGURE 4. Source of bridge voltages. 
Junction box (within broken lines) connected by coaxial cables to an “active direct 
reading ratio set” [7]. 


Thus, the error in the measurement of 


rected for by changing the currents at the terminals 
of connectors a, 6, c . g and observing the change 
in the measured value of Yioo. Since the current at 
the inner or outer terminal of connector a was nearly 
proportional to the internal load, Y; + Y4, and similarly 
for the currents at the terminals of b,c, d,. . . g, the 
following procedure was used: (1) Y4, Ys, Ye . . . Ye 
were adjusted so that each of the seven internal loads, 
Y:+ Ya, Yo+Yp. Yg3+YVe.. . ¥2+Yo, were increased 
by a known factor, (2) the change in the measured 
value of Yioo was used to calculate a third value of Yi00 
corresponding to the internal loads being zero. The 
above procedure amounts to an extrapolation of the 
effect of internal loads to zero. 


The internal loads were measured by disconnecting 
A,B,C ...G froma, b, ¢ . g, shorting the detec- 
tors, and measuring the two-terminal admittances be- 
tween the inner and outer terminals of connectors a, 
b,c... xg. Since the lead impedances within the 
junction box were small, low accuracy measurements 
were sufficient for purposes of extrapolating the 
internal loads to zero. The larger equivalent lead 
impedances in that portion of the circuit which pre- 
ceded the junction box (see fig. 4) required that more 
accurate measurements be made for purposes of ad- 
justing Yz, Yo. Yr. and Y¢g so that Y;+Y¥4=Yo2tYzx, 
Y3+VYco=YstYpo, and Y3+Ye=YetYr=Y7+ Ye. 
The accuracy of adjustment needed was determined 
by connecting additional loads between the inner and 
outer terminals of A, B, C . . . G, each in turn, and 
noting the changes in e required to rebalance the 
bridges. It was found that if the internal loads were 
adjusted within +1 X 10-* mho, the voltages at con- 
nector A were, to sufficient accuracy, independent of 
whether 4 was joined to a or to 6, and similarly for the 
voltages at B,C, D.. .G. 


The coaxial chokes used consisted of a number of 
turns of the coaxial cable threaded through a high 
permeability core. The resulting impedance to net 
current in the cable was about 400 0 at 1592 Hz. 
Since the impedance in the path, P; —A—a—Pz, of 
figure 3a was about 0.1 2, net current through the 
coaxial choke was of the order of a few parts in 104 
of the current in the inner or outer conductor of the 
cable [6]. The impedances in the paths, P;—A—a 
and P,—B—b, differed by less than 2 mQ. The 
maximum internal load was about 3000 pF. 

Thus far in the discussion the effects of externally 
induced emfs and external loading have been ne- 
glected. Although these effects were small in the 
practical circuit used, they were not negligible. Ap- 
propriate corrections were determined from additional 
measurements which consisted of recording the un- 
balance at the detectors when the admittance stand- 
ards of each bridge, in turn, were joined to the shorting 
connectors, A’, B’, C',. . . G', (see fig. 4) while the 
other standards were joined to A, B,C,. . .G. Since 
the accuracy of the corrections depended on the size 
of the unbalance at each detector, it was desirable for 
the effects of externally induced emfs and external 
loads to be small. The techniques used to reduce 
these effects will be considered next. 
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FIGURE 5. Bridge 2 redrawn to illustrate external loads and sources 


of induced emf. 


Bridge 2 has been redrawn in figure 5 for purposes 
of illustrating the different types of induced emf's 
and external loads. Since the shielding of bridge 2 
is not at ground potential, the largest external loads 
are associated with this bridge. External loads are 
represented by C,, C2, C3, and Cy in figure 5. Induced 
emfs are represented indirectly by showing areas in 
which a time-varying magnetic flux is assumed to exist. 
These areas are designated A,, Ay, A3, and A, in figure 
5. 

Most of the flux through A; existed in the region 
preceding the junction box and thus had no effect on 
the measured valué of Yioo. Flux through A» was re- 
duced by using coaxial cables of sufficient length so 
that the three-terminal capacitor could be positioned 
at a distance from sources of appreciable magnetic 
field. Flux through A; was reduced by using coaxial 
cables and magnetic shielding in that portion of the 
detector circuit which precedes voltage amplification. 
The effects of flux through A, were reduced by the 
use of coaxial chokes, e.g., if the high permeability 
core positioned between connector c and Py were 
removed from the circuit, the current through Y3 which 
results from flux through 4, would be increased by a 
factor of a few thousand. 

The effects of external loads were reduced by the 
use of coaxial chokes in conjunction with low-imped- 
ance leads which bypass the junction box (hereafter 
termed bypass leads). The lead connecting P3; to P, 
and the lead connecting P; to P; are examples of by- 
pass leads in figure 5. These leads, which have 
impedances of about 0.1 Q, supply most of the current 
to Cy and Cy. Thus, if coaxial chokes and bypass 
leads are properly positioned, voltage drops will be 
reduced in critical portions of the circuit such as the 
junction box leads and the cables of the admittance 
standards and detector circuit. 

The effects of external loads were also reduced 
by the use of an electrostatic shield which is repre- 
sented by a broken line in figure 6. This figure 
represents more closely the actual circuit used. The 
electrostatic shield, maintained at a potential of 
approximately E, with respect to ground, encloses 
those circuit elements which would otherwise have a 
large capacitance to the grounded surroundings. 
Current to the shield is supplied by the bypass lead 
connected to P3. 








FIGURE 6. Bridge 2 redrawn to show the difference transformer 


and additional shielding. 


WIRE (Coincident with axis) 
PRIMARY WINDING (100 Turns) 
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SHIELD NO. 2 
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SHIELD NO. 4 
SHIELD NO. 5 


CROSS-CAPACITANCE AREA 


FIGURE 7. Cross section through the axis of the difference trans- 


former. 


The transformer shown im figure 6 is termed a dif- 
ference transformer? and is used to change the poten- 
tial of the detector cables with respect to ground from 
approximately Eg, to zero. This was necessary in 
order that a phase-sensitive detector having a grounded 
power supply could be used. An adequate signal-to- 
noise ratio was maintained by using a battery-operated 
preamplifier having a voltage amplification of approxi- 
mately 10° (see fig. 6). Coaxial chokes assured that 
the bypass leads in figure 6 supplied most of the cur- 
rent to the large capacitance which existed between 
the shielding of the difference transformer. 

A cross-sectional view of the difference transformer 
is shown in figure 7. The secondary winding is a 


2 The basic ideas in the transformer are similar to those in a transformer built by A. M 
Thompson and termed by him a difference transformer 
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single turn consisting of shield No. 3 connected to a 
wire through the center of the core. With the excep- 
tion of shield No. 4 which is mu-metal, all shielding is 
brass. 

Most of the current resulting from the capacitance 
between shield No. 1 and shield No. 2 does not encircle 
the Supermalloy core. The current which does en- 
circle the core, if not compensated, would induce a 
voltage in the detector winding. The two capaci- 
tances which result in core excitation are associated 
with the corners of the shield gaps. With reference 
to the cross-capacitance area shown in figure 7, the 
two capacitances are (1) from the upper left to the 
lower right, (2) from the upper right to the lower left. 
It can be seen that the two currents encircle the core 
in opposite senses, and hence their effects can be 
made to cancel by proper adjustment of the gap 
geometry. The gap was adjusted experimentally by 
observing the detector voltage with the transformer 
input open circuited and with the bypass leads 
connected. 


As a result of the techniques described above, the 
unbalance at the detectors was small when the ad- 
mittance standards of each bridge, in turn, were 
joined to the shorting connectors, A’, B’, C’ . . . G’, 
while the other standards were joined to A, B, C 

G. Thus, accurate corrections were obtained 
for most of the effects of induced emf’s and external 
loads. Those effects which were not accounted for 
in the corrections will be considered next. 


Since currents were reduced when the standards 
were connected to the shorting terminals A’, B’, 
GC. G', induced emf’s resulting from currents 
in that portion of the circuit shown in figure 4 were 
not entirely accounted for in the above corrections. 
The use of coaxial cables, parallel-strip leads, and 
magnetic shielding assured that the resulting errors 
were negligible. It is interesting to note that if the 
errors had not been negligible, they would have been 
accounted for in the extrapolation of the internal 
loads to zero. The effect of induced emf’s resulting 
from currents within the standards was negligible 
since the phase of the induced emf’s was shifted by 
180° when the standards were permuted. 


Magnetic fields in the vicinity of the junction box 
may induce emf’s which are selective with respect 


to connectors B, C ... G, thus causing errors 
in the measurement of Yioo. Induced emf’s of this 
type were kept small by using closely spaced copper 
strips in the construction of the junction box (see fig. 4). 
The following procedure was used to assure that the 
remaining errors were negligible: (1) A small probe 
was used to measure the magnetic field in the vicin- 
ity of the junction box, including the region where 
changes in geometry occur due to changes in the 
connections of the standards to the junction box, 
(2) the current in a loop of wire was adjusted so as to 
increase the measured field by a factor of approxi- 
mately 100, (3) the changes in e required 1o rebalance 
the bridges were related to changes ir Yioo. 


TABLE 1. 


The external loads, C; and C2 in figure 5, represent 
the capacitances between the grounde d surreundings 
and those cables which are located between the co- 
axial chokes and connectors c and d. The circuit 
was designed so that Cy and Cz would be small and 
nearly equal. Measurements of C; and C, and of the 
lead impedances in the junction box indicated that 
the resulting error in the measurement of Yioo was 
small. This error is accounted for in item 1, table 1. 
Also accounted for in item 1, table 1 are errors re- 
sulting from capacitance between the copper strips 
of the junction box and errors resulting from imperfect 
balancing of the internal loads. 


Uncertainties in the measurement of Yi00 = Yi00 JV 100 
: 
Sources of error 


Systematic errors in| Systematic errors in 


y y. 
100 | 100 


Parts in 10” Parts in 10'° 
Loading 

Contact impedances 2 
Detector zero | 1 
Calibration of g 1 
Drift of admittances and transformer | 
ratio | I 





Sum 


Standard deviation of 5 


Table 1 lists only the significant errors which were 
not corrected for by extrapolating the internal loads 
to zero or by using the shorting terminals, A’, B’, 
C’ ... G', as described earlier. Item 2, table 1, 
accounts for small changes in ¥ contact impedances 
between a, b,c... g and B.4.. . G which 
were found to occur with Freotnorkic use. Item 3, 
table 1, accounts for possible systematic errors which 
result from observer bias in determining the null 
balance condition. 

Item 4, table 1, accounts for errors in the calibration 
of the instrument used to obtain E,. Ex. Ec. and e. 
The essential features of the instrument, termed an 
“active direct-reading ratio set,” are described in 

Le eee 
a i. oa 
were made by adjusting two sets of dials, one cor- 
responding to changes in gq’ and the other to changes 
in gq’. The maximum range of q’ or q” was + 10 ppm 
with a minimum adjustment of 2 < 10 

Item 5, table 1, accounts for errors which result 
from the admittances and transformer ratio having a 
nonconstant drift rate and from observer bias in timing 
the bridge balances which were made so as to obtain 
values of q in the following sequence: qs. q4, 43: Qi. Q2- 
Giz Qs+ G6. zs G6. 9s: Gis G2. Gi: Gs» G4. Gs. The above 
set of values could be obtained in about 12 minutes 
with E,=100 V, rms. The standard deviations of 
Fioo and Y{o9. Obtained from one such set, were about 
6xX10-'©. Slightly larger standard deviations were 
obtained at lower voltages. A sufficient number of 
measurements was made so that the final values of 
Vioo and ¥jo. had standard deviations of about 3 x 10-"". 


reference [7]. The adjustments of 
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4. Difference Measurements 


The difference in voltage dependence between two 
standards whose admittances are approximately equal 
can be determined from ordinary bridge measure- 
ments, e.g., simple substitution combined with a 
change in bridge voltage. The particular circuit 
used for this type of difference measurement is de- 
scribed in section 4.1. 

When the admittances of the two standards are 
radically different, additional measurements are re- 
quired. These measurement. ire described in section 


4.2. 


4.1. Difference Measurements Involving 
Approximately Equal Admittances 

The circuit shown in figure 3a is typical of the 
circuits used for difference measurements. The 
two general balance equations for this bridge have 
already been presented, namely eqs (1) and (2). For 
the case in which the voltage, E4, has two specific 
values, Ey, and E4y, four specific balance equations 
may be written: 
Yu 

Yo, 

Yor 

Yur 
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™ 

) 1H 


Bat el Me Se 
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Exiteu + OnL t+ qu 


E i — C2u 
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1+ 6en+ qu 


Ean T Cin 
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Exuten 1+ 8en+ qu 


Yu =Y¥u+ yi) 
and 


You Yo(1 + y2) 


then 


Vi — V2 = (Gen — Gor) — (qin — qit) + €3 (27) 
where €,is a small correction term consisting of second 
and higher order terms in q and y. 

The method described above was used to measure 
each of the following differences: yioo, — Yoon. Y100n 
—Yio0c. Yirwoc — Yiwoa. Yi0004 — Yio00n. Yi000R — Yi000c - 
and Yio00c — Yioo,. Where the subscripts refer to three 
100 pF air capacitors and three 1000 pF air capacitors 
designated 1OOA. 1LOOB, 100C, and 1OOOA, 1000B. 
1000C, respectively. Four values were obtained for 
each of the above differences corresponding to the 
follow ing voltage changes: 12.5 to 25, 25 to 50, 50 to 100. 
and 100 to 200 V, rms, at 1592 Hz. The discussion of 
measurement uncertainties which follows is appli- 
cable to each of these values. 


Since the differences, g2u— qe. and qin— qi, were 
small and of opposite sign in eq (27), most errors of 
the type described in section 3 were negligible. The 
only exceptions were errors resulting from drift of 
the admittances and transformer ratio. These errors 
were partially reduced by the measuring sequence 
which follows: iL, QiHs YiLs, Q2L, Q2H, 2. Rather than 
attempting to obtain a precise bridge balance at a 
specified time, each balance in the above sequence 
was made as quickly as possible. Nonuniform timing 
was almost inevitable since the detector sensitivity 
was better at high voltage. The resulting errors were 
kept small by injecting a small current to compensate 
the effect of capacitor drift. This was accomplished 
with an auxiliary circuit, termed a “drifter circuit.” 

If no adjustments of e are made, the relative changes 
in capacitance with time result in a time-variation in 
the amplitude of the current at the detector. The 
drifter circuit is used to inject a compensating current 
at the detector. The basic part of the circuit is a three- 
terminal, variable, air capacitor having a maximum 
capacitance of 10 pF. The adjusting shaft of the 
capacitor is connected to a constant-speed motor 
with speed reducer, thus providing a capacitance 
which varies linearly with time. Coaxial cables are 
used to connect the terminals of the capacitor to the 
detector junction and to an inductive voltage divider 
whose voltage is supplied by the “active direct-reading 
ratio-set”” used to obtain E,4, Eg, and e. 

The amplitude of the current injected at the de- 
tector has a first derivative with respect to time which 
can be adjusted by adjusting the inductive voltage 
divider and a second derivative which is zero. Since 
the drift rate of the admittance standards changed 
only slightly during the time needed for a set of detec- 
tor balances, the use of the drifter circuit resulted in 
appreciable reductions in the time rate of change of e 
required to maintain a detector balance. In addition, 
the deviations from constant drift rate, which can 
result in systematic errors, were easily measured. 

A sufficient number of measurements was made so 
that the standard deviation of each of the final values, 
Yiooa — Yioos> Yiooa — Yioon: Vieson ~ Yieeen? etc., was 
3X10-' or less. A systematic error of 2x 10-™ 
was assigned to each of the final values to account 
for incomplete compensation of the admittance and 
transformer ratio drift. 


4.2. Difference Measurements Involving 
Radically Different Admittances 


The circuit shown in figure 8 is used to measure the 
difference, 





pip ee _ Yet yot Yio Yst yer yz 


Ya— Yo = (28) 


3 3 


If connectors Q, R, and S are permuted cyclically 
with respect to connectors E, F, and G of the junction 
box (fig. 4), six different bridge balances can be made. 
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FIGURE 8. Circuit used for difference measurements with 


Y; Ys Y; # Ys Yo Yio. 


Three of the six balance equations have already been 
presented, eqs (7), (8), and (9). The other three 
equations are easily obtained by changing the sub- 
scripts of Y’, y, e, and q in eqs (7), (8), and (9). Solv- 
ing the six equations, we obtain 





+ €4 


Oe a 


< 


(29) 


where €, is a small correction term consisting of sec- 
ond and higher order terms in y and q. 

The method described above was used to measure 
the difference, 


_ Y10004 + Y1000B Le Y1000c 





1000 ie Y100 _— 


3 


__ Y100A 5 Y100B + Y100c 
3 





(30) 


Four different values were obtained for the above 
difference, corresponding to the following changes in 
voltage: 12.5 to 25, 25 to 50, 50 to 100, and 100 to 200 V, 
rms, at 1592 Hz. The discussion of measurement 
uncertainties which follows is applicable to each of 
these values. 

Corrections were obtained for the effects of induced 
emf’s and external loads (see section 3) by joining 
connectors Q, R, and S to the shorting connectors, 
E', F', and G’ (fig. 4). It was not necessary for the 
internal loads to be equal or to extrapolate the internal 
loads to zero. Errors resulting from admittance and 
transformer ratio drift were reduced by obtaining 
values of q in the following sequence: qs, gs, 47, qe, 
G5, 98, G9, Yio, 99, 98, 95, G6, W7, J6, Qs- A sufficient 
number of measurements was made so that the stand- 
ard deviation of each of the differences, ¥ioo9 -- Yio 
and Yioo0— Yioos Was reduced to about 5X10". In 
addition, an allowance of 5X 10-' was assigned to 


each of the differences to account for systematic 
errors which may be described as similar to items 3, 4, 
and 5 in table 1. 


5. Results 


The values in tables 2, 3, and 4 were calculated 
from the results of measurements described in sec- 
tions 3 and 4. Capacitors 100A, 100B, and 100C are 
of cylindrical construction, whereas capacitors 1000A, 
1000B, and 1000C are of parallel-plate construction. 
Most of the voltage dependence of capacitors 1000A, 
1000B, and 1000C is believed to be caused by elec- 
trode deflections [3]. Most of the voltage dependence 
of capacitors 100A, 100B, and 100C is believed to be 
associated with the electrode surfaces. 

The values in table 4 were obtained either imme- 
diately after assembly of the capacitors, one month 
after assembly, or eight months after assembly. The 
values in table 2 were obtained eight months after 
assembly. A comparison indicates that the values of 
Yiooa, Yioop, ANd Yiooc were probably smallest imme- 
diately after assembly. The changes in y are believed 
to have been caused by structural changes on the 
surfaces of the brass electrodes or by the migration 
of thin films onto the electrode surfaces from regions 
of the capacitor which were difficult to clean. A con- 
sideration of typical dimensions leads to the conclu- 
sion that sizeable effects could result from thin films, 
e.g., the ratio of molecular diameters to typical elec- 
trode separations is considerably larger than the 
desired accuracy of a few parts in 10°. 

Measurements, similar to those in table 4 and with 
comparable accuracies, were also made on capacitors 
1000A, 1000B, and 1000C. The largest change during 
the eight month period was 4X10°'. Capacitors 
1000A, 1000B, and 1000C had been in use for a num- 
ber of years prior to the measurements described. 
The electrodes were gold plated. These capacitors 
and capacitors 100A, 100B, and 100C were evacuated 
with an oil sealed mechanical pump trapped with 
dry ice. 


TABLE 2 | oltage depende nce of capacitors 100A. 100B. and 100¢ 


Values 
changes in) dissipation factor 
ix 10-* 


oly represent tional changes in capacitance (AC ( 


Values of 
SMandard deviation = 3 « 10 s 


ystematic 


Change in voltage 


Volts irm 
100 te 200 
te LOO 
25 to SO 


12.5 te 25 


Change in voltage 


Volts (rms Parts in 10" Parts in 10° Parts in 10° 
100 te 200 +0.) 04 | O.0 
30 te 100 0 + 2 2 
25 to 50 } 
12.5 to 25 i) 
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TABLE 3. Voltage dependence of capacitors LO000A, 1000B, and 


1000C 
Standard deviation = 6 * 


10°) Systematic error 





Change in voltage 


Volts (rms Parts in 10° 
100 te 200 +708 


0.4 
ite LOO + 18.2 +149 
25 4.3 

1.2 


Parts in 10° Parts in 10° 


te 50 +50 
+20 


Change in voltage 


Parts in 10° Parts in 10° 
‘ , + 3 


7) +2.2 
+ 1.1 + 1.6 
+03 | +O.8 


} Instabilities in the voltage dependence 
100A, 100B, and 100C 


of capacitors 


hange in voltage from 100 te 200 V. rms 
2*™ 10 1 


Standard deviation 
0 Systematic error 


Parts in 10° Parts in 10° 
+09 0.4 
) 


1.0 
Time “ A y" y 
5 100) 00) tow 1004 1004 


Vonths Parts in 10° Parts in 10° Parts in 10° 
0 0.2 +04 0.2 
l +1 | 0 | 
8 + 2 4 =. 





The values of y’ in tables 2 and 3 can be used to 
determine the functions which relate capacitance to 
voltage. If we assume the function 


Cy=CA(1+kvV™) (31) 


where Cy is the capacitance at voltage V and Co is the 
capacitance at zero voltage, then the values of y’ in 
each column of tables 2 and 3 will be in the ratio of 


23m; 22m:2™:1. Values of m were calculated from the 
measured values, Yio0a, Yioop, Yioooa, etc., and also 
from the measured differences, yioo, — Yi00p, Y10004 
— Y1000p, etc., which had smaller uncertainties of essen- 
tially different origin (see sec. 4.1). The results of 
the two types of calculations differed only slightly. 
The calculations yielded m ~ 3/2 for capacitors 100A, 
100B, and 100C and m~2 for capacitors 1000A, 
1000B, and 1000C. 

As a partial check on the overall measuring system, 
the value of ¥V1000 = (Y1000a + Y10008 + Y1000c)/3 from 100 
to 200 V was determined by two methods: (1) Indi- 
rectly by measuring the difference, Y:o00o— i100, as 
described in section 4.2, (2) directly by using capac- 
itors 1000A, 1000B, and 1000C in bridge 3 of the 
three-bridge-circuit (sec. 3). The results differed by 
less than 3 X 10-"°. 


Capacitors were chosen from more than two dozen 
capacitors on the basis of stability ascertained by 
difference measurements, similar to those described 
in section 4.1. The voltage dependence of capacitors 
100A, 100B, 100C, 1000A, 1000B, and 1000C was found 
to be considerably more stable than that of the other 
capacitors. The measurements were limited to 100 pF 
capacitors and 1000 pF capacitors. The use of the 
drifter circuit (sec. 4.1) eliminated the need for spe- 
cial temperature control of the capacitors. 


The earliest measurements involved three commer- 
cial 100 pF air capacitors of parallel plate design. 
The values obtained from difference measurements 
were observed to depend not only on time but also on 
the voltages which had previously been applied to 
the capacitors. In addition, a sudden, nonreversible 
change in capacitance of 14 ppm was observed in one 
of the capacitors when the applied voltage was in- 
creased beyond 50 V rms. The cause of the unusual 
behavior was found to be metallic whiskers on the 
tin-plated housing of the capacitor. 

Difference measurements were also made on a set 
of ten 100 pF capacitors of cylindrical design. The 
construction of these capacitors is described in ref- 
erence [2]. The values obtained from the difference 
measurements, corresponding to a change in voltage 
from 80 to 160 volts, were observed to change by as 
much as 3 X 10-8 over a period of afew days. Changes 
as large as 3 X 10-® were also observed when the tem- 
perature of the capacitors was increased by approxi- 
mately 10°C. 

A different phenomenon was observed when dif- 
ference measurements were made on a set of three 
100 pF capacitors of cylindrical design which were 
evacuated during measurement. The time rate of 
change of capacitance was found to be a function of 
voltage. The phenomenon is believed to be related 
to outgassing of the epoxy which was used to insulate 
and support the electrodes. Measurements of the 
time rate of change of capacitance as a function of 
voltage were made as follows: (1) The bridge was bal- 
anced and the drifter circuit was adjusted so as to 
obtain a steady bridge balance, (2) the detector was 
shorted, and the voltage applied to one capacitor was 
reduced to zero for 30 sec, (3) the original conditions 
were restored and the change in bridge balance was 
recorded. With 160 V originally applied to the ca- 
pacitors, the changes in bridge balance were equivalent 
to proportional changes in capacitance of about 1 x10-*. 

Difference measurements, corresponding to a 
change in voltage from 100 to 200 V, were made on 
three commercial 1000 pF capacitors of parallel! plate 
construction. The changes in the values obtained 
from the difference measurements varied from a few 
parts in 10° to 2X 10-8 over the period of one week. 


The basic design of capacitors 100A, 100B, and 
100C is shown to approximate scale in figure 9. Crit- 
ical radii differ by approximately 3.2 mm. Line and 
detector electrodes are separated from the ground 
electrodes by four sets of three glass spacers. The 
force of the compressed springs is about 45 N. 
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FIGURE 9. Cross section through the axis of capacitor 100A, 100B, 


or 100C. 


It was found that after the capacitor was assembled 
and connected into a measuring circuit, very small 
radial displacements of the line electrode could be 
obtained by sharply tapping its base radially. This 
procedure was used to adjust the capacitance very 
close to a minimum, thus assuring that the radial 
components of the electrical forces were small. The 


design of the capacitor is such that it can easily be 
disassembled for cleaning or for experimental pur- 


poses. A grounded housing (not shown in the figure) 


is used to evacuate the capacitor. 


6. Conclusion 


The basic purpose of the present work has been to 
develop an accurate system for measuring both the 
relative and absolute voltage dependence of admit- 
tance standards. The causes of voltage dependence 
have been pursued only to the extent necessary to 
reduce instabilities. Additional measurements will 
be required to evaluate the instabilities which remain. 


Since the mechanisms which cause instabilities in 
the voltage dependence of capacitors may also result 
in general instability, the present work has pointed 
out certain problem areas which should be considered 
in the construction of precision standards of capaci- 
tance. In addition, the present work has served as a 
proving ground for bridge measurements and voltage 
ratio measurements at the higher accuracy levels. 


The author acknowledges the many suggestions of 
R. D. Cutkosky and the assistance rendered by D. N. 
Homan and Lai H. Lee in obtaining the numerical 
results of this paper. 
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Single-Crystal X-ray Diffraction at High Pressures 
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Construction and operation of a high-pressure single-crystal x-ray diffraction precession camera 


is described. 


single crystal studies of ice VI and ice VII. 


lhe instrument appears applicable to single crystal studies to pressures as high as at 
least 30 kilobars and temperatures of at least 250 °C. 


Application of the instrument is illustrated by 


Unit cell parameters are given for these two forms of ice. 


Further modifications and applications are discussed. 


Key Words: High pressure, single crystal, x-ray diffraction, ice VI, ice VII, water, calcite. 


1. Introduction 


X-ray diffraction of materials under high pressures 
presents obvious experimental problems arising from 
the necessity of using a vessel which is transparent 
to x rays, but sufficiently strong to contain the pres- 
sures desired. Jacobs [/]|' fabricated a steel vessel 
equipped with beryllium windows and used helium as 
the pressure medium in his pioneering studies, but 
little additional progress was made until Lawson [2] 
conceived and used the diamond pressure vessel. 
Imposing progress has been achieved in recent years, 
and the problem of x-ray diffraction appears to be 
solved in principle up to pressures as high as several 
hundred kilobars with accompanying temperatures as 
high as 1000 °C [3, 4, 5]. In practice, however, the 
diffraction data are limited in diffraction angle and 
further to those reflections sufficiently strong to over- 
ride the unavoidable background arising from scatter- 
ing by the pressure cell itself. All significant work, to 


date, involves powder x-ray diffraction and, in view of 


the present limitations of the methods, there is little 
promise that unequivocal structures of any high pres- 
sure phases, except those of the utmost simplicity, 
can result from such methods. 

Some high-pressure polymorphs can be stabilized 
at 1 bar by quenching to low temperatures while still 
under pressure. Both powder and single-crystal x-ray 
diffraction studies [6-1/2] have been made on several 
such materials in the metastable state. As it appears 
that only a limited number of presently known high 
pressure polymorphs are quenchable, there is a need 
for producing and studying single crystals of high 
pressure phases in their ranges of stability by means 
of x rays. This paper is a preliminary report of tech- 
niques developed for single-crystal x-ray diffraction 
at high pressures. 


' Italicized figures in brackets indicate the literature references at the end of this paper 


2. Apparatus and Experimental Method 


2.1. Apparatus 


The instrument used for single crystal studies is 
the diamond-anvil high-pressure cell developed at 
NBS. This cell has been used here and in other 
laboratories for spectroscopy, microscopy, and pow- 
der diffraction studies [/3, 14, 15]. A technique for 
mounting a metal gasket between the diamond anvil 
faces was developed by Van Valkenburg [/6] which 
permits the application of pressure to liquids. Using 
this technique it is possible to grow single crystals 
from liquids under pressure or to subject a single 
crystal immersed in the liquid to purely hydrostatic 
pressure. The modifications of the cell and the ordi- 
nary precession camera to permit high-pressure single 
crystal studies are as follows: 

The original cell designed for spectroscopy [13] 
has been improved and is shown in figure 1. The 
instrument has opposed diamond anvils, A, to gen- 
erate pressure. The anvils are fabricated from 
gem quality diamonds by grinding small flat surfaces 
approximately 0.6 mm and 1.2 mm in diameter parallel 
to the table faces. The diamonds are seated on their 
table faces in recesses in metal plates, B, that are 
held to cylinders, C, by means of two screws 180° 
apart. Two other screws at an angle of 90° to the 
holding screws permit alinement of the diamond 
anvil surfaces. The cylinders fit snugly into bear- 
ing, D, which has a shoulder at one end to support 
the thrust applied to the cylinder. The thrust which 
forces the anvils together is supplied by plate, E, by 
means of spring, F, acting through top plate, G, and 
lever arms, H. The diamond support plates, B, have 
always been made of stainless steel, but, for reasons 
to be discussed later, have been replaced by high 
purity beryllium plates. In both cases, a central 
hole is drilled in each plate to permit microscopic 
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FIGURE 1. 


observation of the specimen. To permit the entry 
or exit of x rays at higher angles, it was necessary to 
machine the stainless steel plates out at an angle 
of about 45°. Beryllium metal is relatively trans- 
parent to the x rays used here, and angular relief is 
not required. Due to spatial limitations the usual 
knurled knob on the screw, was replaced with a disk, 
I, which could be turned by means of a removable 
crank. Elevated temperatures were obtained by 
means of a small furnace, J. The furnace was made 
by winding Nichrome wire on a ceramic sleeve which 
fit snugly over the outside of D and inside the block, 
K. Temperature was recorded and controlled by 
means of a thermocouple inserted through a small 
hole in one of the cylinders, C, and plate, B, so as to be 
in direct contact with the metal gasket placed between 
the diamonds. An automatic recording-controlling 
unit was used to monitor temperatures. 

It was necessary to modify the commercial preces- 
sion camera to sustain the weight of the diamond cell. 
The existing crystal bearing assembly was eliminated 
by removing the mounting pillar and replaced by a 
sturdier unit constructed of aluminum to minimize 
its weight. An aluminum alloy shaft 1 in. in diameter 
was coupled to the film precession shaft as in the 
original camera. A large graduated flange was fitted 
over the end of the shaft and locked in place by means 
of set screws. This flange replaced the original “dial.” 
A modified goniometer head was fastened to the other 
end of the flange. This provided only a single arc cor- 
rection in an approximate horizontal plane. The two 
sections of this arc were also locked together with 
set screws. The customary vertical arc could not be 
accommodated in the available space and was elimi- 
nated. A stout horizontal bar clamped to the arc 
correcting unit was used to mount the diamond cell. 
This bar could be moved on a track with positioning 
screws to permit adjustment in the direction parallel 
to the x-ray beam. Vertical motion of the cell was 


Construction details of diamond cell (not to scale). 


achieved by means of a plate fastened to the diamond 
cell with slotted screw holes and sliding on a vertical 
track. The cell was clamped to the horizontal bar by 
screws into this plate. The screw holes through the 
horizontal bar were also slotted to permit the cell to 
slide on a track on the bar. This assembly permitted 
translation along three mutually perpendicular axes, 
and rotation about two mutually perpendicular axes. 
Because of interference between the cell and the x-ray 
tube or the collimator support, and limitations imposed 
by the metal of the cell cutting off the x-ray beam, 
maximum angular corrections of the order of + 15° 
were permitted. 


On assembly it was necessary to aline the new pre- 
cession axis and the diamond cell to this precession 


axis. Most of the operations required are obvious 
and only the operations required to locate the speci- 
men on the precession axis will be described. It 
was found that vertical and lateral adjustments were 
most easily initiated by visual alinement of the gasket 
hole with the collimator with Z=0. The gasket hole 
can be observed through the diamond. Refinements 
in these adjustments were made by maximizing the 
intensity in the spot produced by the x rays on a 
fluorescent screen placed behind the diamond cell. 
The intensity was first maximized at 7@=0 and final 
adjustments were made at higher values of &. The 
adjustment of the cell in the direction parallel to the 
x-ray beam was achieved by eliminating doubling of 
diffracted spots. 

In use, it was expedient to remove the cell frequently 
for microscopic observation, and it was found to be 
essential to provide a means of ensuring that the cell 
was replaced in the same position to eliminate the 
need for realinement. This was achieved by means 
of the track on the horizontal bar which fixed the 
vertical and one lateral position and by means of 
metal stops on the bar which fixed the other lateral 
position. With these guides the cell could be removed 
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FIGURE 2. 


and replaced in a few seconds with no alinement ad- 
justments required. The mounted cell is shown in 
figure 2. 

For the diamond cell and camera used here it was 
found that with dial and arc readings near zero, the 
maximum useful value of & was about 20° with steel 
plates and 30° with beryllium plates. As the dial 
and arc values increased the maximum useful f value 
decreased. With the beryllium plates, however, val- 
ues of 20° were easily obtained at the maximum arc 
and dial settings available. 


2.2. Sample Preparation 


The initial step in sample preparation involves fab- 


ricating the metal gasket. A thin sheet of metal 
(0.010 in. X 0.25 in. X 0.25 in.) is subjected to an initial 
compression between the diamond anvils to produce a 
well-defined indentation of the smaller of the two dia- 
mond surfaces. The pressure required for this 
operation must be adjusted to the properties of the 
gasket material and may vary from 5 kbar-60 kbar 
for different metals. Under a low power microscope 
a small punch mark is made in the center of the in- 
dentation and a small hole (No. 80 drill) is drilled 
through the gasket using the punch mark for centering. 
For pressure work above about 15 kbar it is essential 
that this hole be centered in the indentation to mini- 
mize the failure of the gasket; at lower pressures the 
centering requirement is less stringent. The burrs 
formed in drilling are removed with a small jeweller’s 
file and the gasket is placed over the small diamond 
using the indentation as a guide. In this operation the 
cylinder containing the small diamond should be in 
position in the diamond cell with the diamond anvil 
face uppermost. If a liquid is to be studied a small 
drop of liquid is placed on top of the gasket in contact 
with the hole, the cylinder with the larger diamond 


Assembled single crystal precession unit. 


inserted and the cell carefully assembled to avoid 
displacing the gasket. If a solid is to be studied, the 
selected single crystal may be placed in the hole and 
oriented on the diamond surface with a small needle 
before the hydraulic liquid is applied or the liquid 
may be applied first and the crystal inserted through 
the liquid. In either case difficulties are to be ex- 
pected in maintaining proper alinement of the crystal 
during the assembly of the cell and the initial applica- 
tion of pressure. One technique used here involves 
cementing the gasket and the crystal in position by 
means of thin Canada balsam before assembly. How- 
ever, cements are considered to be undesirable be- 
cause they interfere with microscopic observations 
and may prove to be a source of pressure gradients. 

Liquids which have been used successfully are 
water, glycerine, methanol, and Varsol. Both water 
and Varsol crystallize at moderate pressures, but 
glycerine and methanol or glycerine-water mixtures 
apparently do not crystallize to about 30 kbar. Glyc- 
erine proves to be sufficiently viscous to prevent ap- 
preciable motion of a crystal during the assembly 
process. Of these liquids, water, glycerine, and meth- 
anol are relatively incompressible, a matter of some 
importance because the crystal is crushed if it comes 
in contact with both diamond anvil surfaces simul- 
taneously. 

Some iron and nickel alloys serve to make satis- 
factory gaskets provided pressures above about 15 
to 20 kbar are not required. Brass and aluminum 
alloys are not useful to these pressures but might 
prove useful at lower pressures. Above about 20 
kbar problems involving gasket blowout are ex- 
perienced, presumably because the strength of the 
metal and the diamond-metal friction is exceeded 
by the pressure exerted by the liquid. Two solutions 
are indicated, one to increase the strength of the metal 
gasket, the other to increase the surface area of the 
diamonds. Both techniques appear to be effective. 
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When molybdenum metal gaskets were used no 
blowouts were experienced. Molybdenum appears 
to be more difficult to handle mechanically but is 
very ductile under pressure (see Bridgman [/7]). 
When the diameter of the smaller diamond was in- 
creased from 0.6 mm to 0.8 mm no blowouts were 
experienced using a nickel alloy gasket. The opti- 
mum gasket thickness for the diamonds used here 
appears to be about 0.010 in. Thicker gaskets have 
been used but do not perform satisfactorily, pre- 
sumably because of the nonuniform shape of the 
smaller diamond. This diamond penetrates the 
gasket and appears to produce nonuniform stresses 
in the metal resulting in distortion and frequent blow- 
out of the gasket hole. Gaskets 0.005 in. in thickness 
perform satisfactorily but provide a smaller available 
volume. 

Initally, experiments were conducted using Mo Ka 
radiation. Under these conditions orientation photo- 
graphs were obtained in Y4-hr with ordinary x-ray 
film. With this radiation, however, it was found that 
cone axis, and upper level photographs were not 
satisfactory because the diffracted beams at higher 
angles were blocked by the cylinders or the instrument 
itself. A marked improvement in quantity of data 
obtained has resulted from using Ag Ka radiation 
which has a shorter wavelength (0.5609 A versus 
0.7107 A). With Ag radiation useful cone axis and 
upper level photographs as high as the second level 
have been obtained with Fd* ~ 4. Using Ag radiation, 
however, some sacrifice in exposure time is required. 
By using fast film (such as Kodak No Screen or DuPont 
Cronex II) in conjunction with fluorescent screens 


(Cronex CB 2, Ilford, or Radelin TF)? placed behind 
the film, orientation photographs can be made in Y/a-hr. 
Exposures are increased to at least 1 hr without the 
fluorescent screen. 

It should be noted that the problem of collimation is 
simplified by use of a metal gasket that is opaque to 


the x rays. The hole in the gasket functions as a 
second pinhole and the external collimator need only 
contain a single pinhole. Use of the gasket as part of 
the collimating system produces diffraction rings from 
the metal in the resulting pattern. These rings have 
not caused any difficulties in interpreting the films. 
Extensive collimation of the beam before it is incident 
on the specimen and centering of the beam in the 
gasket hole can be used to eliminate diffraction rings 
from the metal. However, this technique results in 
a marked increase in exposure time which does not 
appear to be justified by the results. 


3. Results 


Initial studies using this instrument were made 
on materials which were crystalline solids at 1 bar and 
had known transitions at elevated pressures. The 
experimental technique used was to orient a crystal- 


Certam « 
adequately specify the experimental procedure In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards. nor does it imply 


that the material or equipment identified is necessarily the best available for the purpose 


ommercial materials and equipment are identified in this paper in order to 


lographic plane of the low pressure form using the 
precession camera. The pressure cell was then 
removed from the camera and mounted on a polarizing 
microscope. The pressure was increased slowly 
until the transition was observed by noting an abrupt 
change in birefringence or dimensions or a fracturing 
of the crystal. If the crystal remained intact in the 
high pressure form, the cell was returned to the cam- 
era for precession photographs. Under the condi- 
tions used to date, only one material of those studied 
appears to remain in single crystal form through the 
transition. The transition CaCQs (calcite) ~CaCOs II, 
[18] occurs with a marked change in birefringence 
and the crystal frequently remains intact. The transi- 
tion CaCOsll —CaCOslII always results in a_ poly- 
crystalline specimen. Transitions in KNOs; [19] and 
the alkali halides [20] invariably produced _poly- 
crystalline products and no transition could be de- 
tected in KIO, [2/]. 

In only one experiment with calcite was it possible 
to obtain a photograph of a densely populated plane 
of the CaCOslIl form. Zero level and cone axis photo- 
graphs could be interpreted on the basis of a cell with 
dimensions a=4.97 A, b=7.72 A, and c=13.25 A 
with the a and 6 axes orthogonal. No information on 
any other angles or on upper levels could be obtained. 
Characteristic absences h+k=2n+1 were observed 
on the zero level photograph. Davis [22] has reported 
powder diffraction data on CaCOslll and noted that 
CaCO Il gave diffraction data similar to ordinary 
calcite except that the 113 reflection was missing in 
CaCO,II. With the parameters found here and 
assuming an orthorhombic cell, it is possible to index 
almost all the lines reported for CaCQOs;I (calcite), 
CaCO,II, and CaCO, within the probable errors 
of the data. Further study is necessary before this 
situation can be cleared up. However, calcite 
crystals present problems in orientation because of 
their cleavage characteristics and it was decided to 
defer further studies until an improved instrument 
becomes available. 

From the limited studies conducted here on solid- 
solid transitions at room temperature in single crystals, 
it appears that rentention of the single crystal form 
through a high pressure transition will not occur as 
frequently as originally anticipated. Other techniques 
for production of high pressure polymorphs in single 
crystal form are available, such as crystallization from 
solution or from a melt under conditions of stability 
of the desired polymorph. Both techniques are pos- 
sible with the present instrument which has an upper 
temperature limit of operation above 250°C. Neither 
method has been attempted because of the great ease 
with which single crystals may be grown from ma- 
terials normally liquid at 1 bar and room temperature 
and current interest in the structure of such crystals. 


3.1. Ice 


Bridgman [23, 24] has shown that six forms of ice 
are formed .under pressure. Of these six forms only 
two, ice VI and ice VII, are formed from water at or 
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above room temperature. Ice VI crystallizes from 
water at room temperature and approximately 9 kbar 
and ice VII forms at pressures above about 22 kbar 
and 82 °C—the VI-VII-L triple point. Both ice VI 
and ice VII can be grown as single crystals in the dia- 
mond cell, ice VI growing readily but ice VII with 
considerable difficulty. 


3.1.1. Ice VI 


Powder diffraction data for ice VI have been re- 
ported by Bertie, Calvert, and Whalley [8] at 90 °K 
and | bar and by Kamb and Davis [25] at —50 °C and 
8 kbar. Both reports agree reasonably well on the 
powder lines and relative intensities when allowance 
is made for the different conditions of thie two measure- 
ments. From single crystal studies here a prelim- 
inary report has been given for a unit cell which can 
be used to index the observed powder diffraction 
data [26]. 

Typical precession photographs of a single crystal 
of ice VI are shown in figures 3-7, with figures 3, 4, 
and 5 being made with steel plates on the cylinders 
and figures 6 and 7 with beryllium plates. A typical 
orientation photograph using unfiltered Ag radiation 
with a fluorescent screen is shown in figure 3. The 
exposure time is Y2-hr with @=10°. The intense 
streaks and loops are contributed by the diamonds, 
the diffraction rings by the metal gasket, and the nar- 
row streaks by ice VI. A zero level photograph using 
Pd filtered Ag radiation is shown in figure 4. This 
is a 72 hr exposure with B= 16° and a layer level screen 
11.5 mm in diameter at a spacing of 40 mm. It is 
apparent that the filter and the layer level screen 
eliminate most of the diamond interference. It will 
be observed that diffraction spots in the 4th quadrant 
of this photograph are weaker than those in the 2d 
quadrant as a result of shadowing by the instrument 
itself. This occurred because of the angular correc- 
tion required to orient the crystallographic plane. 
The first level spots are diffracted at higher angles 
than those on the zero level and the shadowing by the 
instrument makes it much more difficult to photograph 
these levels as shown by figure 5. This is the Ist 
level corresponding to figure 4 and is a 72 hr exposure 
with @=10°. There are several rows of diffraction 
spots in this photograph concentrated principally 
in the 2d quadrant. The third quadrant which was 
slightly shadowed in the zero level photograph is now 
completely empty of diffraction spots. The shadowing 
restricts both the diffraction angle and the precession 
angle and for this reason it was decided to replace 
the steel plates, B, of figure 1 with beryllium plates. 
The effects of the substitution may be ascertained 
by comparing figures 4 and 5 with figures 6 and 7 which 
made use of the beryllium plates. The latter two 
photographs are 48 hr exposures with &#=22° and —& 
=10° respectively. A higher background is noted 
in the latter two photographs presumably arising from 
incoherent scattering by the approximately 1/2-in. 
of polycrystalline Be traversed by the x-ray beam. 


FIGURE 3. Typical orientation photograph (unfiltered Ag radiation 
using fluorescent screen; Z=10°; 1/2 hr exposure. 


FIGURE 4. Ice VI—hkO photograph with steel plates (Pd filtered 


Ag radiation; {= 16°; 72 hr exposure). 


However, no appreciable shadowing is apparent on 


either photograph. With the beryllium plates shadow- 
ing becomes important on the 2d level of this crystal 
but even on this level several reflections were obtained 
whereas none were obtained with the steel plates. It 
is of interest to note the threefold diamond axis in 
figure 6. 

Final answers on the unit cell and symmetry for 
ice VI have not been obtained to date. The unit cell 
dimensions at room_temperature and approximately 
9 kbar are a=8.38 A, b=6.17 A, c=8.71 A with un- 
certainties of 0.05 A in all cases. It is assumed 
that these axes are orthogonal. From the AsO films 
systematic absences are hkO, h=2n+ 1, and k=2n+1. 
From hkl, systematic absences are h+/=2n+1. 
From upper level photographs and upper level loops 
seen on zero level photographs it appears that a is 
a true twofold axis but 6 is not. To give an unequivo- 
cal space group requires more facility of orienting 
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FIGURE 5. Ice VI—hk1 photograph with steel plates (Pd filtered 


Ag radiation; {= 10°; 72 hr exposure). 


FIGURE 6. Ice VI—hkO photograph using beryllium plates (Pd 


filtered Ag radiation; fi= 22°; 48 hr exposure). 


the crystal and less difficulty with shadow produced 


by the cell. A new camera now under construction 
is expected to facilitate solution of this problem. 


3.1.2. Ice VII 


Bertie, Calvert, and Whalley [8] and Kamb and 
Davis [25] have reported powder diffraction data for 
ice VII. The former authors worked at 90 °K and 
1 bar and the latter at —50 °C and 25 kbar. Bertie, 


FIGURE 7. Ice VI—hkl photograph using beryllium plates (Pd 


filtered Ag radiation; ~7= 10°; 48 hr exposure). 


Calvert, and Whalley reported more reflections 
than Kamb and Davis, but noted that their strong 
reflections could be indexed in the cubic system. 
Kamb and Davis concluded that all their reflections 
were consistent with a body-centered cubic cell with 
a=3.30 A, and they have proposed this structure 
for ice VII. Single crystal studies of ice VII in this 
laboratory arrived at identical conclusions at about 
the same time. 

Single crystals of ice VII were grown in the diamond 
cell at approximately 90 °C and 25 kbar. The crystal 
grew as an anhedral spherical globule and was cooled 
to room temperature for study. It was found that very 
slow cooling was essential to prevent fracturing the 
crystal. Three different orientations of ice VII have 
been studied showing reciprocal lattice, planes per- 
pendicular to the [100] and the [111] axes, and a plane 
containing the [110] and [121] axes. This latter plane 
was derived from the orientation perpendicular to the 
[111] axis by angular corrections consistent with the 
body-centered cell proposed by Kamb and Davis [25]. 
There are few reflections observed and the intensities 
decrease very rapidly with increasing order as noted by 
Kamb and Davis [25]. No weak reflections indicative 
of lower symmetry or superlattice have been observed, 
although it should be noted that the small dimensions of 
the unit cell precluded upper level photographs. From 
the present data, it is concluded that the cell is cubic 
with a=3.40+0.05 A. Systematic absences h+k 
+ 1=2n+1 indicate a body-centered cell. 

Bridgman [27] observed that water at 25 °C under- 
goes a change in volume of 0.359 cm‘/g between 1 
bar and 25 kbar. Therefore, the measured density 
of ice VII at 25 kbar and 25 °C is 1.56 g/em*. From 
the dimensions of the unit cell obtained here contain- 
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ing two molecules of water the calculated density 
is 1.51 g/em*. The single crystal data, therefore, 
corroborate the body-centered cubic structure pro- 
posed for ice VII by Kamb and Davis [25]. 


4. Conclusions 


The present studies have shown that successful 
single-crystal x-ray diffraction work at high pressures 
is possible to pressures of at least 30 kbar. Success- 
ful structure analyses and space group determinations 
require more data and greater precision and facility 
of operation than the present precession camera can 
provide. A new goniometer head providing all the 
usual angular and translational motions and designed 
to support the weight of the diamond cell is being 
constructed. In addition, a new diamond cell con- 
structed entirely of high purity beryllium metal is 
being fabricated. 


Note added in proof: In relation to the structure of ice 
VI (see section 3.1.1) Kamb has recently published 
(Nature 150 205, (1965)), x-ray diffraction data on ice 
VI obtained at 1 bar and—175 °C. He concluded that 
the cell was tetragonal with dimension a=6.27 A and 
c=5.79 A. This cell differs from the preliminary cell 
presented here although the two cells are related. 
However, there is sufficient difference that the cell 
proposed by Kamb cannot index both the 3.4 A and 
3.6 A lines reported in the powder pattern by Bertie, 
Calvert, and Whalley [8] which the present cell can do. 
We have recently obtained an improved cell and orient- 
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ing device. The final answer to the cell of ice V1 under 
equilibrium conditions has not been obtained, but, by 
reorienting about our b axis we have obtained a zero 
level pattern which shows a four-fold axis of 3.11 A 
which is compatible within the limits of error with 
Kamb’s 6.27 A (2 X 3.11 A=6.22 A). Further studies 
are in progress. 


5. References 


. B. Jacobs, Phys. Rev. 54, 325, 468 (1938). 
A. W. Lawson and T. Y. Tang, Rev. Sci. Instr. 21, 815 (1950). 
J. C. Jamieson and A. W. Lawson, J. Appl. Phys. 33, 776 (1962). 
J. D. Barnett and H. T. Hall, Rev. Sci. Instr. 35, 175 (1964). 
». A. Perez-Albuerne, K. F. Forsgren, and H. G. Drickamer. 
Rev. Sci. Instr. 35, 29 (1964). 
R. L. McFarlan, J. Chem. Phys. 4, 60, 253 (1936). 
J. E. Bertie, L. D. Calvert, and E. Whalley, J. Chem. Phys. 38, 
840 (1963). 
[8] J. E. Bertie, L. D. Calvert, and E. Whalley, Canad. J. Chem. 42, 
1373 (1964). 
9| R. H. Wentorf and J. S. Kasper, Sci. 139, 338 (1963). 
0| F. P. Bundy and J. S. Kasper, Sci. 139. 340 (1963). 
1] C. Bates, F. Dachille, and R. Roy, Sci. 147, 860 (1965). 
2] B. Kamb, Acta Cryst. 17, 1437 (1964). 
3] C. E. Weir, E. R. Lippincott, A. Van Valkenburg, and E. N. 
Bunting, J. Res. NBS 63A (Phys. and Chem.) No. 1,55 (1959). 
[14] A. Van Valkenburg, Rev. Sci. Instr. 33, 1462 (1962). 
[15] G. J. Piermarini and C. E. Weir, J. Res. NBS 66A (Phys. and 
Chem.) No. 4, 325 (1962). 
[16] A. Van Valkenburg, unpublished. 
[17] P. W. Bridgman, Large Plastic Flow and Fracture (McGraw- 
Hill Book Co., New York, N.Y. 1952). 
[18] P. W. Bridgman, Am. J. Sci. 237, 7 (1939). 
. W. Bridgman, Pro. Am. Acad. Arts & Sci. 51, 581 (1916). 
. W. Bridgman, Pro. Am. Acad. Arts & Sci. 74, 21 (1940). 
. W. Bridgman, Pro. Am. Acad. Arts & Sci. 72, 45 (1937). 
. L. Davis, Sci. 145, 489 (1964). 
. W. Bridgman, Pro. Am. Acad. Arts & Sci. 47, 441 (1911). 
. W. Bridgman, J. Chem. Phys. 5, 964 (1937). 
. Kamb and B. L. Davis, Pro. Nat. Acad. Sci. 52, 1433 (1964). 
5. Block, C. Weir, and G. J. Piermarini, Sci. 148, 947 (1965). 
. W. Bridgman, Pro. Am. Acad. Arts & Sci. 74, 399 (1942). 


| 
| 
1 
l 
1 


| 
| 
[ 
[ 








JOURNAL OF RESEARCH of the National Bureau of Standards—C. Engineering and Instrumentation 
Vol. 69C, No. 4, October-December 1965 


The Sondheimer-Wilson-Kohler Formula in 
Platinum Resistance Thermometry 


R. J. Corruccini 


Institute for Materials Research, National Bureau of Standards, Boulder, Colo. 


(June 7, 


1965) 


An attempt was made to determine the ideal resistivity versus temperature function for platinum 
by a parameter-variation method using the Sondheimer-Wilson-Kohler formula for the deviation from 


Matthiessen’s rule. 


At the same time a test of the applicability of the SWK formula for practical 


thermometry was made by examining the accuracy with which it could fit the resistance-temperature 
functions of selected thermometers, the characteristics of which were restricted in a manner believed 


to favor their accurate representation. 


No ideal resistivity function could be found such that most 


of the thermometers could be accurately represented. 


Key Words: Cryogenic, electrical resistivity, low temperature, platinum, resistance thermom- 
etry, temperature measurement, thermometry. 


1. Introduction 


The resistivity of a metal can be expressed by the 
formula 


pr= pot pir+ Ar (1) 


where po is the imperfection-dependent resistivity at 
T=0, pir is the temperature-dependent “ideal” resis- 
tivity of the hypothetical perfect lattice, and Ar is a 
small term resulting from the interaction of defect and 
lattice scattering of the electrons. Sondheimer and 
Wilson [1]' showed on the basis of a two conduction- 
band model that 


PiTPo 
Ap=—_ ee ‘ 
; apir + bpo (2) 


and Kohler [2] obtained the same formula on more 
general grounds. In either treatment, the coefficients 
a and b are positive quantities and may be functions 
of temperature, but their forms have not been estab- 
lished from theory. 

Equation (2), with a and 6 taken to be constants, 
was first applied to platinum resistance thermometry 
by Schultz [3]. It has recently received considerable 
attention as a device for establishing a platinum resist- 
ance scale of low temperatures [4]. A very extensive 
analysis of its application to this field has been made by 
Berry [5] who showed that it did not apply with suf- 
ficient accuracy to most of the thermometers con- 
sidered or even to most of those of high purity; i.e., 
having high values of the coefficient, a =(R373.15 
— Roz3,15)/100Re273.15, or low values of the reduced 
residual resistance, Wo = Ro/Ro7z3.15. 


' Figures in brackets indicate the literature references at the end of this paper 


In this paper we attempt to determine the ideal 
resistivity function by a new method; also we test 
the applicability of the SWK formula to a rather homo- 
geneous group of high-purity thermometers that have 
been previously correlated by other means [6]. Some 
preliminary work along the lines of the present paper 
was reported earlier [6]. 


2. Procedures 
2.1. Reduced Form of the SWK Equation 
Equations (1) and (2) with a and 6 taken to be con- 
stants are first transformed into relations involving 
resistance ratios. We follow Schultz [3] but use the 
notation of Berry [5] where applicable. Equation (1) 
becomes 


Wr=Wirt+ Wo — Wirhi+Tr) (3) 


where 


Wr = Rr/Roz3.15 
W ir = Rir/Ri 273.15 
Wo = Ro/Roz3.15 


and 
Py = (Ar — Aoz3.15W ir)/ p273.15V (1 — Wir). 


Introduction of eq (2) leads to 
(Ur)swx = AW ir|(1 + CWir) (4) 


in which A and C are constants which can be evaluated 
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if the resistance is known at any two temperatures, 
T, and 7>, other than 0 and 273.15 °K. The following 
expressions are found for A and C: 
Pr, —T2Wii/ Wir 

(T2—T)Wis 


A=T2(1+ CW i2)/ Wiz. 


(= 





The subscripts 1 and 2 refer, respectively, to 7; and 
T:. The constants A and C are related to the SWK 
parameters by the following close approximations: 


a~Cl/A 
b ] /A Wo. 


In the calculations to be described, 7; and T2 were 
always taken to be 20 and 90 °K; i.e., approximately 
the boiling points of hydrogen and oxygen, respectively. 


2.2. The Thermometers 


These consisted of 34 commercially made capsule- 
type thermometers plus thermometer T4 of the 
National Research Council, Canada. They are more 
fully described in the earlier paper [6]. All of the 
calibration data were expressed on the NBS 1955 
temperature scale. The group is of high purity as 


indicated by a> 0.003924 and Wolestd.) < 7X 10-4. 
The success with which they were correlated by the 


simpler empirical method of reference 6 suggested 
that they were rather uniform with regard to the forms 
of their deviations from Matthiessen’s rule and thus 
would provide as favorable a practical test of the SWK 
formula as one could reasonably hope to obtain. The 
thermometers are referred to by identification numbers 
from | to 35 assigned in order of increasing W 29.263. 


2.3. The Residual Resistance Ratio W, 


Experimental values of Ws. were available for five 
thermometers. These were adjusted to Wo using 
the data in Berry’s figure 1. For the remaining ones, 
W. was estimated from data in the region 10 to 20 °K 
using Method 1 in Berry’s appendix III. 


2.4. The Ideal Resistance Ratio Function W;7 


The function Wir for platinum has been estimated 
by various methods. A summary and discussion was 
given by Berry who indicated objections to the various 
earlier methods. He derived a new function, slightly 
different from the others. In the region 10 to 90 °K, 
this was found by extrapolating plots of W;— Wo versus 
W,. for a large number of thermometers to Wo=0. 
Since the values of W;—W, unfortunately lay in a 
band rather than on a line, only a few selected ther- 
mometers having low values of W;— Wo were used in 
making the extrapolation. Correspondingly, when 


Berry fitted the SWK formula using this Wir function 
to the various thermometers, the fit was satisfactory 
only for those few that had been used to derive the Wir 
function. 

The validity of the SWK model for the present group 
of thermometers is here tested by a different method, 
as follows. One may derive various arbitrary Wir 
functions from the Wr function of any one thermom- 
eter via eqs (3) and (4) by assuming various values of 
its parameters A and C. Equations (3) and (4) can 
then be fitted to all of the remaining thermometers at 
20 and 90 °K using these various Wir functions in turn, 
and the errors in the representations of Wr at other 
temperatures can be examined. If it is supposed that 
the SWK formula holds for the group, then it should 
be possible to find in this way a particular Wir function 
such that all of the remaining thermometers except 
the one used to generate the Wir functions are accu- 
rately fitted. Furthermore if various thermometers 
are used in turn to generate such an optimum Vir 
function, then the various optimum Wir functions 
should be found to be all identical and therefore equal 
to the “‘true’” Wir function. 

To carry out this test, nine thermometers were 
selected which, on the basis of earlier work [6], were 
judged to be representative of the set of thirty-five. 
The five of these for which experimental values of W4.2 
were available were used to generate the trial Wj7 
functions. 


3. Results 
3.1. The Function Wir 


When the procedure just described was carried out 
it was found that there existed a spectrum of suitable 
Wir functions. In going from one to another of these 
functions, the identities of the thermometers for 
which the fit was least satisfactory changed, as did the 
temperatures at which the errors were largest, but no 
uniquely superior fit of all the thermometers could be 
found. In figure 1 are shown some of the better Wir 
functions referenced to Berry’s Wir function as a base 
line. 

The fitting errors (see the next section) were roughly 
an order of magnitude larger than the precision with 
which the resistance thermometers are believed to 
have been calibrated. Consequently, the variation in 
Wir found by the present method is an artifact of the 
SWK model and is not chargeable to experimental 
error. 


3.2. Accuracy of the SWK Equation With Constant 
Coefficients 


The quality of the better fits obtained using Wir 
functions derived by the method of section 2.4 may be 
summarized by stating that the maximum temperature 
error of the nine representative thermometers in the 
interval, 20 to 273.15 °K, was never much less than 20 
mdeg, and the average of the maximum errors of the 
thermometers was never much less than 10 mdeg. In 
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FIGURE 1. Selected useable Wi functions plotted relative to that 


of Berry {5}. 


(The letters, a-b-c-d-e, identify functions referred to in table 1.) 


table 1 are summarized the errors of the SWK formula 
resulting from use of five of the better fitting Wir 
functions of figure 1. For comparison, Berry’s Wir 
function was also used. In the last two rows are given 
the corresponding errors of the simple empirical 
method evaluated earlier [6]. 


TABLE 1. Accuracy of the SWK formula compared to an empirical 


method of interpolation 


| 


Method of interpolation | Max. temp. error | 


} 


Average |ATna,\* 


Range examined 


SWK eqs (3-4) using 
the following # 
func¢tions 

u. Figure | 

b. Figure | + 20.3 12.4 
c, Figure | +21.3 10.6 
d. Figure | 8.6 
e. Figure | 2 10.6 
Berry [5] + 27.1 12.8 
Berry [5] | +27 12.7 


mdeg 


as 


oe 
wwwunww 


NNN NN 


= 


Empirical method ' 


Do 


* Average of the maximum errors of the nine representative thermometers 
© Reference 6, eq (10) using 20. 90. and 273.15 °K as 


calibration points 

All earlier studies of the resistance-temperature 
relation of platinum have encountered exceptional 
difficulties in representing it accurately in terms of 
temperature errors below 20 °K. It might reasonably 
be supposed that the SWK formula would be advan- 
tageous here because of its theoretical foundation 
and the constraint imposed on its behavior at the lowest 
temperatures by use of the residual resistance. How- 
ever, its behavior below 20 also is unsatisfactory 
inasmuch as the error at 10 °K was usually of the order 
of 0.1 deg. 


3.3. Further Characterization of the Thermometers 


Following Berry, we may characterize the thermom- 
eters via the function (Frexpu. calculated from eq (3) 
using the experimental values of W7, i.e., 


Wr-Wir _ 
(1—Wir)Wo 


(I Texptl sa 


This function, if it is to be fitted by eq (4), must rise 
monotonically—at first in proportion to Wir, i.e. 
roughly as 7°—and approach an asymptotic limit. 
Berry found that the thermometers which would be 
well represented could be identified as those with (a) 
Poo < 0.35 (C99 approximates the high-temperature 
limit of T'7), (b) Wo < 1 X 10>, (c) monotonically rising 
curves of 7. Only about 20 percent of the thermom- 
eters examined by him met all three requirements. 

It will be recognized that conditions (a) and (c) test 
only the uniformity of the group. They are not 
meaningful criteria of “normal” versus “anomalous” 
behavior inasmuch as the behavior displayed depends 
entirely on the form of the adopted Wir function. 

The present group of thermometers was examined 
in the above way using Berry’s Wir function. All 
of the thermometers met criterion (b). Eighteen out 
of the 35 met criterion (a). Fifteen of these 18, as 
well as two of the others, failed to strictly meet cri- 
terion (c). However only five of these (consisting of 
the thermometers with the lowest values of I'99) showed 
marked maxima in (I7)expu. ‘If a requirement based 
on calibration at common fixed points were to be 
adopted for the purpose of excluding thermometers 
having maxima in Iz (e.g., P90 > P54 > P27 > Peo), only 
these five would be excluded. Had members of this 
set of five been excluded from the nine used for most 
of the calculations, the SWK fits would have been 
improved for certain of the Wir functions derived and 
selected by the parameter-variation method, but would 
still have been inferior to that given by the empirical 
method of reference 6. The average fit of the SWK 
formula using Berry’s Wir function would actually 
have worsened. ( 

Figure 2 shows a sampling of the (Ir)expy curves. 
Numbers 17, 31, 5, 34, 23, and 11 are from the se- 
lected group of nine used for most of the calculations. 
Number 29 also was a typical curve: however, number 
12 was unique. Numbers 11 and 29 provided refer- 
ence functions used in the empirical method of ref- 
erence 6. Numbers 23, 31, and 34 were used in a 
preliminary study [7] in which it was found that they 
could be represented within about 1 mdeg by the 
method of reference 6. 


Table 2 lists the SWK parameters of the selected 
group of nine representative thermometers used in 
the derivation of Wir and in examining the fit of the 
SWK equations. It may be compared with Berry’s 
table IV. The range of the parameters a and b is 
very nearly the same as that covered by Berry although 
the range of Wo is much narrower. 
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FIGURE 2. Curves of (V'r)ex 9, for selected thermometers calculated 


using Berry’s Wy function. 


(The curves are labeled with the identification numbers of the thermometers.) 
TABLE 2. SWK parameters of the selected thermometers 


T T 


1 
Thermometer Wo X 10 A | , -|(- 


) | °( 
number | a a 
E 
| 
| 


4 
z 


] 
cw.) 


427 | 146 | 0 
512 |} 83 113 
666 | 86 166 
468 | 251 
681 135 430 
655 134 529 
485 231 926 
| 398 146 662 
| 498 | 386 1980 


1 


*Thermometer T4 of National Research Council, Canada 


4. Discussion and Conclusions 


The application of the SWK formula to platinum 
resistance thermometry has been based on the rea- 
sonable expectations that, because of its theoretical 
foundations, it should accurately represent a wider 
range of platinum purities and a wider range of tem- 
peratures than could purely empirical correlations. 
The work of Berry and others has already dashed 
these expectations. The failure of the present study 


to obtain a unique ideal resistivity function from the 
SWK model with constant coefficients as well as the 
failure to fit the model to any major part of the data 
used within the precision of those data indicate that 
the model has only qualitative significance. This 
supports the conclusions already reached. 

Thus the present study is of use mainly in examining 
the question: Would the SWK formula have some 
special usefulness if the characteristics of the ther- 
mometers were narrowly restricted? The present 
thermometers have narrowly restricted a and Wo 
values and are of uniform origin and style of con- 
struction. It has been shown that they can be ac- 
curately represented from 20 to 90 °K by a simple 
empirical method [6] using the same calibration tem- 
peratures as the present study except that the residual 
resistance is not required. Thus by these various 
criteria the present thermometers appear to be a quite 
homogeneous group. A_ correlation scheme that 
required any further restrictions on acceptable ther- 
mometer characteristics would not be justifiable on 
practical grounds. Yet the present study shows that 
the thermometers considered are a diverse group when 
categorized according to any of the following: (1) shape 
of the Ty function obtained using Berry’s Wir function, 
(2) magnitude of the SWK parameters, such as a, and 
(3) accuracy with which the SWK formula fits their 
resistance-temperature functions. Consequently, it 
seems unlikely that the SWK formula has any prac- 
tical utility for accurate thermometry or for classifica- 
tion of thermometers into acceptable and unacceptable 
types. 


The assistance of William J. Hall and Hans M. Roder 
in programming the calculations for automatic com- 
putation is gratefully acknowledged. 
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A liquid hydrogen cooled, tape-wound, axially segmented, high-purity aluminum electromagnet 
has been built and tested to 95 kilogauss with a magnet power consumption of 22 kilowatts. Formulae 
for the axial. radial. and circumferential stress components in this type coil are derived under the 
assumptions of isotropy and homogeneity. Both plane strain and plane stress analyses are given. 
The hoop stress is also derived from the “floating shell” concept using thin shell theory. The for- 
mulae from these analyses and thick-wall cylinder theory are used to calculate the hoop stresses in 
the experimental coil, and the results are presented for comparison. The theoretical burst field of a 
monolithic cylindrical coil disk is derived. 


Key Words: Axial stress, burst field, circumferential stress, cryogenic electromagnets, electro- 
magnets, magnetic body forces, magnetic field, plane strain, plane stress. 


1. Introduction 


Many current research programs require the application of high magnetic fields. Plasma 
containment, heating, and acceleration are examples. Cryogenic electromagnets are rapidly 
becoming the most convenient, economical means of producing these high fields. Purcell and 
Payne [1|' have discussed the advantages of cryogenic electromagnets and described in detail 
the coil construction treated herein. 

The maximum field attainable with an electromagnet is limited by the mechanical strength 
of the coil. The field strength obtained with a given coolant system and magnet power supply 
may also be limited by the heat transfer efficiency and heat capacity of the coolant and the con- 
ductor resistivity. Therefore, optimum coil design requires simultaneous consideration of stresses, 
heat transfer, and conductor resistivity. The first of these is of sufficient importance to be treated 
separately in this paper. 

In high-field magnets, the mechanical stresses may burst the coil. The coil may also fail by 
localized yielding which can result in increased coil resistance and electrical shorting. There- 
fore, it is necessary to determine the stress field for a specific magnet construction. The coil 
construction of particular interest here is the stack of tape-wound epoxy-bound cylindrical disks. 
The coil disks are alternately stacked with fiber disks which provide electrical insulation and sector- 
like openings for radial coolant passages. The stack of coil disks and separators is shown on figure 
1. This photograph was taken following electrical and intermediate structural failure of the magnet 


*This study was conducted on Contract AT-(49-2)-1165 for the U.S. Atomic Energy Commission, PROJECT SHERWOOD 
' Figures in brackets indicate the literature references at the end of this paper 


FIGURE 1. Stack of magnet coil disks and disk 
separators with top four coil disks removed (after 
electrical and intermediate structural failure of 
the magnet at 95 kG). 





at 95 kilogauss (kG). The foil conductor is assumed to be tightly wound with alternate layers of 
paper insulation and then vacuum impregnated with epoxy resin to form a monolithic cylindrical 
disk. The paper and epoxy combined are not permitted to exceed 5 percent of the disk volume. 
The epoxy-paper tensile and compressive yield strengths exceed those of the conductor. This 
construction technique yields to the simplified analysis of an isotropic, homogeneous cylindrical 
disk. Although the epoxy-bound laminated construction cannot be perfectly isotropic, the assump- 
tion of isotropy is considered representative of the coil construction. These cylindrical coil disks 
are analyzed herein and formulae are developed for the axial, radial, and circumferential stress 
components. The radial and circumferential stress components are derived under the assumptions 
of both plane stress and plane strain. The hoop stress is also developed from the “floating shell” 
concept which assumes that each turn of the coil disk is an individual thin wall cylinder. This 
type analysis should be appropriate for loosely wound, poorly bonded coil disks. The formulation 
for predicting the burst field of a monolithic cylindrical coil disk is also developed. Stress and 
burst field design curves based on the analyses are presented. 


Magnet stresses have been given considerable attention in the literature. For the most part, 
only the circumferential (hoop) stress component has been examined. This is reasonably justi- 
fiable since the circumferential stress is the largest component. Post and Taylor [2] point out that 
the magnetic pressure is proportional to B?/87, where B; is the axial field component. They show 
that the main incremental body forces are proportional to the radial derivative of the magnetic 
pressure. Bitter [3] extends the magnetic pressure concept to a very useful form by assuming the 
coil to be a thick wall cylinder with internal pressure B?/87 (Bo is the field at the geometric center of 
the coil). This calculation scheme is quick and easy to use and, as will be shown later, provides 
a rough estimate of the coil hoop stresses. However, it should be noted that this calculation 
scheme is valid only for the hoop stress (a) component. If a shear failure theory is to be assumed 
for the coil, it is necessary to also know the axial (o;) and radial (o@,) stress components. — If the 
maximum normal stress failure theory is employed, then a; will suffice. It will be shown herein 
that o, is not negligible and can cause coil failure or excessive deformation under the disk 
separators. 


Detailed stress analyses for various type magnets have been given by Cockcroft [4], Daniels 
[5], Giauque and Lyon [6], and Wells et al.. [7]. None of these are deemed applicable to the 
magnet construction of interest for the following reasons: Cockcroft [4] assumed that the mag- 
netic body forces were balanced at the outer periphery by an external pressure. This may be 
true for coils rigidly reinforced at the outer edge but no reinforcement was used for the coils 
described here. Cockcroft also made the simplifying assumption that the solenoid was a ho- 
mogeneous, isotropic body, but neglected the largest stress component (¢-) in his analysis. In 
addition, the formulae developed are dependent upon the calculation of the mutual inductance of 
a solid coil and a circle in the end plane. Although tables are given for specified variations in 
coil geometrical parameters, the calculation is laborious. Daniels [5] also utilizes the mutual 
inductance of one turn with the entire coil to derive stress formulae for coils fabricated by stacking 
tape-wound cylindrical disks. Since the disks were not considered to be structurally reinforced, 
the problem is similar to the one under consideration in this paper. Daniels [5] considered the 
problem as two dimensional and ignored stresses in the axial (Z) direction. The windings are 
assumed to be homogeneous and anisotropic because the coil was fabricated with nylon insulation 
between adjacent turns, and the composite was not integrally bonded in any way. It is assumed 
that the coil cannot unwind due to friction between adjacent turns. Daniels [5] includes the all- 
important hoop stress and the conditions of stress compatibility in the analysis. It is concluded 
that each turn of the coil is in equilibrium under its own body forces and that axial compression 
does not affect the magnitude of the stresses in the radial and circumferential directions. Daniels 


[5] also concludes that evaluation of the formulae developed is very laborious and suggests experi- 
mentation with models as an alternative to calculating the stresses. Giauque and Lyon [6] treat 
the problem of helical-wound cylindrical coils. Each helical layer of the solenoid was separated 
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by an electrical insulating material which also defined the cooling fluid annuli. Thus, the coil 
could not provide uniformly distributed resistance to radial forces and was prevented from un- 
winding only by frictional drag or deformation of the conductor. Each helical layer of the coil 
is assumed to be self-supporting until the stresses in the conductor exceed the conductor yield 
strength. Giauque and Lyon [6] also use the mutual inductance method in the development of 
formulae for the axial and circumferential stress components. A simplified expression for the 
calculation of the mutual inductance between current element and coil is given for nonhomoge- 
neous coils. This simplified expression is reported to be useful for the quick evaluation of maxi- 
mum stresses in a design feasibility study. Under the foregoing assumptions Giauque and Lyon 
[6] neglect the radial stress component, assume a maximum normal strain failure theory, utilize 
the mathematical expression for calculating mutual inductance, and derive very useful equations 
for the hoop, axial, and failure stresses in unbonded helical-layer cylindrical coils. Wells et al., 
[7] have extensively examined the distribution of magnetic forces in toroidal solenoids. Wells [8] 
has also developed very useful formulae for the computation of hoop and radial stresses in Bitter 
magnets (current density is inversely proportional to the coil radius). The stress formulae devel- 
oped herein are dependent upon the magnetic flux density (B). The recent data of Brown et al., 
[9, 10] provide numerical values of B for any solenoid. Therefore, stress calculations using the 
derived expressions are simplified. 


2. Analysis 
2.1 Nomenclature 


B= vector magnetic flux density 
Bo = magnetic flux density at the geometric center of the coil 
C=constant of integration 

E = Young’s modulus of elasticity 

F =force generated by magnetic field 

]=current per turn of conductor 

J =current density in conductor 

L=half length of cylindrical coil 


m= absolute value of the slope of B,/Bo versus y plot 
V = number of turns in coil 


Pm= Magnetic pressure 
r=radius at some turn in the coil (radial coordinate) 
ri= inside radius of coil 
ro= outside radius of coil 
t= wall thickness of a thin cylinder 
u= measure of radial displacement 
V = conductor volume 
y= absolute value of B,/Bo at r=ri{y=1) 
Z = denotes distance along axis from center of coil (axial coordinate). 


a. Greek 


a= ratio of outer to inner coil radius (= ro/ri) 

B=ratio of coil half length to inner coil radius (= L/r;) 
y=ratio of turn radius to coil inner radius (= r/r;) 
€=unit strain in conductor 

0= circumferential coordinate 

= Poisson’s ratio 

o =unit stress in conductor 





oa, = yield stress of conductor 


= magnetic flux geometry parameter 


{=£ 1 ee } 


b. Subscripts 


c= denotes circumferential direction 
r=denotes radial direction 
z=denotes axial direction 


c. Superscripts 

= denotes plane strain condition 

= denotes maximum (peak) value 

= denotes average value. 
The analysis is performed under the following general assumptions: (1) The tape-wound, 
epoxy-bound axially segmented cylindrical coil is considered to be homogeneous and isotropic 
in each disk. (2) The radial and circumferential stresses are derived under the assumptions of 
plane stress and plane strain. The coil is comprised of a stack of relatively thin disks which are 
separated by a disk-shaped insulating material with sector-like openings for coolant flow. There- 
fore, it would be extremely difficult to exactly determine the effect of axial compression on the 
radial and hoop stresses. The assumption of plane stress is judged most applicable since the 
major portion of each disk face is exposed to coolant. Axial stresses due to axial magnetic force 
are transmitted from disk-to-disk via the spokes in the disk separators. The axial compression 
is magnified by the reduction in compression area under the disk separators. 


(3) The coil is sym- 
metrical with symmetrical flux distribution and uniform current density. 


(4) Frictional drag 
under the spokes of the disk separators does not resist radial deformation of the conductor. 
(5) Buckling or warping of the disk does not occur. This also means that no bending takes place 
between spokes of the disk separators (due to axial compression). (6) The stress formulae devel- 
oped are applicable only so long as the conductor stress remains below the proportional limit. 
(7) The maximum normal stress failure theory is chosen because the hoop stress component dom- 
inates. 


The maximum normal strain or maximum shearing stress failure theory is more appli- 
cable to the conductor material immediately below the separator spokes; however, the material 


under the spokes comprises a small portion of the disk volume, and the high compressive stresses 
can be redistributed to the unstressed conductor. 


Consider the cross section of a coil with a magnetic flux linkage as shown in figure 2a. The 
horizontal (radial) component (1) of the vector field induces vertical (axial) forces on the conductor. 
The axial component (2) of the magnetic field induces radial forces on the conductor. The dis- 
tribution of the axial field component across the solenoid radius is shown on figure 2b. The field 
variation along the solenoid axis is shown in figure 2c and the radial field distribution on figure 2d. 
The axial field creates the bursting effect and varies almost linearly across the coil radius. The 
radial field is nearly sinusoidal over the coil radius and increases with distance along the axis (Z) 


for Z<L. 


Now the incremental force on a conductor segment in a magnetic field is given as 


WV 10 (1a) 
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dF _ BJ. 








FIGURE 2. (a) Cross section of solenoid showing 
magnetic flux components and associated body 
forces. (b) Distribution of the axial field com- 
ponent over the coil radius. (c) Distribution of 
the axial field component over the coil axis. 
(d) Variation of the radial field component with 
radial and axial position. 





The current density for the solenoid in figure 2a is 


VI 


~ 21 (ro — rj) 


J 


Also, the field at the geometric center of the coil may be written, 
_ 4nNI 


Bo= 10(2L) 


®. (le) 


Combining eqs (1a), (1b), and (1c) we obtain the magnetic body force per unit volume of conductor; 
FB es 
dV Att(ro— ri)®P 

The radial force component is 
dF, = B.By 
dV 4a(ro — r)®’ 


and the axial force component is 


dF; BrBo__ 
dV Amro — ri) 
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FicurE 3. Distribution of magnetic body forces on a coil 
disk at midplane (Z/L=0). 























Equations (2b) and (2c) are valid for any point in the coil cross section. Equations (2b) and (2c) 
indicate that the radial and axial forces are linearly proportional to the field components. Refer- 
ring to the field distributions in figure 2, we can sketch the loading imposed on a typical disk in 
the coil. The distribution of the magnetic body forces is illustrated on figure 3. 


: Bio, : 3 : . : : dF, 
Radial equilibrium of a differential element of conductor with magnetic body force Ww 
is satisfied by 

do, __ dF, 


Oc—Or—Tr i, a (3) 


See Timoshenko and Goodier [11]. 


2.2. Plane Stress 


For the condition of plane stress o-=0. Hooke’s law for the condition of plane stress may 
be written [12] in terms of stress and displacement as follows: 


- (l—p*) (r Vedr 


— | {t+ “I. 
= (1 — pw?) (dr a 


Combining eqs (3) and (4) we obtain the equilibrium equation in terms of displacement 


age dr re 


du , ldu 5--()5* 


Equation (2b) must now be evaluated in order to solve this differential equation. It will be shown 
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fa Ff 


later that B-= Bo ( 


) provides fair results in computing the stresses at the midplane; how- 
ro Ti 


ever, to maintain generality we shall let 
B:=Bo[y+m( —y)]. (6a) 


where y=the value of B,/Bo at r=rj (fig. 2b), and m= absolute value of the slope of B./Bo (fig. 2b). 
Combining eqs (2b) and (6a) we obtain 


2 P — 
dF, _ Boly+ mil y)). (6b) 
dV 4a(ro— ri)P ' 





Substitution of eq (6b) into eq (5) yields 


fe] .. <=K ly+m(1—y)]. 
a 0 O6rdadr OTF 


where 


47 Er(a— 1) 7b 


The left side of eq (7a) is recognized as Euler’s equation and the complementary solution is written 
. C2 

u, =Cywt+— (8a) 
" 


The particular solution to eq (7a) is 


(8b) 


ay) mr 
Me =H 5 — —}, 
- 3 8ri 


and u=u;+us. The constants of integration C; and C. are evaluated trom the boundary conditions 
o,=0 at the inner and outer radii of the disk. 


Invoking these boundary conditions and combining 
eqs (4b) and (8) we obtain 


- KQir; 4 KQ2r3 
O=-" =- 


where Q; = (-- ———: {Qo(1 — w) + BV + 2)+ m(7 + Sy)}, and 
+ bh 
02 = = - (—) {8y{ w+ 2) (1 — a) +m[(7 + 5p)— Bal uw + 2) + 0°(9 + 3y)]}. 
‘ et Na 


Substituting into eqs (4) we obtain the following formulae: 
Bi | (1 —p) 
——m  1()(1 +) + Oo —3— + Bmy(1+3p)— ; 2 O. 
V6 7a — lb |<! ae Q» y my*(1 bh) 8y(y + m) (1+ b)}. (9a) 


Be | , : 
ie O11 + p)— Qe “— or + 3)— 8yy + m) (e+ 2)}. (9b) 


The stresses are in dyn: 


s/em? when Bo is expressed in gauss. To convert to psi multiply B? 
x 14.5 x 10-8, 
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2.3. Plane Strain 


For the case of plane strain €,=0. Pursuing the method used in deriving the plane stress 
formulae we obtain for plane strain, 
B (1 — 2p) | 


Te = yf =~ Ge ———— + Oyl y+ a) (1 + a) — Savy] + 2m) 7, (10: 
967r(a — 1) (u— 1)® 1Qr—Q: y Ww lad f BS) a) 





B {ox— oe CA) ay-+m) (2—p)—3my°3—2p)| 
967(a — 1) (u—1)® (~' . yy 2— ph): 2m). 





(10b) 


where 


Ot =—8y(2—) (4+ 1) + m(7— 2p) ( = + me {8a(2 — ) 


at] a’ -- : 





a“ —]| 


2 
ar 





+ 3a°(2u —3)}, and QF = {8y(2 — pw) (1 — a) + m[(7 — 2p) — 8al(2 — p) 


(a? — Qa —1) 
— 3a°?(2y — 3)]}. 
2.4. Simplified Plane Stress 


If the hoop stress at the midplane only is desired for design calculations some simplifying 
assumptions can be made. 


Set B.= By (2—), 


To— Ti 


which requires that y= 1 and m= in eqs (9). The resultant equations are 


] 
a—] 
B2 


l— 3y?(1+ 3p) a 
= et ap {Hl tm) +e (SE) + aya )at+20}> ae 





> 


1— 3y2(u + 3) a 
= Sema yp {Ml +H) ve ( a+ i — By ( i) m+ 2)} ina 





where Wi; = Cres |e 1—p)+8(u+ 2)+ || 7+ 5u)|, 





and ws = (— (— -) || {8a( + 2)— a%(5u+7)—3u—9}. 


2.5. Bitter Method 


Bitter [3] has suggested the use of formulae for the thick-wall cylinder subjected to internal 
" Ss . 2 eRe . . . u . . : 
pressure 83/87. Utilizing this method we obtain for stresses at the coil midplane, 


_ Bil t+(aly? 


Oc ; 
87 a?—] 


_ Bl- (aly? 
8r a] (12b) 
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As mentioned earlier, eq (12b) is not applicable because the boundary condition at the inner radius 
is not satisfied. 


2.6. Floating Shell Method 


The magnet fabricated at NBS admitted annular cracks of light upon completion, indicating 
imperfect epoxy impregnation. Thus, it seems worthwhile to calculate the hoop stresses as if 
each conductor turn is self-supporting and unrestrained, but closed as in a thin wall cylinder. The 
magnetic pressure may be written 


_ dF, _ (dF, 
dpm= "a4 = (Gp) a 


Restricting this calculation to the midplane we set 


B:=Bo(2—") 


"Fe 
dr=t, dpm=Apm, and combine eqs (2b) and (13a) with thin shell theory to obtain 


B2y(a— y) 
a = (13b) 


This hoop stress reaches a maximum at the radius r=ro/2._ Then by the floating shell method 
the maximum hoop stress 


B2 


‘ . > 
(Ochy=2 = ( ): 


al 


-  167r@ 


FicurRE 4. Equilibrium of a coil disk when the entire ; 
cross section of the disk is stressed to the yield point. Oy (ro-ri)dz Oylro-rj)dz 
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2.7. Burst Strength of Monolithic Disks 


It is assumed that failure occurs when the coil material becomes plastic throughout. Although 
it has not yet been demonstrated, a, and o; are relatively low compared to o-, and therefore the 
maximum normal stress failure theory is selected. Referring to figure 4, equilibrium is satisfied 
when 


20,(\To— Fi de = | dF, cos @. 
. (15a) 


Choosing unity thickness (dz= 1) and combining eq (6b) and (15a) and performing the integration 
we obtain the theoretical bursting field 


73(a — 1) 
B2= 1.73(a— 1Poy® 





om ’ (15b) 
(y+ m) (a?—1)—=> (a? — 1) 


where By is in kG and a, in pounds per square inch. 


If for the sake of simplicity we substitute y= 1: m=1/(a— 1) into eq (15b), the bursting field 
at the midplane is 
5.2(a— lhoyP 


(a+ 2) 





By= 


2.8. Axial Compression 


The axial magnetic body forces are transmitted from disk to disk by the spokes of the disk 
separators. Thus, the forces are additive from the ends to the center so that the disk in the mid- 


plane must support the axial body forces of the entire coil. Axial equilibrium of a differential 
element require® that 


do. dF; 


dz dV (16) 


In order to obtain a we must evaluate B,(r, z) in eq (2c). It will be shown in the sample calcula- 
tions which follow that the radial field component is approximated by an expression of the form 


oe a a . —] 
Be, Kot (K— Ko) {sin 7 (~—) , 


(17a) 


where Ky = (Br/Bo)=r “ Cites for each Z/L of a fixed coil geometry, and K= (7) (peak 
0 = 
Y 





= —ar+l 


9 


value of B,/Bo for each Z/L of a fixed coil geometry). Ko and K — Ky can be linearly related to the 
axial dimension Z by: 
Ko=m2Z/L, K—Ko=m2Z/L, (17b) 


where m;, is the slope of the Ko versus Z/L plot, and mz is the slope of the (K — Ko) versus Z/L plot 
for a fixed coil geometry. Combining eqs (17a) and (17b) we obtain for the radial field, 


By _Z F ( 
By L m,; +> Mz Sin 7 - 
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Combining eqs (2c), (16), (17c) and performing the integration the axial stress is derived, 


_ Be 
8rLr{a— 1) 


i 


: ak. , 
{m+ ms sin (2) +fin+C. (18a) 
nea 
Since external forces are not being considered C=0 and f(r) is obtained from the boundary con- 
dition (o,)z-,=0. Performing the indicated algebra we arrive at a solution for the axial com- 
pressive stress 


BBL? — Z*) =) 
= — 1m, +m 8 fj}. 
O: Bulan le | tm sin wl (18b) 


This expression is applicable throughout the axial dimension of the coil. The axial stresses are 
largest at Z=0, and 
= BBS 
(Oz)z=-0 = 


“i m,+ my sin 7(%=) (18e) 


\a— ] 
a+] 
2 


rigid plate so that the axial stress is equally distributed we find the average axial stress (@-) from 


The radius at which (o,)z-9 becomes maximum is y= Assuming that each disk acts as a 


—_ ] [° 
o:= ody. (19a) 
a-—l1j, 


Performing the indicated integration on eqs (18b) and (18c) we obtain 


ote =e" 12m 
8rLa— ly® |" ’ 


2 


a = BB? 2m2 
(@ne-0=— see pm tt (19¢) 


In order to obtain the maximum average compressive stress under the spokes of the disk separa- 
tors at the midplane we need only to multiply the results of eq (19c) by the area magnification ratio 
(A aisk/A spokes): 


2.9. Sample Calculations 


In order to compare the various formulae and illustrate their use, sample calculations are 
performed for the coil geometry of the NBS magnet. 

a=" = 3.67; p=== 1.87; re=1.5 in.: ro=5.5 in.; L=2.8 in.; w= 0.33 (aluminum); o, = 4500 

i i 

psi; ®=0.64. To evaluate eqs (9), (10), and (18c) we must determine the values of y, m, m;, and 
m». All calculations will be performed for the midplane (Z=0) where the stresses are largest. 
All data for the radial and axial field components are taken from the tabulations of Brown et al., 
[9]. Brown and Flax [10] recently presented a graphical method for determining the magnetic 
field components of any finite coil by superposition of four semi-infinite solenoids. The axial field 
component (B,/Bo) is plotted for various coil parameters (a, 8) in figure 5. The dotted line indi- 
cates the extrapolated curve chosen to fit the coil parameters of the NBS magnet. Since B is 
shown to have little effect on the slope, the selection of the dotted line is based on linear extra- 
polation between the curves for a=3 and a=4. From figure 5 we obtain y= 1.04 and |m| =0.45. 
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FIGURE 6. Variation of end-plane radial field 
with coil geometry 








3 
yy er7c7ry 


FIGURE 5. Variation of midplane axial field with 
coil geometry 


On figure 6 is plotted the radial field component for various coil parameters. It is difficult to extend 
these curves to intermediate values of a and B. The effect of smaller B is to increase the magni- 
tude of B,/Bo with little phase shift. Thus we plot on figure 7 the B,/Bo data for a= 4,8 =2, and 


various Z/L. We note that the curves are essentially sinusoidal from y=1lto y=4=a. Recalling 
that 


B, 
+ 
Jaa (ay 


¢ = (2: 
oie (Be) aac’ 


we plot Ko and (K- Ko) against Z/L on figure 8. From the Ko plot we obtain the slope m; = 0.213, 
and from the (K— Ko) curve we obtain the slope m2=0.15. Substituting these values of m; and 
mz into eq (17a) we check on the validity of the assumptions. The B,/Bo is calculated from eq (17a) 
and plotted on figure 7 for comparison. The calculated values plot high for lower values of Z/L: 
however, recall that the effect of a lower B (fig. 6) is to increase the B,/Bo. Also, the fit is good 
for higher Z/L where the largest body forces are generated. The various formulae are evaluated 
and plotted fof’@omparison on figure 9. From eqs (15b) and (15c) we calculate the bursting field 
as 93 kG and 84 kG, respectively. Using eq (19c) and an area magnification ratio of 4.4 we obtain 
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FIGURE 8. Ko and (K— Ko) as a function of axial 
y 2e/tj position for a fixed coil geometry (a= 4, B= 2). 


» 7. Radial field as a function of axial position 
for a fixed coil geometry (a= 4, B= 2). 


the average axial stress (psi) under the separator spokes as — 0.843 B2 (By in kG). The maximum 


hoop stress [(o-)y=1 from eq(11a)] is plotted for various coil parameters (a, 8) on figure 10. Figure 
1] shows the bursting field [Bo from eq (15c)] as a function of a and B. 


2.10. NBS Magnet Failure 


The NBS tape-wound aluminum magnet failed electrically with intermediate structural dam- 
age at 95 kG. The magnet failure is attributed to excessive deformation at the outer radii which 
normally precedes bursting. Therefore, the failure field should not be much different from the 
actual burst field. The failure field agrees well with the 93 kG burst field calculated from eq (15b). 
The actual burst field is expected to exceed the theoretical value due to the frictional restraint of 
the axial compression (which was neglected in the analysis). Eight of the fourteen disks suffered 
circumferential cracks at radii of about 3.75 in. indicating excessive deformation at the outer radii 
and structural inhomogenieties. Some of the cracks were 0.125 in. wide. There is evidence of 
some radial deformation over the entire coil radii. The deformation of individual coil disks varied 
with coil axial position. The disks nearest the coil midplane were exposed to the largest magnetic 
pressure and suffered greater damage. The pronounced failure at an intermediate radius has been 
attributed to an intrinsic weakness in the coil disk construction. Several of the disks were rather 
loosely wound, and consequently some sizeable voids in the outer turns of the disk were epoxy- 
filled. All of the disks that failed possessed these voids from which the circumferential cracks 
propagated. Other disks in the coil had been tightly wound and suffered little or uniform stress 
damage. It would be quite difficult to determine what type structure the loosely wound epoxy- 
bound disk construction represents; however, it is estimated that the stresses should be interme- 
diate to floating shell and plane stress analyses. The hoop stresses for the floating shell (eq (13b)) 
and simplified plane stress (eq (11a)) conditions are superimposed on the conductor yield strength 
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FIGURE 9. Circumferential, radial and axial stress distributions at the midplane 


for the NBS magnet (a= 3.67, B= 1.87) as computed from the various theories 
disc ussed 


on figure 12. The variation in radii at which failure occurred is also indicated on figure 12. The 
conductor yield strength (0.2% offset) was experimentally determined to be 3300 psi at room 
temperature. The yield strength was then inferred to be 4500 psi at 20 °K (the operating temper- 
ature of the magnet). The actual strength of the material throughout the coil was unknown because 
the magnet had been tested at fields up to 70 kG several times prior to the 95 kG test. Thus, the 
material in the inner turns had been strain-hardened. Increased magnet resistance substantiated 
strain-hardening. The cyclic effect of strain on resistance and strength is under study; however, 
it is well known that the yield strength of pure aluminums can be increased three-fold at 20 °K by 
strain-hardening. A three-fold increase in yield strength of the inner turns of the coil disks would 
explain the circumferential cracks at intermediate radii when the magnetic field was increased 
from 70 to 95 kG. The epoxy voids provided the inherent weakness to permit the outer turns to 
expand radially outward from the strain-hardened inner turns. From figure 12 one could only 
conclude that the loosely wound turns of the epoxy-bound disks conformed more closely to the 
floating shell analogy. The shaded area on figure 12 indicates the region where the actual stresses 
for the loosely wound coils may plot. Some of the disks near the midplane were also badly dam- 
aged by axial compression. Extrusion of the conductor into the flow passages varied from 0.001 
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FiGURE 10. Maximum hoop stress as a function of FIGURE 11. Burst field as a function of the coil 
the coil parameters, a, B. parameters, a, B. 


to 0.014 inches depending upon the disk axial position and the apparent aluminum-epoxy com- 
posite strength. Although the disks are separated by 0.025 inch fiber spacers a 0.004 in. feeler 
gage could not be inserted in some of the flow passages upon post-run inspection. 


3. Conclusions 


For quick design computations the method of Bitter, eq (12a), though nonconservative, is by 
far the most convenient. If a more precise analysis of both circumferential and axial stresses is 
desired eqs (9a) and (18c) are recommended. Equation (lla) is more convenient than eq (9a) and 
provides conservative results for hoop stress calculations. 

Equation (13b) is suggested for hoop stress calculations on coils which are loosely wound and 
not integrally bound. The radial stress component is so small that it can be neglected in the coil 
construction described herein. Equation (11b) is recommended for those coil constructions where 
the radial stress is pertinent. Equation (15b) predicts the burst field of the tape-wound coil disks 
with reasonable precision. Equation (15b) and the simplified eq (15c) predict burst fields lower 
than those to be expected in practice. The actual burst field of a coil may be increased by strain- 
hardening the conductor through cyclic testing to the requisite fields. 
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FIGURE 12. Hoop stresses (as computed from the 
most dissimilar theories) and yield strength are 
superimposed on the failure radii of the NBS 
magnet. 
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Centerline Correction for Precision Roughness Specimens 
J. L. Chamberlin 


(May 24, 1965) 


In graphic determinations of roughness height (a measure of surface texture) a specific centerline 


about which the total area is measured is required. 
centerline is subjective and a repositioning was generally required. 


The first approximation of the position of this 
Here presented is a method for 


determining the roughness height without resorting to an additional area measurement. 


Key Words: Centerline corrections, precision roughness specimens, roughness height. 


Roughness height, one measure of the texture of a 
surface, is the arithmetic average deviation of the 
undulations of the surface, and is expressed by 


1 L 
Ra=t [lvls (1) 


where 


R, is the roughness height, 
L is the length over which the roughness height is 
measured, 
y is the amplitude of the undulations measured nor- 
mal to their centerline, and 
dx is the infinitesimal length measured along the 
centerline. 


For uniform repetitive wave forms. like those com- 
prising a precision roughness specimen, the center- 
line is taken as a line parallel to the general direction 
of the profile section of the surface tested, and so 
positioned that the area contained between it and those 
portions of the surface’s profile which lie above it 
are equal to the area between it and the profile por- 
tions lying below it. 

The .length over which the roughness height is 
determined is generally taken as equal to a roughness- 
width cutoff or some whole multiple thereof. The 
roughness-width cutoff is the greatest spacing of repeti- 
tive surface irregularities which are to be included in 
the measurement for roughness height. 


In the calibration of precision roughness specimens 
a graphical interpretation of roughness height is often 
employed. This results from performing the indi- 
cated integration of the roughness height defining eq 
(1). This integrated form is 


] 
Ra=7 (Ar+ Ap) (2) 


where 


Ar is the area between the centerline and that portion 
of the profile lying above it; 

Ax is the area between the centerline and that portion 
of the profile lying below it. 


In practice the surface to be measured is recorded 
using a surface finish analyzer. Such analyzers slowly 
traverse a sharp stylus across the surface. The 
undulation of this stylus relative to a straight datum 
is amplified and recorded, resulting in a magnified 
profile section of the surface. A centerline is posi- 
tioned on this recording by visual estimation, and the 
length over which roughness height is to be determined 
is marked off to the nearest whole cycle of the undula- 
tions. The areas above and below the centerline are 
then planimetrically determined. If the measured 
areas are the same then application of eq (2), being 
mindful of the appropriate recorder magnifications, 
yields the roughness height. If the measured areas 
are not the same then it has been general practice to 
reposition the centerline and redetermine the areas. 

If three conditions are met. this remeasurement 
may be circumvented by making a correction for the 
malplacement of the centerline. These conditions 
are as follows: First, the profile traced is comprised 
of uniform, repetitive wave forms; second, a minimal 
amount of waviness is present; and third, the difference 
between measured areas is small compared to their 
sum. Waviness is the long wavelength undulations 
of the surface upon which the repetitive wave forms 
are superimposed. 

For consideration as a precision roughness specimen 
the first and second conditions must exist. The third 
condition follows as a natural consequence of the first 
two conditions since the more uniform the specimen, 
the easier it is to position the centerline. On typical 
specimens it is common to have the difference in areas 
less than 5 percent of the total area on the first visual 
estimation of the centerline position. 
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The required correction also is obtained as a con- 
sequence of the above conditions. That is, the rough- 
ness height of a specimen may be obtained from the 
measured areas and a knowledge of other surface 
parameters. In this development the measured values 
are primed and the desired values are not. Thus, the 
total measured area, A’, is 


A'=A;,+ Ap, (3) 
and the desired total area for the roughness height, A, is 
A=Art+ Ap. (4) 


The error, E, in determining the correct area for appli- 
cation of eq (2) is 


E=A—A'=Ar+Ap—(Ay7+Ap). (5) 


The total correction is equal to the sum of the indi- 
vidual corrections of each cycle making up the length 
tested, or n times the average correction where n is 
the number of cycles in the length investigated. In 
this development figure 1 is to represent the average 
cycle and the corrections are developed for this case. 

For the small correction to the total area involved 
it is sufficient to consider areas a’abb’ and b'bcc’ as 
trapezoids with altitude Ay. 

The area above the centerline is for a single cycle 


a : 
ar= ar —> [(b'—a')+(b—a)| (6) 
and below the centerline 
. Ay 
ag = apt > [(c’ —b')+(c—))| 


where Ay is positive when 
> 42. 


(The lower case represents a measure referred to the 
single average cycle.) The error is thus 


E=n\—4Y 


> (b'—a’)+(b—a)| 
Ay 


> | 


+ 


(c’—b')+(c—)]t- (8) 


FIGURE 1. Representation of the average cycle of a precision rough- 
ness specimen as measured about a malplaced centerline. 
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(b—a)=(c—a)—(c—)) 
(b' —a')=(c' —a’')—(c' —b') 
(c'—a')=(c—a)=l'=L'/n 


(c—b)=(c' —b')+ Ayvicot 0. — cot Oy). (9) 
Thus, the error becomes, upon substitution of eqs (9), 
E=n{—Ay[l' —2(c' — b') — Ayicot 0, — cot Oy)] }(10) 


By definition 


4; — Ay = AA’. 
Therefore the measured area difference is 


AA’ ="d) [(b' —a’)+(b—a)+(c' —b')+(c— b)] 


or upon reduction 


AA’=L' Ay. (12) 
The area required for the roughness determination 
is therefore 


te oth ~ 90? —§') 


L’ 


_ BA cot A — cot Oy)). (3) 


L’ 


The desired roughness height based on measure- 

ments made with a malplaced centerline is therefore 
AA’ 

Ra=Ra —F- 


L'/n ‘& 





[ 2(c' — b’) 


AA’ (cot 6. — cot I. (14) 


L'/n 


In practice the parameters of the average cycle 
used in the adjustment of the roughness height are 
the averages of individual values measured at points 
randomly chosen along the trace. Because of the 
uniformity of most precision roughness specimens it 
is satisfactory to use the cotangent of the average 
angle rather than the average cotangent of the angles. 
For example, the cotangent of the average angle and 
the average cotangent of the angles for four precision 
roughness specimens that were being tested during 





the writing of this paper are shown in table 1. The 
values given are for five individual measurements 
made along the recorded trace. The recorded traces 
contained eight cycles for specimens 1 and 2 and ten 
cycles for specimens 3 and 4. 


TABLE 1 


Cot Oy 








Specimen cot By Range of 6, cot B- Range of 6- 





— 0.896 
— .B88 
— .202 


Degree 
1.155 5 
1.108 
0.233 

















— 186 2.3 246 


To illustrate the use of the correcting equation, the 
roughness height of a trace similar to that shown in 
figure 2 was measured using a visually estimated cen- 
terline and four malplaced lines. The malplaced 
lines were approximately 0.05 and 0.10 in. above and 
below the centerline. The results are tabulated in 
table 2. 


TABLE 2 


Position 4, An (ec —b) cot By 


sq. inches| sq. inches| Inches Inch 
1.378 2.123 | 7.523 
1.612 1.872 | 7.521 48: AS 3211 
1.775 1.704 | 7.524 463 : 3213 
1.932 1.567 | 7.527 Ai A 3271 
2.131 1.383 | 7.526 y 1567 3325 


0.10 in. above 
0.05 in. above 
Centerline 
0.05 in. below 
0.10 in. below 


























Each of the average values is based on eight indi- 
vidual measurements. The trace contained eight 
cycles. 


FIGURE 2. Typical trace of the precision roughness specimen used 


in illustrating the correcting equation. 


To obtain the specimen’s roughness height, the 
roughness height obtained from measurement on the 
recorded trace needs to be divided by the vertical 
magnification of the recorder. In this example, the 
vertical magnification was 2112 x. 

In this example 24 min were required to position 
a line and measure the areas about it. It required 17 
min to measure the various parameters and to calcu- 
late the correction. A time savings of 29 percent could 
have been achieved in this case over repositioning 
for a malplaced centerline. 


The use of this correcting equation has eliminated 
the centerline repositioning for graphic roughness 
height determinations. The time required to measure 
the necessary parameters is less than that required 
to reposition the centerline and make a new plani- 
metric determination. Depending upon the complex- 
ity of the trace being measured, time savings up to 
75 percent are possible. 
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Electric Currents and Potentials Resulting From the Flow of 
Charged Liquid Hydrocarbons Through Short Pipes 


M. R. Shafer, D. W. Baker, and K. R. Benson 


(June 21, 1965) 


The electrical currents and potentials produced in pipes of intermediate and very high resistivities, 


by the flow of a charged liquid hydrocarbon have been investigated. 
to the ground were in the range 1 to 6 microamperes. 


The maximum pipe currents 
Depending upon the electrical resistance of 


the pipes, these currents produced potentials ranging from essentially zero to values in excess of 
30,000 volts which were sufficiently severe to cause electrical breakdown and arcs within some of the 


pipes under investigation. 


it is concluded that hazardous pipe potentials, resulting from static elec- 


tricity, can be eliminated in practical applications if the electrical resistance from each and any portion 
of the interior surface of the pipe to the ground does not exceed about 107 ohms. 


Key Words: Breakdown, charged liquid hydrocarbons, charging tendency, charge separation, 
current, filters, inner liner, petroleum liquids, potential, relaxation, resistance, static 
electricity, tetrafluoroethylene tubing. 


1. Introduction 


The petroleum industry has long been concerned 


with the problem of static electricity resulting from 
the flow of relatively nonconducting liquid hydrocar- 
bons through pipes, filters, and other components. 
Considerable literature is available which considers 
methods by which such charges are produced and the 
hazards arising from the introduction of electrically 
charged petroleum liquids into receiving tanks, as 


frequently 
operations. 

Problems related to the electrical currents and po- 
tentials produced in pipes and hoses through which 
charged petroleum liquids flow have not been investi- 
gated thoroughly. Carruthers and Marsh [1]! and 
Carruthers and Wigley [2] considered this problem 
recently and developed equations relating the currents 
and potentials produced with the electrical resistance 
of the pipe, the electrical characteristics of the fuel 
and its velocity. Experimental checks were _per- 
formed using pipes of essentially infinite and zero 
electrical resistance. Although these conditions 
were well suited for investigations of charge relaxa- 
tion, neither condition provided experimental veri- 
fication of the potentials produced in pipes having 
large, but not infinite, values of electrical resistance. 

Under the sponsorship of the Department of the 
Navy. Bureau of Naval Weapons, the National Bureau 
of Standards has conducted an investigation to deter- 
mine the desirable electrical resistance charac- 


occurs during refinery and fueling 


' Figures in brackets indicate the literature references on page 316. 


eristics of hoses used to interconnect the various fuel 
handling components of aircraft and other internal 
combustion engines. It is believed that some of the 
results of the investigation are of sufficient general 
interest for presentation here as they provide further 
experimental verification of the theory and ideas pre- 
sented in [1] and [2] and their Discussions. 

The experimental investigations to be described 
were concerned primarily with a determination of 
the magnitude of electrical currents and voltages 
encountered in pipes of intermediate and very high 
resistivities through which electrically charged liquid 
petroleum flows. The pipes had a nominal inside 
diameter of 0.875-in. and lengths in the range of 20 
to 48 in. Flowrates of 5 to 30 gal per minute were 
used with fuels having conductivities of 8 to 200 
picomho per meter. Pipe currents to the ground 
encompassed the range 10°§ to 6X10 A. Depend- 
ing upon the electrical resistance of the pipes, these 
currents produced potentials ranging from essentially 
zero to values in excess of 30,000 V, which were suf- 
ficiently high to cause electrical breakdown and arcs 
within some of the pipes under test. 


2. Theoretical Considerations 


Consider a pipe of length L, as shown in figure 1, 
through which a relatively nonconducting liquid is 
flowing at a mean axial velocity v. The liquid enter- 
ing the section may or may not contain a significant 
initial charge gi per unit length of pipe which cor- 
responds to a convection current /;= qv entering the 
section under consideration. Likewise the liquid 
leaves the pipe with a charge q, per unit length of 


307 











FiGURE 1. Electric currents resulting from the flow of a fluid. 


pipe causing an exit convection current /,=quv. As 
discussed in [1], the value of g, and hence /, is depend- 
ent upon: the magnitude of Jj, the resistivity of the 
fuel, the residence time, t, during which the fuel re- 
mains within the length of pipe, and the amount of 
charge separation which occurs at the fuel-pipe 
interface. 

The equilibrium current /, =1;—J,, flowing from the 
pipe to the ground, is a direct measure of the net rate 
at which charges are deposited on the interior surface 
of the pipe during passage of the fuel. It is conven- 
iently measured by electrically insulating the pipe 
section from the remainder of the flow circuit and 
from the ground, and electrically connecting each 
end of this section to the ground through a suitable 
ammeter. 

Convenient instrumentation for the direct measure 
of J; is not available. However, when a _ charge- 
producing component such as a microfilter is placed 
immediately upstream of the test section, the current 
I; between the filter and the ground is approximately 
equal to /;. This filter current is conveniently meas- 
ured by an ammeter if the filter is electrically insu- 
lated from the remainder of the flow circuit. Thus, 
I;=1; provided the fuel entering the filter is essen- 
tially uncharged. 


The relation between /, and /;, as developed in the 
appendix, is 


Ip=(i—Ip) 1 —e A) (1) 


where / is the convection current that would emerge 
with the liquid from an infinitely long pipe under the 
same conditions. If /; is caused by a filter placed 
immediately upstream from the test section and is 
large compared to Jo, then the relation 

Ip ~ [41 —e-) (2) 
should be correct for the more severe charging 
conditions. 

The term A is the ratio of resistance time, L/v, of 
the liquid in the pipe to relaxation time, 7, of the 
liquid. Relaxation time is the time required for the 
charge existing in the stationary fuel to be reduced 63 
percent in magnitude by electrostatic forces. 

Considering a pipe grounded at each end, an ap- 
proximate relation for the maximum potential devel- 


oped on the pipe wall, as a result of current /,, is 


Vmax =0.1 Rep, (3) 
where Ry is the end-to-end resistance. Commencing 
with an expression from [1] for the voltage developed 
on the pipe wall at any distance from the inlet, this 
relation has been developed in the appendix subject 
to the following assumptions: 

(a) The ratio A has numerical value less than 6. 

(b) The electrical resistance Rr of the pipe is small 
compared to that of the liquid contained within the 
pipe. 

The constant 0.1 of eq (3) varies through the approxi- 
mate range 0.08 to 0.13 for values of \<6. Values 
of A >6 are encountered only with long pipes, fuel of 
very high conductivities, or very low fuel velocities 
for which conditions nearly all of the separation of 
charges, qi, occurs near the entrance of the pipe. 
Values of A < 6 are those normally encountered during 
the flow of liquid petroleum products, and a constant 
factor of 0.1 can be selected for this range of A with 
little sacrifice in accuracy. 

The assumption that pipe resistance Ry is small 
compared to that of the liquid within the pipe is valid 
for pipes having resistances no greater than 10! Q 
per linear inch of length even when fuels having con- 
ductivities as large as 200 pmho/m are involved. 
Pipes having higher electrical resistances are not 
considered desirable because of the possibility of 
electrical discharges resulting from the extremely 
high voltages produced. Thus, it is believed the as- 
sumptions involved in deriving the approximate 
relation expressed in eq (3) are valid as far as’ practical 
applications are concerned. 

It is of interest to note that eq (3) is essentially 
identical to the relation V,,,=0.106/;—Io)Rr which 
was presented by Dr. Klinkenberg in the Discussion 
of [2]. 

Data were obtained during these experiments which 
verifies eqs (2) and (3) using pipe specimens having 
Ry values in the range zero to 10'°Q. Also, the in- 
fluence of the electrical resistance of the pipe on the 
current J, was investigated; and the electrical break- 
down phenomena which occur when large values of 
I, are encountered in pipes having extremely high 
values of Rr were observed. 


3. Experimental Apparatus 


The flow circuit, shown in figure 2, included a sup- 
ply tank containing 50 gal of naptha and a centrifugal 
pump of 30 gpm capacity. The flow then passed 
through a suitable flowmeter into the microfilter used 
for the generation of the initial charge gi. A by-pass 
valve was installed to provide a method of directing 
the flow through or around the filter as desired. The 
pipe sample under test was connected immediately 
downstream from the filter and the fuel then re- 
turned through a flow control valve to either the supply 
or receiver tank as desired. Each end of the pipe 
sample, as well as the entrance connection to the 
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filter, was electrically insulated from the flow system 
by short lengths of nonconductive tetrafluoroethylene 
(TFE) tubing. Short adapters, of metal tubing, were 
used as connectors between the ends of the test 
pipe and the insulating tubing. They were inserted 
into the ends of the test pipe to a depth of about one 
inch, thus forming internal electrodes. These two 
electrodes were electrically connected to each other 
and to the ground through an ammeter for the meas- 
urement of Jp. Likewise, the filter was connected 
to the ground through another ammeter for the meas- 
urement of /;. All other portions of the flow system 
were constructed from conductive materials and 
grounded as in normal service. 

The liquid used throughout the tests was a water- 
white, hydrocarbon naptha designated as MIL-F- 
7024A, Type II. This naptha has a 60/60 °F specific 
gravity of 0.771, a kinematic viscosity at 80 °F of about 
1.2 centistokes, an initial distillation point above 
300 °F, and a flash point above 100 °F. The electri- 
cal conductivity of the naptha as received, in 55-gal 
drums, at our laboratory is about 1 pmho/m. Con- 
taminants were added to obtain the appreciable filter 
charging effects required. 

Two different filters were used. Each had an ex- 
ternal metal body containing a phenolic resin impreg- 
nated paper element. One filter had a 25-um size 
element of a pleated design containing a total of 35 
pleats l-in. deep by 4.5-in. long giving a total filter 
surface area of about 300 in.22) The nominal flow rating 
of this element was 25 gpm with a fuel of 1 centi- 
stroke viscosity. The second filter had a 10-um size 
element of a stacked-paper washer design having a 
length of 4 in., an outside diameter of 1.5 in., giving 
a total filter surface area of about 18 in.2 This ele- 
ment had a nominal flow rating of 5 gpm with a fuel 
of 1 centistoke viscosity. Both of these filters and 
their elements were identical in size and design to 
those commonly used in fuel systems of some internal 
combustion engines. 

The ammeters used for the measurement of J; and 
I, were multipurpose electrometers having a range 
10°’ to 3 A full scale. The voltage drop of these 
instruments never exceeded 1 V during the measure- 
ment of current. 











SUPPLY 
TANK 


RECEIVER 
TANK 


FIGURE 2. Diagram of the flow circuit. 


Pipe potentials of less than 100 V were measured 
with a multipurpose electrometer having a resistance 
to ground of the order of 10'* Q. Voltages in the 
range 100 to 10,000 V were measured with electro- 
static voltmeters having resistances to ground of the 
order of 10'6 Q. No attempts were made to measure 
potentials exceeding 10,000 V as corona effects with 
the resulting breakdown of resistance were always 
encountered in this region. 

Pin-type probes were constructed for the meas- 
urement of the potentials developed on the inner sur- 
faces of the pipe walls. As will be discussed in detail, 
the test specimens were tetrafluoroethylene (TFE) 
tubing containing an inner liner of intermediate re- 
sistivity and having a total wall thickness of 0.040-in. 
The probes contained pin points about 0.06-in. long 
and were inserted through the tubing from the outside, 
thus puncturing and making electrical contact with 
the inside conductive liner of the pipe. Two pins 
pierced the wall at each probe location. 

Fuel conductivities were determined by a conduc- 
tivity cell constructed in accordance with figure 8 of 
reference [3]. The actual plate separation was 7.9 
mm. A constant test potential of 150 V de was used 
resulting in a test potential gradient of 19 V/mm. 
nificant polarization was obtained during the deter- 
mination of fuel conductivities as evidenced by a 
considerable and continuous decrease in current with 
time under the condition of a constant applied voltage. 
For comparative purposes, the current used to com- 
pute fuel conductivity was that occurring 30 sec after 
applying the test voltage to the cell. Experience 
with this and other methods of conductivity measure- 
ment indicates that the determination of fuel con- 
ductivity is not an exact art. 


Sig- 


4. Test Specimens 


The pipe samples used during the investigation had 
nominal inside diameters of 0.875-in. and their lengths 
were in the range 20 to 48 in. One sample was a 
22-in. length of stainless steel tubing which for pur- 
poses of this investigation is considered to have zero 
electrical resistance. All of the other samples were 
formed from the thermoplastic tetrafluoroethylene 
(TFE) which in a pure state has a volume resistivity 
>10'§ Q-cm. Each sample of TFE pipe had a wall 
thickness of about 0.040-in. The dielectric strength 
of the pure TFE was such that electrical breakdown 
would not occur when a potential of 30,000 V de was 
applied across the wall. 

Experimental samples of intermediate resistivity 
TFE tube were produced for use in this investigation 
by members of the aircraft hose industry. These 
differed from conventional, pure TFE in that each 
contained an inner layer to which small amounts of 
carbon had been added. This resulted in an inner 
liner of moderate conductivity. about 0.010-in. thick. 
extending the full length of the pipe and attached to 
and surrounded by an outer wall of nonconductive 
TFE about 0.030-in. thick. Thus, electrical conduc- 
tion was in the longitudinal direction only. 
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FIGURE 3. End-to-end resistance of two samples of TFE pipe with 
carbon inner liners showing the variation in resistance with 
current. 


Sample A @ Sample B 


Pipes made in this manner have a voltage depend- 
ent, nonlinear, electrical resistance which tends to 
decrease appreciably with increased potential gradient. 
Thus, in performing tests to determine the ohmic 
resistances of the samples, it was found that current 
is a nonlinear function of the applied voltage. Typical 
results of end-to-end resistance tests of two of the pipe 
samples, 36 in. long, are shown in figure 3. These 
resistance values were obtained by applying selected 
constant test potentials across the pipe sections and 
measuring the resulting current. Unless stated other- 
wise, Rr values listed in this paper are those obtained 
with a test current of lwA. 


5. Tests and Results 


5.1. Fuel Contamination and Filter Charging 
Tendency 


Initial tests were conducted to investigate the rela- 
tion between fuel contamination and the charging 


tendency of paper microfilters. The flow circuit, as 
shown in figure 2, utilized a two-tank system to enable 
use of one tank for supply and the other as a receiver. 
This system is desirable for experiments with fuels 
of low conductivities as it ensures uncharged fuel 
entering the pump inlet line. Thus, service operating 
conditions are duplicated closely. In these tests, a 
48-in. length of semiconductive TFE tube (Rr=50MQ) 
was used as the pipe sample, and the filter current, 
I;, and pipe current /, to the ground were measured. 

The naptha had an initial conductivity of about 
1 pmho/m. The two contaminants investigated were 
asphaltenes and a commercial antistatic compound 
(ASC). Fuel conductivities in the range 2 to 200 


pmho/m were obtained by varying the concentration 
of these two materials. The asphaltenes were dis- 
persed quickly and evenly throughout the test naptha 
by first dissolving in isopropylbenzene (Cumene). 
In order to generate large charging currents, these 
initial tests were performed with the 5-gpm stacked 
paper filter at a flowrate of 8 gpm through the filter. 
This corresponds to an average velocity of 4 ft/sec 
across the section of the 0.875 i. d. TFE pipe specimen. 

The results of the stacked-paper filter charging 
tendency tests are shown in figure 4, a plot of filter 
current versus the fuel conductivity. Contamination 
with asphaltenes gave very significant filter charging 
currents with fuel conductivities above 10 pmho/m. 
The maximum charging tendency obtained by contam- 
ination with the antistatic compound (ASC) was only 
about 1/25 of that obtained with the asphaltenes under 
identical test conditions. Also, with the ASC con- 
taminant, charging tendency decreased rapidly to a 
low value as conductivity exceeded 100 pmho/m. One 
additional point on figure 4 shows the effect of adding 
a small amount of the ASC to the asphaltenes mixture 
having an original conductivity of about 200 pmho/m. 
A considerable decrease in charging tendency was 
noted with a relatively small increase in fuel conduc- 
tivity. 

The charging tendency of the napthas containing 
the asphaltenes and ASC contaminants were also 
studied using the apparatus described on page 50 of 
reference [3]. Briefly, this provides a method of 
determining the amount of electric charge produced 
as a fuel flows by gravity through a steel capillary 
having an inside diameter of 3 mm and a length of 
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FIGURE 4. Charging tendency of a stacked-paper filter at a con- 
stant flowrate of 8 gpm showing the influence of fuel conductivity 
with two different fuel contaminants. 
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0.5 m. No relationship was found to exist between 
charging tendencies obtained with the paper filter as 
compared to those obtained using the steel capillary. 
In fact, the ASC contaminant gave approximately the 
same charging tendencies in the steel capillary as 
was obtained with asphaltene contamination. 

Also in tests of other fuels having conductivities 
of 84 and 165 pmho/m caused by unknown contam- 
inants, relatively insignificant charging tendencies 
were indicated by tests with the steel capillary while 
extremely large charging currents were obtained when 
the same fuels flowed through paper microfilters. 

Thus, the investigation of the relation between 
contamination, fuel conductivity, and charging 
tendency showed that the charging tendency of a given 
liquid hydrocarbon is not related to fuel conductivity 
alone. Rather, as might be expected, it is dependent 
upon both the material and the amount of the con- 
taminant in the liquid and upon the wall material at 
the liquid-wall interface. The investigation also 
showed that maximum charging tendency in micro- 
filters will likely occur when fuel conductivity exceeds 
10 pmho/m. 


5.2. Preliminary Tests With Various Pipe Samples 


Preliminary experiments were conducted with var- 
ious pipe samples to determine: (a) the magnitude of 
pipe currents to be expected from pipes installed im- 
mediately downstream from the microfilters; (b) the 
influence of flowrate upon these pipe currents; and 
(c) whether the electrical conductivity of the pipe 
had a significant influence upon these currents. Fuel 
having a conductivity of 165 pmho/m was used through- 
out the preliminary tests. Pipe current, J), was meas- 
ured, but unfortunately, the filter current, J;, was not 
measured during the preliminary tests. However, an 
effort was made to keep the test conditions the same 
as the different pipe specimens were investigated. 
Thus, it can be assumed that the convection currents, 
Ii, entering the different test specimens were approx- 
imately equal at identical flowrate conditions. Only 
the pleated-paper filter was used during these tests. 
The filter by-pass valve in figure 2 was closed; thus all 
flow passed through the filter. 

The pipe specimens selected for the preliminary 
tests all had a length of 22 in. These included a stain- 
less steel tube and two samples of TFE pipe contain- 
ing carbon inner liners having Rr values of 10° and 
10° Q. Voltage measurements were not made and 
the pin-type probes were not applied to the test 
samples. 

The results of the preliminary tests are shown in 
figure 5. It is seen that J, increased significantly 
with flowrate and that the largest values observed 
were in the range 2to6 wA. Also, the pipe currents 
were about equal in magnitude for all values of pipe 
resistance investigated. Thus, it was verified (as 
might be expected) that pipe resistance has very 
little influence, if any, upon the rate of relaxation of 
charges from the charged fuel to the walls of the pipe. 
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FIGURE 5. Current to the ground from pipes of different resistiv- 
ities, as influenced by flowrate (each pipe 22-in. long installed 
downstream from a pleated-paper filter, fuel conductivity 165 
pmho/m). 


5.3. Electrical Breakdown in Nonconductive Pipes 


When a nonconductive TFE pipe sample 22-in. 
long, was inserted in this same flow circuit, electrical 
sparks and other discharges occurring from the sample 
were observed. Also pinhole punctures appeared 
in the walls of the sample through which fuel slowly 
seeped to the atmosphere. The conditions leading 
to the electrical breakdown of the pure TFE pipe 
follow. 

An unpunctured, translucent sample of bare TFE 
pipe 22-in. long, wall thickness 0.040-in., was installed 
downstream from a paper microfilter. The sample 
was electrically insulated from the remainder of the 
flow system by short lengths of TFE tubing. Short 
metal adapters were inserted between each end of 
the sample and the insulating tubing. These two 
adapters were electrically connected to each other 
and to the ground. A probe for measuring electro- 
static potential was installed on the outside at the 
center of the sample. This consisted of a TFE in- 
sulator bushing, 0.25-in. wall thickness, having an 
aluminum foil strip wrapped around the exterior and 
connected to an electrostatic voltmeter. There was 
no other grounded metallic object close to the TFE 
pipe. 

With the room darkened so electrical discharges 
could be seen more readily and using a fuel having a 
conductivity of about 165 pmho/m, flow through the 
system was quickly adjusted to 15 gpm at 25 psig. 
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FiGuURE 6. Electrical discharges in a pure TFE pipe connected 
downstream from a pleated-paper filter (lowrate 15 gpm, fuel 
conductivity 165 pmho|m). 


The electrostatic voltmeter rose rapidly to a maximum 
of about 8 kV and then oscillated about this value 
indicating that some electrical discharges were likely 
occurring. These preliminary discharges were not 
visible in the darkened room. 

After a period ranging from about | to 4 min at this 
constant flow condition, a very bright electrical dis- 
charge took place inside of the pipe. It consisted 
of bright flashes running lengthwise inside of the tube 
with numerous lateral flahses branching off as shown 
in the photograph of figure 6. After this initial dis- 
charge, the voltmeter oscillated about the 5 kV point, 
and the wall of the test pipe contained numerous pin- 
hole punctures through which fuel slowly seeped. 

Light for the photograph of figure 6 was obtained 
entirely from the electrical discharge which took place 
within the translucent pipe sample. The actual por- 
tion of the pipe shown is the upstream part with the 
entrance to the left and the TFE insulator bushing 
of the voltage probe at the right. The scale of the 
picture may be determined from the approximate tube 
diameter of 0.9-in. 

Three pinhole punctures caused by the electrical 
discharge may be noted in figure 6. One is on the 
top of the pipe near the TFE insulator bushing. The 
second is the bright spot in the center of the pipe: and, 
the third is a bright spot at the bottom of the pipe 
towards the entrance end. None of these punctures 
were present prior to the initial breakdown. 

When the TFE pipe was cut open to display the 
inner surface after the initial breakdown, the paths 
over which the numerous spark discharges occurred 
were observed. Rather wide and deep transverse 
and sometimes lateral paths were found to be cen- 
tered about each pinhole puncture. Lighter paths 
branched off from these to cover practically the whole 
interior surface of the pipe with a rectangular net- 
work pattern in which the distance between paths 
was of the order of a few millimeters. The impres- 
sion gained from this network is that of electrical 
charges draining along tributary paths to a central 
stream and thence to a puncture point. Also, it is 
believed that the paths traced by the spark discharges 
are proof that discharges such as are shown in fig- 
ure 6 occur along the inner surface of the pipe and not 
in the liquid itself. 

As flow continued after the initial breakdown, a 
glow discharge and sometimes small sparks were 
noted within the pipe centered about one or more 
of the punctured points. Electricity in the form of a 
purple corona discharge was seen streaming into 


the punctured points from the air. This corona dis- 
charge continued so long as the fuel flowed and its 
effect was to enlarge or erode the pinhole punctures. 

Although not known definitely, it seems reasonable 
that the potentials existing just prior to breakdown 
may have been in the range 50 to 100 kV as this would 
likely be sufficient to cause the dielectric failure of 
the pipe wall. The electrical breakdown of noncon- 
ductive TFE pipe as described above is repeatable. 
It was observed in the same manner with numerous 
pipe samples, and other photographs of these initial 
breakdown flashes have been obtained. 


5.4. Current-Voltage-Resistance Relationship 


Three samples of TFE pipe with carbon inner liners, 
designated as A, B, and C, each having a length of 
36 in., were selected for the investigation of the cur- 
rent, voltage, resistance relationship as expressed 
in eq (3). The end-to-end resistances Rr of samples 
A and B are given in figure 3. The resistances of 
samples A, B, and C were 14,000, 600, and 12 MQ. 
respectively, with a test current of 1 wA. Using the 
flow circuit of figure 2, pipe current, /,, filter current, 
I;, and pipe potential, Vax, were measured. Pin- 
type probes described previously were placed on the 
test samples at distances of 9, 18, and 27 in. from the 
upstream end for the measurement of Vmax. The 
potential between each of these three probes and the 
ground was measured with high impedance volt- 
meters for each flow condition established. 

The hydrocarbon naptha had a conductivity of 84 
pmho/m throughout these tests. Flow rates in the 
range 5 to 26 gpm were selected. These rates, in 
combination with the particular fuel conductivity and 
the inside pipe diameter of 0.875-in., resulted in values 
of A in the range 1 to 6. Pipe currents, Jp, ranging 
from 2X 10-" to 3X10-® A could be obtained con- 
veniently by conducting tests with the filter by-pass 
valve both open and closed over the range of flowrates 
selected. é 

The results of these tests are summarized in figure 
7, a log-log plot of the maximum pipe potential, Vmax. 
determined from measurements at the three probes, 
as a function of the current, /,, from the pipe to the 
ground. All of the pipe currents in figure 7 having a 
magnitude less than 10°7 A were negative. These 
were obtained with the filter by-pass valve open, a 
condition of minimum filter charging. It is believed 
these smaller currents all resulted from the separation 
of charges at the liquid-pipe interface. The pipe 
potentials corresponding to these small negative cur- 
rents were also negative with the single exception of 
sample A at 510°’ A where positive readings of 
voltage were obtained. No explanation can be sug- 
gested for this single deviation, but it is of interest 
to note this test point does not lie on the straight line 
defined by the other points. 

All J, currents greater than 10-7 A were positive and 
were obtained with flow through the filter, i.e., the 
by-pass valve completely closed. These larger cur- 
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FIGURE 7. Maximum potential versus pipe current in pipes with 
carbon inner liners (length 36 in., flowrate range 5 to 26 gpm, fuel 
conductivity 84 pmho|m). 


rents were the result of the separation from the fuel 
to the wall of the pipe of charges produced during 
flow through the filter. Also, all potentials resulting 
from the positive currents were likewise positive. 

Still referring to figure 7, it is seen that one straight 
line for each pipe sample defines the test points quite 
closely and the slopes of these three lines are approxi- 
mately equal. Also, the slopes are considerably less 
than 45°, indicating that the resistances of the three 
pipe samples were not constant, but tended to decrease 
with increasing voltage. 

The experimentally determined maximum pipe 
potentials of figure 7 can be compared to those pre- 
dicted by eq (3) provided Rr values corresponding 
to the pipe currents under consideration are used. 
Figure 3 gives these values for samples A and B. In 
selecting a value of Rr from this figure which cor- 
responds to a designated J, current, it is assumed 
that the net deposition of charges is relatively uni- 
form lengthwise along the pipe. Thus, the average 
current existing in the wall of the pipe can be con- 
sidered to be approximately 0.25 J,. This is the 
current which has been used in obtaining values of 
Rr from figure 3. 

The comparison of calculated and measured maxi- 
mum pipe potentials for samples A and B is presented 
in table 1. Considering the practical difficulties 


TABLE I. 
voltages 


Maximum pipe voltage 
R 
(hg. 3) 


Sample | 
Calculated Measured 
Vmax = O.1R rl, (hig. 7) 


Amperes Ohms Volts | Volts 
10 x 1H 10 | 700 
Ca | x 10 2.200 | 2.700 
10 2x 10" 12,000 11.000 


10-6 x 10" | 90 | 100 
10-5 7x 10° | 470 490 


Comparison of calculated and measured maximum pipe 


involved in measuring Ry and Va, accurately it is 
believed that good agreement between the calculated 
and measured values was obtained. Certainly agree- 
ment such as this is quite adequate for the practical 
purpose of estimating potentials to be encountered 
in pipes through which a known current /, to the 
ground is encountered. However, it must be realized 
that eq (3) applies only to the conditions in which the 
conductive portion of the pipe is in actual contact 
with the fuel throughout its length and the pipe is in 
electrical contact with the ground only at the two ends. 
Grounding of the conductive portion of the pipe at 
many locations throughout the length will of course 
result in lower voltages. Insulating the conductive 
portion of the pipe from the fuel will likely result in 
potentials having magnitudes sufficient to 
electrical breakdown. 


cause 


5.5. Charge Relaxation in Short Pipes 


Both filter current, J; and pipe current. /,. were 
measured during the experiments reported in figures 
4 and 7. During that portion of these tests in which 
the fuel passed through the filter. the values of /; always 
exceeded 0.5 uA. These are very large compared 
to the currents generated at the fuel-pipe interface 
when the fuel bypassed the filter. Thus. for flow 
through the filter, it can be assumed that the measured 
I, values are equivalent to the charges which migrated 
to the pipe wall during passage of the charged fuel 
through the pipe. Thus, data are available to verify 
eq (2). 1, ~ [41 —e-). 

The results of this check of the relation /, ~/4 -e~) 
are shown in figure 8. Each plotted point is the ratio 





} 
LENGTH} 
IN 
36 
8 
@ 6.0x10 36 
X 12x10’ 36 
D 5.0x107 48 








4 


FIGURE 8. Experimental verification of the relation \,/1,=(1 —e>*) 
for various pipe samples, flow rates, and fuel conductivities. 
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I,/l; as calculated from the measured values of these 
two currents. The solid line represents the function 
(1 —e~) versus A, computed from the relation A= L/vr. 
Data corresponding to the three pipe samples 36-in. 
long were all obtained with a constant fuel conductivity 
of 84 pmho/m, and flowrate in the range 5 to 26 gpm. 
Data for the one pipe sample 48-in. long were obtained 
at a constant flow of 8 gpm with fuel conductivities 
in the range 11 to 200 pmho/m using asphaltenes as 
the contaminant. The few experimental points ob- 
iained for the A range 5 to 10 are not included in this 
figure because these measured values of /,/l; were 
all within a few percent of unity. It is believed that 
the agreement shown in figure 8 is reasonable over the 
range of conditions encountered. However, some of 
the points indicate relaxation may at times occur more 
rapidly than predicted. Certainly these results help 
to verify that relaxation of a charged liquid flowing 
in a pipe is essentially independent of the electrical 
resistance of the pipe for resistance values in the range 
zero to about 10° Q. 

It is believed, but not definitely known, that at least 
some of the deviations shown in figure 8 are related 
to the determination of fuel conductivity. Perhaps 
the effective conductivity of the charged fuel in the 
pipe can be influenced by the potential gradients 
existing. Thus, it is possible that present methods 
of determining fuel conductivity are not completely 
adequate for the experimental verification of the 
theoretical relations. 


6. Discussion 


The purpose of this section is to discuss and further 
emphasize those results leading to the desirable resist- 
ance characterisitcs of pipes and hoses through which 
charged liquid petroleum products may flow. It is 
believed that the experimental conditions of flowrate, 
and of pipe diameter, length, and installation down- 
stream from microfilters duplicate closely conditions 
frequently encountered, especially in fuel systems 
of the larger size internal combustion engines. Ad- 
mittedly, this investigation was concerned primarily 
with the more severe conditions rather than with 
average conditions. Also, the fuel conductivities of 
84 and 165 pmho/m used during many of the tests 
are higher than those normally encountered. How- 
ever, experiments with fuels of lower conductivities 
show that maximum filter charging effects usually 
occur in the conductivity range 10 to 200 pmho/m 
depending upon the filter design and the rate of flow. 
Experience has shown that conductivities in the range 
10 to 30 pmho/m are encountered rather frequently 
in liquid petroleum fuels flowing into the engines. 

Perhaps one point should be emphasized concern- 
ing fuel conductivity. In investigations concerned 
with the more severe filter charging conditions and the 
resultant currents and potentials in pipes installed 
downstream, it has been found, as shown in figure 4, 
that greatest electrification effects are obtained with 
fuels of moderate conductivities, range 10 to 200 
pmho/m, rather than those having low conductivities 


of the order of one pmho/m. Therefore, elaborate 
precautions to remove absorbed water and as much 
of the contamination as possible are not necessary 
during investigations involving the generation of large 
filter charging currents. Nevertheless, in this con- 
nection some effects of water should be mentioned. 
As discussed in [3] and investigated briefly by the 
authors, ordinary tap water, which has relatively very 
high conductivity, does not increase fuel conductivity 
even when present to the extent of 5 percent or so. 
In fact, water frequently tends to leach out and reduce 
the amount of contaminant in a fuel with the resultant 
effect of a possibly large decrease in conductivity. 
Finally, water mixed with fuel, especially when dis- 
persed in the form of droplets, usually increases elec- 
trification considerably during flow through pipes and 
filters. 

It has been shown that electrical breakdown can 
occur when a charged fuel flows through nonconductive 
pipes. These phenomena are perhaps similar to 
the charging of an electrical condenser until the result- 
ing potential gradients are sufficient to cause break- 
down. Thereafter, depending upon the current and 
the particular materials involved, either a continuous 
discharge will be obtained through the breakdown 
path, or a series of intermittent discharges will result. 

The question arises as to what happens when the 
nonconductive pipe has a sufficiently high dielectric 
strength so failure through the wall does not occur 
at the time of the arc discharge. An example of such 
conditions is described below. 

Variable-area meters are used frequently in the 
laboratory for the flowrate measurement of hydro- 
carbon liquids. The metering tubes of these instru- 
ments are moulded of thick heavy glass, have a length 
in the range | to 3 ft, and frequently have no internal 
bonding wire or its equivalent. Occasionally such 
meters are connected downstream from microfilters 
which charge the fuel significantly under certain op- 
erating conditions. 

It is rather well known that electrical discharges 
occur at times in glass-tube flowmeters so installed. 
Photographs of such discharges are very similar in 
appearance to figure 6 except that puncture or break- 
age of the glass pipe does not occur. Rather, the 
strong lengthwise spark terminates at any convenient 
metallic ground within the end of the glass tube. 
Numerous lateral flashes are noted also. Although 
not definitely known, it appears from photographs, that 
the electric arcs occur on the internal surfaces of the 
glass pipe. Arcs extending five or more inches up 
into the glass tubes have been observed under the 
more severe conditions. 

The use of an internal bonding or ground wire in 
nonconducting pipes should reduce maximum _ po- 
tentials significantly and thus eliminate the possible 
hazard arising from electrical arc discharges. This 
internal bonding wire will likely have little, if any, 
influence upon the rate at which separation of charge 
from the highly charged fuel occurs. 

The possibility of placing a conductive skin around 
the outside of nonconductive hoses was mentioned 
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in the Discussion of reference [2]. This idea may be 
practical provided the transverse leakage of elec- 
tricity from the interior of the hose to the conducting 
skin does not puncture the wall. Experience has 
demonstrated this is not practical when applied to 
hoses containing a thin tube of pure TFE surrounded 
by an exterior, grounded, metal braid. In such hoses, 
discharges of electricity similar to that of figure 6, 
frequently puncture the TFE tube resulting in liquid 
leakage through the wall. One solution to the prob- 
lem is the development of a practical inner liner or 
coating of moderate resistivity which would be 
grounded at each end of the hose. 

Considering now the flow of a highly charged liquid 
hydrocarbon through a pipe of moderate resistivity. 
the experimental results help to verify that the rate 
of charge separation from the fuel to the pipe wall is 
essentially independent of the pipe material or of its 
electrical resistance. The relation /,=/{1—e~*) is 
in reasonable agreement with the experimental results 
provided J; is large compared to Jo, where [o is the 
convection current flowing out with the fuel from an 
infinitely long pipe of the same material. 

The maximum magnitudes of the currents /, en- 
countered in the investigation were in the range | to 
6 uA. These are of the same order of magnitude and 
are in reasonable agreement with a possible maximum 
value of 10 wA cited at times in the literature. Ifa 
safety factor is desired, a maximum value of /,= 100 
pA will perhaps cover all operating conditions which 
could be encountered in practical applications. 

The rather limited experiments indicate that the 
maximum voltages occurring in pipes can be predicted 
with reasonable accuracy by the relation V max = 0.1 R rl, 
provided A<6. According to eq (5) in the appendix 
and as shown in figure 9, for any given value of Rr 
and /,, maximum potentials will be less for values of 
\ considerably greater than 6. However, such con- 
ditions are not usually encountered in practical appli- 
cations. If potentials existing within pipes are not 
to exceed an arbitrary value of 1,000 V for a possible 
maximum pipe current of 100 wA, then the end-to-end 
resistance of the pipe may be computed as Rr= 1,000/ 
(0.1 X 100 x 10-8) = 108 . 

Thus. a value of Ry not to exceed about 108 may 
be suggested for pipes and hoses intended for use 
with liquid hydrocarbons. It is believed this upper 
limit is sufficiently low to eliminate the hazard of spark 
discharges occurring either within or through the 
walls of hoses and pipes through which a liquid hydro- 
carbon flows. However, this suggested maximum 
value for Rr is subject to three conditions: (a) the pipe 
must be uniformly conductive throughout its length; 
(b) the conductive portion of the pipe must be in con- 
tact with the fuel throughout the length of the pipe: 
and (c) the conductive portion of the pipe must be 
electrically connected to the ground at each end of the 
pipe. 

When conventional end fittings of aircraft hose are 
used as electrodes, the end-to-end resistance R; is 
not a reliable parameter for specification purposes 
because it does not assure that the interior surface 


of the pipe is conductive and electrically connected 
to the ground. However, this requirement will be 
checked if the resistance from each and any portion 
of the interior surface of the pipe to the ground is 
measured instead of the end-to-end resistance. Con- 
sider an empty pipe grounded at each end. An in- 
ternal electrode located at the center of this pipe has 
the effect of placing two Ry7/2 resistances in parallel 
resulting in an effective resistance to the ground of 
R7/4. With allowance for this, and considering orders 
of magnitude only, it is concluded that hazardous 
pipe potentials resulting from static electricity can 
be eliminated if the electrical resistance from each and 
any portion of the interior surface of the pipe to the 
ground does not exceed about 107 Q. The use of an 
internal electrode in measuring the resistance of pipes 
and hoses will be inconvenient at times. but this 
method does provide an electrical path to the ground 
similar to that seen by the liquids contained within 
the pipe during normal usage. 

It is acknowledged without reservation that the 
relations and ideas presented in [1] and [2] and their 
Discussions provided nearly all of the background 
material required for this investigation. It merely 
remained necessary to demonstrate that such rela- 
tions and ideas were applicable to the particular 
problem under consideration. Thus, it is gratifying 
to be able to present additional experimental data 
which, for practical purposes, is in good agreement 
with the work of previous investigators. 


7. Conclusions 


1. High electrical charge densities can be produced 
in petroleum liquids having moderate to high electrical 
conductivities (10 or more pmho/m) when such liquids 


flow through conventional microfilters. However. the 
charging tendency of the liquid is not dependent upon 
fuel conductivity alone. Rather. it is dependent upon 
the amount and material of which the contaminant is 
composed and the material at the liquid-solid interface. 

2. Maximum currents in the range one to 6 wA were 
obtained from the pipe to the ground when highly 
charged petroleum liquids flowed through short pipes 
or hoses. These comparatively large currents were 
caused by the migration to the wall of the pipe of the 
electrical charges contained within the charged fuel. 

3. The rate of separation of electrical charges from 
the charged fuel is essentially independent of the pipe 
material and its electrical resistance. 

4. When the rate of charge separation is sufficiently 
high in pipes constructed from nonconductive ma- 
terials, electrical discharges in the form of sparks 
and ares will occur. These discharges will originate 
on the inner surface of the pipe. They may be con- 
tained entirely within the pipe or may cause puncture 
or other failure of the pipe wall depending upon the 
pipe material, the wall thickness, and the ground 
locations involved. 

5. Maximum potentials produced by static electric 
phenomena within pipes, grounded at each end, may 
be predicted with reasonable accuracy by a simple 
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relation provided the electrical resistance of the pipe 
and its current to the ground are known. 

6. Considering pipes or hoses through which a 
highly charged petroleum liquid may flow, an elec- 
trical resistance from each and any portion of the 
interior surface of the pipe to the ground not exceeding 
about 10° © is adequate to limit pipe potentials result- 
ing from static electricity to less than 1,000 V. 
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9. Appendix 


The object of this appendix is to develop an approxi- 
mate relation for the maximum electrical potentials 
developed on the inner wall surface of pipes as a re- 
sult of static electric phenomena caused by the flow 
of hydrocarbon liquids. Essentially this theory is a 
simplification of that developed by Carruthers and 
Marsh [1] subject to the assumptions that the elec- 
trical resistance of the pipe is low compared to that 
of the liquid and that the ratio of the residence time 
of the liquid in the pipe to the relaxation time of the 
liquid has a value less than 6. Both of these assump- 
tions appear to be valid as applied to the more severe 
conditions likely to be encountered in practical 
applications. 


I;= current between the ground and the filter 

1;= qiv=convection current entering the pipe due 
to the flow of liquid containing an initial 
charge qi per unit length of pipe 

/,= current flowing from the pipe to the ground 

1)>= convection current which would emerge with 
the liquid from an infinitely long pipe 

L=length of pipe 

m=x/L 

R,=L/B=end-to-end resistance of the pipe 

Vm=the potential on the inner wall of the pipe at a 
position m= x/L from the pipe entrance 


V max =the maximum potential 
wall of the pipe 
V,=the potential on the inner wall of the pipe at a 

distance x from the pipe entrance 

v= mean axial velocity of the liquid 

x= distance along pipe from the entrance 

8B=conductance per unit length of pipe wall as- 
suming the conductivity of the pipe wall is 
considerably greater than that of the liquid 

A=L/vr=ratio of the residence time for the liquid 
in the pipe to the relaxation time of the liquid 

7=relaxation time of the liquid. Its value is com- 
puted from eeo/k where k is liquid conduc- 
tivity, e the dielectric constant of the liquid, 
and eo the electric constant of free space. 
When & is expressed in mho/m, 7 in seconds, 
and e is taken as 2, T=17.7 X 10-!2/k. 


existing on the inner 


Carruthers and Marsh [1] in a consideration of the 
potentials produced by the flow of hydrocarbon liquids 
have developed the expression 


> UT = oni x —Livr 
.=(1;—I) = er = —( ee (1) 
V-=([i— 1) B (l—é« ) L" e ] 


4 


for the potential at any point of the cross section of 
liquid and pipe. 

If it is assumed that the electrical resistance of the 
pipe is small compared to that of the liquid contained 
within the pipe, Rr=L/B, and by letting m=x/L and 
A=L/vr, eq (1) may be written as 


Vn = (Uo) SE ((—e-™)— m( —e-) (2) 


Also in [1], the expression 


ly = lie A+ ](1—e A) (3) 


is developed for the convection current, /,, emerging 
with the liquid when the pipe is of finite length. This 
assumes both the relaxation of /; and the generation 
of Io at the fuel-pipe interface are exponential proces- 
ses. As shown in the discussion of figure 1, /,=/i—TIp. 
and from eq (3) the currents /,, /;, and Jo are related 
by the expression 


1, = (Ui — To) A — e~). (4) 


Substituting for (/;— Jo) in eq (2) another expression 


— R; (1 —e-™) 
Vale [Gmc r) 


~m|=KIpRs (5) 


is obtained for the potential at any point along the pipe. 

Figure 9 shows K plotted as a function of A for 
selected values of m. It is seen that K varies between 
0.13 and 0.08 for values of A less than6. As A becomes 
greater than 6, K decreases rapidly becoming less 
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FIGURE 9. The relation K=> doen mm) asa function of 
ae 


for selected values of m. 


than 0.02 for A=50. Values of A exceeding 6 are 
encountered only with long pipes, fuels of very high 
conductivities, or very low fuel velocities: conditions 
in which nearly all of the charge relaxation occurs in a 
short upstream portion of the pipe. Values of A less 
than 6 correspond to the more severe conditions en- 
countered in practical applications. Therefore, with 


little sacrifice in accuracy an average value of K=0.1 
can be used to obtain an approximate expression for the 
maximum potential produced in a pipe, and 


J ‘max =0.1 Ril, (6) 


where Rr is the end-to-end resistance of a pipe, 
grounded at each end, from which a current /» is 
flowing to the ground. 
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A study of the aging behavior of low-power germanium alloy switching transistors has revealed 
a relationship between small changes in junction leakage current, induced by a brief aging stress, and 
later deterioration in performance. This relationship may provide the basis for a nondestructive 
screening procedure which would serve to identify germanium alloy transistors likely to deteriorate 
through excessive growth of junction leakage current. The proposed screening procedure irvolves 
the determination of relatively small changes in junction leakage current, increases of the order of 
15 percent or more, associated with 1,000 hours of aging at a shelf (bake) stress of 100°C. Because 
the leakage current changes of interest are small, relatively high demands are placed upon measure- 
ment repeatability. There is evidence that transistors bearing identical type numbers but of different 
manufacture respond differently to the same screening procedures and would therefore require dif- 


ferent screening limits. 


Key Words: Aging behavior, alloy transistors, failure prediction, leakage current, screening pro- 
cedures, temperature stress, transistor, transistor failure, transistor screening. 


1. Summary 


The data from a 20,000-hr transistor aging behavior 
study indicate that small changes in junction leakage 
current, induced by baking, provide a means for iden- 
tification of units likely to deteriorate through growth 
of leakage current. A subsequent independent ex- 
periment, which included transistors of different manu- 
facture, provided further evidence that this behavior 
may be generally characteristic of germanium alloy 
transistors. The second experiment also examined 
the effects of more stringent screening stresses. 


2. Introduction 


Studies of aging behavior of transistors are never 
free from many of the confounding effects which tend 
to degrade measurement repeatability. The data 
from recent experiments indicate that, generally, 
junction leakage current tends to decrease during 
early stages of aging. As aging continues, or as aging 
stresses are increased, the downward trend ultimately 
reverses. For some transistors, junction leakage 
currents were observed to increase from the outset. 
Those units exhibiting the largest increase formed a 
group which contained all (or nearly all, depending 
upon performance criteria) of the transistors which 
deteriorated more rapidly as aging continued. It 
appears that for these transistors, early increases in 
leakage current are a manifestation of a continuous 
process during which leakage current continues to 
increase, leading eventually to “failure.” Since leak- 
age-current degradation is not uncommon, careful 
observation of early changes in leakage current shows 
promise as a basis for a practical means of screening 
transistors to further improve transistor reliability. 


The fractional changes in leakage current of interest 
are relatively small, and therefore all sources of 
variability in the measurement process are of concern. 
The inherent behavior of the device itself is one such 
source, and in some instances is the major source 
of variability. The measurement of transistor param- 
eter values is not only dependent upon bias conditions 
and ambient temperature, but also upon the time dura- 
tion of the measurement and upon the bias and 
temperature histories resulting from prior measure- 
ment and aging stresses. For example, voltage 
breakdown measurements using current-limited values 
even as low as 10 wA have been observed to commonly 
cause marked disturbances in the leakage current 
values, requiring several weeks for the effects to sub- 
side. Voltage breakdown measurements should be 
avoided wherever sensitive leakage current measure- 
ments are being made. As another example, the 
transfer of a transistor from an aging to a testing 
environment is accompanied by a shift in parameter 
values which usually requires several days to approach 
a new center value. This stabilization effect ' appears 
repeatable for a given transistor but is not consistent 
among transistors. The above mentioned disturbances 
cannot, in general, be eliminated, but their effects 
may be substantially reduced by careful choice of 
operating procedures and test methods. 

Similarly, the instability of measurement equipment, 
especially long-term drift, as well as variation in am- 
bient temperature during measurement, introduces 
errors which cloud the meaningfulness of aging be- 
havior data. 

The National Bureau of Standards in cooperation 
with the Department of the Navy is conducting 


' Zierdt, C. H., Jr., On importance of operating tests as compared to storage tests of 
transistors, The Solid State Journal 2, No. 9, 21-27 (1961) 
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transistor aging studies in which every reasonable 
effort is being applied towards minimization of un- 
wanted variability in the aging data through careful 
instrumentation, selection of test methods and _ pro- 
cedures, and design of experiments. A key feature 
of the instrumentation for measurement of d-c param- 
eters, for example, is close ambient-temperature 
control. The transistors are immersed in a “‘con- 
stant-temperature” recirculating air stream. The 
air stream temperature is centered at 25 °C and fluc- 
tuations are regulated within +0.01 °C. Further 
precautions include replication of measurements 
coupled with computer analyses programmed to detect 
and identify questionable agreement between repli- 
cations. This enables us to repeat questionable 
measurements, if necessary. 

These and other precautions have helped to main- 
tain a high level of equipment performance and 
measurement repeatability. A case in point is the 
record of temperature-sensitive leakage current 
measurements accumulated during the course of a 
three-year aging study wherein each of 10 “control” 
transistors was measured a total of over 3,000 times. 
Among these measurements are four sets of junction 
leakage current measurements, each set consisting 
of 300 recordings. These recordings were found to 
have gradually decreased approximately 14/2 to 2 
percent during the 3-year period. The major part of 
this decrease is traceable to creep of 0.08 °C in the cali- 
bration point of a thermometer due to stabilization 
of the glass bulb. Corrections for the gradual change 
in ambient temperature between each measurement 
period resulted in an estimated coefficient of variation 
of less than 1 percent. This performance record is 
consistent throughout the data and demonstrates 
the stability of the transistor-measurement system 
combination and serves as a measure of confidence 
which may be placed in the transistor aging data dis- 
cussed in this paper. 


3. First Experiment 


The first of two experiments described in this paper 
provided a survey of the behavior of approximately 
500 type 2N396 germanium alloy transistors in re- 
sponse to a wide range of aging stresses. The second 
experiment used approximately 100 transistors to 
explore more fully behavior in the 100 to 125 °C region 


and to introduce different types and manufacture of 


germanium alloy transistors into the study. All 
transistors had met stringent military requirements 
and were considered “high reliability” units. 

A full reporting of the first experiment, in which 
over Ya million measurements have been recorded, 
would be prehibitively lengthy. A brief description 
of the entire experiment, however, is given in the 
appendix. Attention is focused on that part of the 
experiment involving the four groups of transistors 
which were subjected to four different combinations 
of aging. Each combination was calculated to produce 
a junction temperature of 100 °C (identified as aging 


conditions 13, 14, 15, and 18 in table 3). This tem- 


perature lies roughly in the center of the aging-stress 
spectrum for the total experiment. 

The 100 °C aging conditions are of particular interest 
because the small but significant changes in leakage 
current, occurring within the first 1,000 hr of aging, 
take on meaningful significance when compared with 
later aging behavior. 

The aging conditions at 55 °C and below produced 
too little deterioration in 20,000 hr of aging to be of 
much value in testing relationships between early 
change and aging. Comparison of results for the 
different aging temperatures, however, were not in- 
consistent with the expectation that failure rates (and 
chemical activity) double for each approximately 10 °C 
increase in junction temperature. The volume of data 
from the lower stress conditions, 55 °C and below, pro- 
vided assurance that the measurement system, as well 
as the transistors, was quite stable and factors affecting 
repeatability were in close control. 

The aging conditions at 145 °C and above, on the 
other hand, caused relatively rapid deterioration, and 
in a few cases, total destruction occurrell prior to the 
first measurement period, i.e., during the first 340 
hr of aging. Further aging showed conclusively that 
the failure rates were too high to expect the data for 
these conditions to shed much light on the meaningful- 
ness of small early changes. 

For the reasons given above, the remainder of the 
discussion of the first experiment will be devoted to 
the behavior of the four groups consisting of 160 tran- 
sistors aged with the junction temperature maintained 
at 100 °C for a total of 20,000 hr. 


4. Data 


Each of the four groups of 40 transistors was sub- 
jected to a different combination of ambient tempera- 
ture and applied bias power. The four aging conditions 
are identified in table 3 of the appendix. Three of 
the conditions were essentially shelf aging at an 
ambient temperature of 100 °C; condition number 13 
provided no bias power, condition 15 included a 
nominal 3 mW of power dissipation, and condition 18 
subjected the collector-to-base junction to a reverse 
bias voltage stress of —24 V. The fourth aging con- 
dition provided a combination of 70 °C ambient tem- 
perature and 120 mW of power dissipation which was 
the power required to raise the junction temperature 
to an estimated 100 °C. | 

The general effect of the four 100 °C 


aging condi- 
tions was to cause only a very small and gradual in- 
crease in the dispersion of the electrical parameters 


with no perceptible trend in the modal values. For 
those transistors which eventually exhibited substan- 
tial degradation, the changes were mainly gradual and 
thus presented the problem as to when performance 
should be judged “unsatisfactory.” Since no _par- 
ticular circuit applications with associated perform- 
ance limits were involved, the basis for such a decision 
is necessarily arbitrary. The criteria selected are 
similar to those appearing in the specification for the 
environmental acceptance tests, IEBO (—10 V) or 
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ICBO (— 10 V) greater than 10 wA or hFE (10 mA, 
—5 V) less than 20. Using these criteria, the per- 
formance of a total of 10 transistors became “unsat- 
isfactory” during the 20,000 hr of aging at 100 °C. 
(The total elapsed time of the experiment, including 
the measurement periods, exceeds 3 years.) 

The records of behavior of these 10 transistors have 
been scrutinized and a summary is given in table 1. 
On the left, the individual transistors are identified 
along with their respective aging conditions. The 


TABLE 1. 


0-hr values of Early changes 


T 
| 


Ratio values of IEBO or ICBO, 


whichever is larger 


| 


Transistor | Aging 


and plug No.| condition IEBO ICBO 


340-hr 1,000-hr 1,000-hr 


340-hr 
1.061 
2.182 
1.426 
1.054 
1.342 
1.024 
/1.049 
1.049 
1.005 
1.265 


1058 
1066 
1040 
1052 
1064 


1.645 
1.152 
1.612 
1.496 
1.489 
1050 "2.873 
146 1.206 
1068 | 7854 | 
1070 17 
1.215 | 


3.578 
1/1.039 
1.218 
3.658 
1/1.106 
1.041 
11.013 


2.101 
1.737 
3.854 
1.213 
1.066 
1/1.062 
1.064 
1/1.082 

2.190 


wd O-hr | 0-hr | 
57 | 3.795 | 
| 
| 
| 
| 


1O70 
i | 
‘Case leak detection test was applied after 20 khr to all of above transistors except 


1050 and 8-- 1046. Only 1— 1070 leaked: rate was ~ 200 * 10-* 
Variable 


“Only 0-hr measurement was variable 


std. cm/s 


(') Wafer was found to be cracked 


@ 1 arge and abrupt increase in emitter current 


“Indicates emitter and “C™ collector junction 


Returning to table 1, the initial or 0-hr IEBO (— 10 
V) and ICBO (— 10 V) values listed on the left were ob- 
tained before application of the aging stresses. The 
distribution of these values is not perceptibly different 
from that of the larger group of “satisfactory” tran- 
sistors. Even the single maximum value listed would 
have been ineffective for screening, because at least 
10 “satisfactory” transistors had larger values. We 
conclude therefore that these 0-hr values provide no 
effective means for identification of the “‘unsatis- 
factory” transistors. 

Variability of the parameter values during measure- 
ment was another property of interest. A few of the 
0-hr measurements in table 1 bear superscripts in- 
dicating that these measurements were classified as 
variable on the basis that the automatic digital volt- 
meter, which has a 2- to 3-sec conversion time, would 
not balance within 50 sec. This variability test was 
applied uniformly to all transistors throughout the 
experiment. There is no evidence that 0-hr measure- 
ment variability is uniquely related to deterioration. 
Over an extended period of aging, however, the time- 
series measurements indicated that consistently 
variable units exhibited more pronounced deteriora- 
tion. The association was not so close (nor the 
technique amenable) as to favor the use of this cri- 
terion to screen potential failures. 


i i 


code numbers for the aging conditions are fully de- 
scribed in table 3 of the appendix. “Unsatisfactory” 
transistors appeared in all four conditions, and dif- 
ferences in the number of “unsatisfactory” units 
corresponding to each condition are not sufficiently 
large to be considered significant. For these reasons, 
and for the lack of any other substantial differences 
in observed transistor behavior between the four aging 
conditions, it was decided to combine the data for the 
subsequent analyses. 


Summary of aging behavior of “unsatisfactory” transistors 


T 


| 


Deterioration 


“Unsatisfactory” Character of 


Hours of deterioration 


aging 


criteria 


IEBO ICBO Ohmic 


leakage 


Avalanche | Other 


| 
| 


In contrast to the ineffectiveness of the 0-hr leakage 
current measurements, changes in the value of junction 
leakage current observable within the first 1,000 hr 
of accelerated aging appear promising as a screening 
index. None of the other parameters in the experi- 
ment, measured or computed as listed in the appendix, 
appears to be as useful. For this reason the remainder 
of table 1 provides only for comparison between frac- 
tional changes in junction leakage current measure- 
ments and deterioration. Ratio values greater than 
unity indicate an increase in leakage current with 
time. 

The ratio values given are the ratio of the 340-hr 
measurement divided by the 0-hr measurement, and 
the ratio of the 1,000-hr measurement divided by the 
0-hr measurement. In each instance the ratio in- 
dicates the fractional change in junction leakage 
current; either the IEBO or ICBO change is listed, 
whichever is the larger. By thus combining the larger 
change into one table, a closer relationship between 
early fractional changes and deterioration in either 
junction is more apparent than if IEBO or ICBO were 
treated separately. 

Although table 1 does not indicate which junction 
current is associated with the larger ratio value listed, 
the data reveal that with one exception, the junction 
associated with the larger ratio value was also the 
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junction which deteriorated and led to “‘unsatisfac- 
tory” performance. The one exception was transistor 
0-1040, where “unsatisfactory” was indicated by 
hFE* becoming less than 20, first observed after 3.4 
khr of aging. In the very next measurement period 
(4.7 khr), the ICBO value exceeded 10 wA, and it was 
the collector junction which exhibited the larger ratio 
value. The consistent correspondence between 


early changes and deterioration in the same junction 
argues for the existence of a determinative relationship. 


5. Character of Deterioration in 
Performance 


The “Hours of Aging” in table 1 indicate when the 
transistor was first observed not to have met perform- 
ance criteria, and was thereby adjudged “unsatisfac- 
tory.” In each instance there was an increase in 
leakage current. The nature of these increases 
suggested the three categories listed under “Character 
of Deterioration.” Avalanche conduction was ascribed 
to two transistors which exhibited a several-fold in- 
crease in (emitter) junction leakage current at 10 \ 
with essentially no associated increase at 2 V. Con- 
duction by current multiplication was clearly present. 
The second, and most common form of deterioration, 
was a gradual increase in leakage current for both bias 
levels. Six transistors displayed this type of deteriora- 
tion. For each measurement period the estimated 
intrinsic saturation current, which was found by 
extrapolation of the current versus voltage curve to 
estimate the saturation current at 0 V, remained rela- 
tively constant, thus indicating the presence of a 
linear (ohmic) relationship between bias voltage and 
observed leakage current. The third characteristic 
of deterioration, designated “other,” applies to two 
transistors which became catastrophic failures. The 
first of these two transistors, number 0—-1050, de- 
veloped a collector-to-base and _ collector-to-emitter 
“short,” which was apparently caused by a fractured 
wafer found by optical examination. The presence of 
the overt fault was observed after this transistor, as 
well as all other transistors under study, was re- 
mounted and soldered to a new type of high-tempera- 
ture plug. The implied causal relationship, damage 
by remounting, is not so clear-cut as it might appear 
in light of the following considerations. First, the 
soldering and handling techniques had been checked 
in a preliminary experiment in which the effects of 
the same soldering techniques were found to produce 
no significant change (less than 1 percent) in leakage 
current. All transistors were monitored before and 
after remounting, and only 0—1050 showed a change. 
Second, earlier behavior of 0-1050 was suspect. The 
initial 2-V ICBO value was 12 wA, and the 10-V value 
was, oddly, only 3 wA. Both measurements were 
variable. Replicate ICBO measurements made later 
in the same day were also variable, but were near 
1 pA. The initial higher values did not repeat. 
Measurements of hFE were variable and large. They 
became normal after 340 hr of aging. In fact, from the 


340-hr measurement period until the time it was with- 
drawn from the experiment at 14 khr, transistor 
0-1050 was one of the best-behaved units in its aging 
condition. The above observation would lead one to 
suspect the wafer was cracked in the beginning. 
Such a suspicion is not unreasonable, in view of the 
fact that six of the original 500 transistors received 
for this experiment were found to have fractured 
wafers. 

The second transistor, 8-1046, exhibited an abrupt 
increase in emitter leakage current between 14 and 
20 khr of aging, with no prior indication of impending 
deterioration. The resulting voltage-current relation- 
ship had the character of ohmic conduction (~ 10 
kQ), but differing from the six other transistors clas- 
sified as having ohmic leakage in that the change 
appeared as an abrupt steplike increase. Unfor- 
tunately, in removing the case for optical examination, 
the transistor was damaged. 

With the exceptions of transistors 0-1050 and 
8-1046, which had previously had their cases removed, 
the “unsatisfactory” transistors were subjected to 
125 °C mineral-oil gross-leak detection test and a 
helium fine-leak detection test. Only transistor 1- 
1070 was found to leak, and at a rate of approximately 
200 < 10-8 std. cm*/s. (He). 


6. Changes in Leakage Current Related to 
Deterioration 


Table 1 lists the values relating early fractional 
changes in junction leakage current to deterioration 
which led to “unsatisfactory” performance. It is of 
interest to note that the two transistors which ex- 
hibited the largest early fractional changes were the 
only two which later developed avalanche conduction. 
The large fractional change in leakage current oc- 
curred within the first 340 hr of aging. 

The significance of the relationships between frac- 
tional changes and deterioration is more fully revealed 
in figures 1,2, and 3. Figure 1 relates aging behavior 
to the fractional increase in junction leakage current 
associated with the first 340 hr of the aging stress. 
Both junction currents were measured, but only the 
larger of the two. fractional changes was _ plotted. 
The distribution of the ratio values for the entire 160 
transistors is shown by the major histogram. The 
smaller receding histograms show cumulatively the 
ratio values for those transistors which had become 
“unsatisfactory” by the time indicated. The cross- 
hatch shows when a transistor was first observed to be 
“unsatisfactory.” 

Similarly, figure 2 relates the change produced by 
the first 1,000 hr of aging to later behavior. A com- 
parison of the two figures shows that, in general, the 
1,000-hr ratio values of the “unsatisfactory” units 
fall more consistently in the upper range of the ratio 
distribution than do the 340-hr ratio values. In figure 
2, the correspondence between relatively high ratio 
values and deterioration is clearly more pronounced 
for those units which deteriorated more rapidly, as 
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FicurRE 1. Histogram of ratio values showing change in IEBO or 


ICBO after 340 hr of aging. 
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FIGURE 2. Histogram of ratio values showing change in 1EBO or 


ICBO after 1,000 hr of aging. 


shown by the histograms through 8.7 khr. The rela- 
tionship becomes more tenuous for those units which 
last longer. as indicated by the histogram distributions 
spreading down scale with increasing time. 
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FiGURE 3. Histogram of ratio values showing change in IEBO or 


ICBO between 340 and 1,000 hr of aging. 


Figure 3 relates the change observed between 340 
and 1,000 hr of aging to later behavior, thereby ex- 
cluding the more pronounced effects associated with 
the first several hundred hours of aging. The distribu- 
tion of values was so much more compact that the class 
interval of figure 3 was reduced to roughly one-third 
that used in figures ] and 2 in order to obtain a reason- 
able representation of the distribution. 


The relationships in figures 1 through 3 suggest 
means for the identification of transistors likely to 


deteriorate. These data are based on 100 °C stress 
for screening as well as for subsequent aging. It 
should be borne in mind that this stress is higher than 
that which would normally be encountered in appli- 
cations. A convenient and reasonable approximation 
may be made of expected aging performance by assum- 
ing that the rate of deterioration doubles for every 10 
°C increase in junction temperature.2. For example, 
the rate of deterioration at 55 °C may be conservatively 
estimated to be Y10 that at 100 °C. Figures in this 
report, therefore, may serve to roughly approximate 
aging under more normal conditions (say 55 °C) by 
multiplying the time scale (associated with 100 °C 
aging conditions) by 10. Reexamining figure 2 with 
the expanded time scale in mind indicates that the 


? Peck, D.S., Semiconductor Reliability Predictions from Life Distribution Data, ch. 5 
of Semiconductor Reliability, edited by Shwop. J. E. and Sullivan, H Reinhold Pub 
Corp., New York, 1961) 
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first three “unsatisfactory” units might have caused 
trouble in an application by not providing “satis- 
factory” service for at least 1 or 2 years. On the other 
hand, the remainder of the transistors would have 
become “unsatisfactory” only after a period of service 
in excess of 10 years. The fact that the first three 
“unsatisfactory” transistors have unusually large 
ratio values strengthens the expectation that the rela- 
tionship between ratio values and later preformance 
may serve as the basis for a practical screening test 
to detect transistors likely to exhibit junction leakage 
current deterioration. 

Since the suggested screening procedures stem 
from an a posteriori examination of all the data, no 
significant tests of the effectiveness of these pro- 
cedures can be made using these same data. A 
second, smaller experiment described later in this 
paper provides a set of independent data which is 
used for this purpose. However, in order to demon- 
strate the observed relationships of the first experi- 
ment in a screening test context, an example of a 
screening analysis is presented. 


7. A Screening Example 


The ideal screening test serves a dual goal: one is 
to reject all units which would, in fact, fail to meet 
ensuing operational requirements for at least a speci- 
fied period of time, and the other goal is to not reject 
any units which would meet such requirements.  Fail- 
ure to achieve the former goal through the erroneous 
acceptance and use of faulty units would lead to unre- 
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liable performance. Because this is of obvious and 
immediate concern to the consumer, this type of error 
is usually associated with a consumer’s risk. On 
the other hand, failure to achieve the /atter goal 
through the erroneous rejection of satisfactory units 
would directly penalize the producer of such units. 
This type of error is therefore associated with a pro- 
ducer’s risk. It is evident that one or the other of 
these errors can be made as small as is desirable at 
the expense of increasing the other. A graphic view 
of the interplay between these two errors for various 
choices of performance time requirements and screen- 
ing criteria is presented in figures 4, 5, and 6. The 
information was derived from figures 1, 2, and 3, 
respectively. The three groups of histograms in each 
figure illustrate the effects of choices of screening 
limits by application of three different magnitudes of 
fractional change upon which to base the decision to ° 
reject or accept individual transistors. The indicated 
sets of limits are arbitrary and were chosen so as to 
provide a reasonable basis for discussion. 


Each bar in figures 4, 5. and 6 is a complete analysis 
in itself, relating screening results with aging perform- 
ance. The screening test identifies those transistors 
to be rejected. The number of rejected units is indi- 
cated by the height of the bar above the axis. Super- 
imposed on the bar is the aging performance record 
whereby the individual transistors are classified as 
“unsatisfactory” if they fail to meet the end-of-life 
requirements within the time period indicated on the 
abscissae. For example, the second bar from the left 
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FIGURE 4. 


Screening analyses using relative change in values of 1EBO or ICBO after 340 hr of aging. 
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AGING PERIOD IN KILOHOURS 


FIGURE 5. Screening analyses using relative change in values of IEBO or ICBO after 1000 hr of aging. 


TOTAL OF I60 TRANSISTORS 





NSATISFACTORY 








UNSATISFACTORY 





40 4IBWNN 


INCREASE = 5%, REJECT 16 


SHOLSISNVYL 








INCREASE = 10% REJECT 8 

























































































AGING PERIOD IN KILOHOURS 


FIGURE 6. Screening analyses using relative change in values of L\EBO or ICBO between 340 and 1,000 
of aging. 
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in figure 4 indicates that after 340 hr of aging the 
junction leakage current (IEBO or ICBO) increased 
15 percent or more in 3] of the transistors out of the 
total of 160. Two transistors out of the 160 were 


“unsatisfactory” after 2 khr of aging at 100 °C, one of 


which was within the group of 31 rejects and is repre- 
sented by the black portion of the column. The second 
“unsatisfactory” transistor was not among the rejected 
groups and is represented by a dotted column extend- 
ing below the axis. There remains, of course, the 
total number (128) of “satisfactory” units accepted 
correctly. The column in our example with the three 
area distinctions thus indicates: 

1. By the height of the bar over the axis, the num- 
ber (31) of transistors rejected: 

2. by the white area, the number (30) of “‘satis- 
factory” units rejected in error, which is related to 
the “producer’s risk”’; 

3. by the black area, the number (1) of “unsatis- 
factory units rejected correctly: 

4. by the dotted area, the number (1) of “unsatis- 
factory” units accepted in error, which is related to the 
“consumer's risk.” 

It is important to recognize that the time scales of 
figures 4, 5, and 6 refer to performance at 100 °C, a 
temperature which can be expected to produce more 
rapid aging than those ordinarily encountered in normal 
usage. Aging corresponding to a junction tempera- 
ture of 55 °C would be more realistic and may be con- 
servatively estimated by multiplying the time scale 
by a factor of ten. The basis for this procedure was 
discussed in the previous section. 

A comparison of figures 4, 5, and 6 sets figure 4 
apart by the relatively larger number of “unsatis- 
factory” units accepted (associated with “consumer's 
risk’’), even though the corresponding errors associated 
with “‘producer’s risk” are roughly comparable. It 
appears that a substantial reduction in the risk to the 
consumer is gained (without unduly increasing risk 
to the producer) by extending the screening stress to 
1,000 hr rather than terminating at 340 hr. 

There are several distinctions which should be noted 
between figures 5 and 6. The difference in screening 
results between the 5 percent steps of figure 6, despite 
the smaller increments, is more pronounced than that 
between the 15 percent steps of figure 5. This may 
have been expected because 340- to 1,000-hr ratio 
values of figure 6 form a more compact distribution 
(shown in figure 3) than the values for figure 5 (shown 
in figure 2). The effective use of the information in 
figure 6 requires a finer measurement distinction. 
The need for finer distinctions is obviously undesirable: 
as the necessary resolution approaches the magnitude 
of measurement error, the effectiveness of fine dis- 
tinctions evaporates. Stated in other words, one of 
the desirable qualities of a screening parameter is 
that it have a broad distribution, and one of the essen- 
tial qualities is that its dispersion be large relative to 
measurement error. For these reasons the criteria 
represented in figure 5 are favored. An additional 
desirable quality is the apparei:t lack of sensitivity of 
screening results to the choice of the ratio change, 
whether it be 30 or 60 percent. 


On the basis of the factors discussed above, it was 
decided to eliminate the 340-hr measurement period 
and use only the observed relative change after being 
subjected to an accelerated aging stress for 1,000 
hr, as illustrated in figure 5. 


8. Second Experiment 


The value of the screening concept discussed in this 
report depends, to a large extent, upon the relation of 
early-changes versus deterioration being a general 
behavioral characteristic, a characteristic which other 
types of transistors and transistors of different manu- 
facture have incommon. The completion of a modest 
second experiment tends to confirm that the relation- 
ship does exist for other types and manufacture of 
germanium alloy transistors. 

As a sequel to the first experiment, a second was 
initiated, wherein approximately 125 °C junction 
temperature screening stress was applied for 1,000 hr 
and then reduced to 100 °C for continued aging. The 
parameters measured were identical, and the form of 
the data was similar to that of the first experiment. 
It was, therefore, a simple matter to apply the screen- 
ing criteria to the data without foreknowledge of the 
actual performance. For this reason, and the fact 
that the data of the second experiment are independent 
of those used (in the first experiment) in development 
of the screening test, the second experiment provides 
a valid basis for evaluation of the effectiveness of the 
screening procedure. 

The second experiment consisted of applying a 
screening test and then aging all units for nearly one 
year. The experiment included a group of 40 type 
R212 transistors of manufacture B, 40 of manufacture 
C, and 40 type 2N396 transistors of manufacture E. 
The 2N396 units are the same type and manufacture 
as those used in the first experiment. Furthermore, 
an effort was made to compare results of two different 
forms of screening stress. Half the transistors were 
exposed to a 1,000-hr bake screening stress of 125 °C 
with no applied power. The remaining half of the 
transistors were exposed to a contrasting screening 
stress having a low ambient temperature combined 
with a high power dissipation calculated to produce 
the same 125 °C junction temperature. The latter 
stress was obtained with an ambient temperature of 
—65 °C (which is a commonly accepted temperature 
test point) combined with the calculated bias power 
of 760 mW. 

The calculations were based on measured values of 
thermal resistance. The thermal resistance of the 
units of manufacture B was somewhat higher than that 
of the other two groups. It was decided, since ratings 
were similar, to apply the same bias power for all 
three groups, even though the estimated rise in 
junction temperature would be approximately 5 percent 
greater for the units of manufacture B. This later 
proved to be a crucial decision resulting in’ rapid 
deterioration of units of manufacture B. 

After completion of the two screening tests, the two 
groups were again subdivided, half of each group being 
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subjected to 6,000 hr of 100 °C shelf aging and the 
remaining half subjected similarly to power aging. 
The estimated 100 °C junction temperature for power 
aging was obtained with the combination of 70 °C a 
ambient temperature and 120 mW of bias power (con- foe, NR UE BR 
dition 14 in table 3). —e = me mane 
Figure 7 summarizes the results of the 125 °C bake et 
stress screening test. Similarly, figure 8 summarizes 
the results of the power stress screening test. The 
criterion for “satisfactory” performance was that IEBO 
and ICBO at 10 V must be no greater than 5 wA, as 
called for in the procurement specification for R212 
transistors. Referring to the histograms, which show 
the distribution of fractional changes in junction leak- 
age (IEBO or ICBO, whichever is the larger) induced 
-by the screening stress, in figure 7 the pronounced 
mode lies below 1.0 and contrasts sharply with the 
bimodal distribution in figure 8, where both modes 
lie above 1.0. This indicates that the 125 °C shelf 
stress produced a reasonably uniform decrease in 
leakage current for all three manufacturers’ tran- 
sistors, whereas the nominal 125 °C power stress pro- 
duced nonuniform increases in leakage current. The be 
nonuniformity was due mainly to the R212 transistors “7 
of manufacture B exhibiting substantial deterioration: HISTOGRAM OF RATIO VALUES SHOWING CHANGE IN 
the leakage current increased in value to the extent pa AFTER 1,000 HOURS OF SHELF AGING 
that it exceeded 5 pA in 7 out of 20 transistors. Since , 
these seven units no longer met the criterion for “‘satis- 
factory” performance, they were removed from the 
experiment. It is likely that the deterioration of these 
units was caused by excessively high junction tempera- 
ture. As was mentioned earlier, the thermal resist- 
ance of the transistors of manufacture 8 was somewhat 
higher than that of the others and was calculated to 
produce a junction temperature more nearly 135 °C GO TRANSISTORS AGED 7K HOURS 
than the estimated 125 °C of the others. In sharp (3Q.UNITS '00%C_SHELF_AND 30 UNITS 00°C POWER 
contrast to the above described deterioration under ee es 
the power screening stress, the transistors of manu- a | a The | 
facture B exhibited the least degradation during the mee? | 5 | [53°] 
125 °C shelf screening stress or the prolonged 100 °C 
shelf or power aging stress. The collective results 
of this experiment indicate that the effects of the 
power screening stress used was critically sensitive 
to the difference in thermal resistance and could have 
led to erroneous conclusions regarding the relative 
merits of the units of manufacture B. These results 
illustrate the possible hazards of applying uniformly, 
to a given type of transistor, excessive power as a 
means for screening or accelerating aging effects. 
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FIGURE 7. Analysis of 125 °C screening test (shelf). 
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Referring again to figures 7 and 8, the relationship 
between fractional changes in leakage current and 
subsequent aging is shown for aging periods extending 
from 2,000 to 7,000 hr, as was done in a similar fashion 
in figures 1, 2, and 3. Although two aging conditions 
were used, a 100 °C bake stress and a 100 °C junction 
temperature obtained with applied power, the limited 
volume of data indicated no significant difference in 
aging behavior and, therefore, data from the two con- 


“0 13 16 19 227 
ditions were combined. The effectiveness of the HISTOGRAM OF RATIO VALUES SHOWING CHANGE IN 1EBO 


OR ICBO AFTER 1,000 HOURS OF POWER AGING AT ESTI- 


screening process is illustrated in the two boxes in MATED 125° JUNCTION TEMPERATURE 


the upper portions of figures 7 and 8. The screening 
criterion of 1.6 ratio or greater, in the upper box, was FIGURE 8. Analysis of 125 °C screening test (power). 
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chosen because it is the largest of those used in the 
first experiment, and therefore allows comparison of 
the results. Since the screening stress was larger 
(125 °C), a second analysis was made using a larger 
value as a screening criterion. The second value, 
2.2, was the largest one available withcut reprocessing 
the data and, fortunately, appears to be a reasonable 
choice. 

The footnote under the boxes of figure 7 indicates 
that one unit of manufacture C failed, i.e., became 
“unsatisfactory,” during the shelf screening test. 
Although such a failure may offhand appear to indicate 
rapid and destructive deterioration at 125 °C shelf 
stress (as was discussed for the 7 units of manufacture 
B which failed under the power screening stress), 
such a conclusion is surely questionable in light of the 
fact that one out of four of the same manufacturer's 
units also failed while subjected to 25 °C shelf storage, 
which was part of an experiment control. 

Comparison of the screening summaries indicates 
superior results for the shelf-stress screen of figure 7. 
Such comparisons are usually difficult because of the 
several factors which should be considered. A better 
intuitive understanding may be reached for judging 
the effectiveness of the different screening tests by 
examining the summary in table 2. This table com- 
pares the relative number of “unsatisfactory” tran- 
sistors present before screening with the number 
present after screening. An indication of the cost of 
the improvement through screening is given in terms 
of the relative number of “unsatisfactory” units 
rejected correctly compared to the total number 
rejected. The latter ratio indicates the (fractional) 
number of units rejected correctly. 


TABLE 2. Summary of effectiveness of screening tests 


T 
| 
Screening test Relative number of * 


(125 °C) units present 


unsatisfactory” 


Rejected correctly 
Total rejected 


T 


Stress 


| 
} 
| Limit Before screening i After screening 
Shelf | 1.6 39 | Yss or 1.8% Yq or 75% 


| 22 | Vse or 1.8% Ys or 100% 
| ¢ sa 1 -q 1 
> 


Power 37 or 2.7% 


Ys2 or 1.9% 


16 or 6% 
“1 or 100% 





This summary shows that for shelf stress screening, 
the two screening limits (1.6 and 2.2) were almost 
equally effective. For the power stress screening, 
however, the results were very sensitive to the choice 
of limit. The 2.2 value was preferred by far because 
the 1.6 limit caused 15 “satisfactory” units to be re- 
jected in error. This is a relatively poor showing 
compared to the other three outcomes. 

The results of the shelf screening test, even though 
the screening stress (125 °C) is somewhat more severe 
than that used in the first experiment (100 °C), support 
the general findings of the earlier experiment. The 
fractional change in leakage current appears to pro- 
vide a useful measure of expected performance. 


9. Conclusion 


A study of the aging behavior of low-power ger- 
manium alloy switching transistors has revealed a 
relationship between small changes in junction leakage 
current, associated with a brief aging stress, and later 
deterioration in performance. Other parameters 
studied were ICES, ICER, ICEO, and 1+ hFB, which 
were measured, and the voltage sensitive component 
of the leakage current, the intrinsic leakage current, 
and hFE*, which were computed. With the possible 
exception of hFE*, none of these parameters appeared 
promising for stress-screening purposes. 

These findings apply to transistors which had been 
accepted for stringent military requirements and were 
considered “high reliability” units. They had con- 
formed to military acceptance specifications which 
had included a series of measurements on each device. 
It was not surprising, therefore, that the NBS 0-hr 
measurements proved to be of little value as a screen- 
ing criterion. Even the variable or noisy character 
of some of the 0-hr measurements offered little hope 
as an indicator of expected performance. Shifts or 
jumps in parameter values while lightly tapping the 
transistor case, however, proved to be quite useful. 
Six type 2N396 transistors having fractured wafers 
or other mechanical faults were clearly identified 
while making preliminary measurements using this 
procedure, and were not included in the aging study. 

The observed junction leakage current relationship 
to aging appears characteristic of several types and 
manufactures of germanium-alloy transistors. This 
relationship may provide the basis for a promising 
nondestructive screening procedure which would 
serve for the identification of germanium-alloy tran- 
sistors likely to deteriorate through excessive growth 
in junction leakage current. Since leakage current 
deterioration has continued to be one of the more 
common forms of degradation in transistor perform- 
ance, a screening procedure which is effective in 
guarding against such a fault would be of substantial 
value. 

The proposed screening procedure involves the de- 
termination of relatively small changes in junction 
leakage current, increases of the order of 15 percent 
or more, associated with 1,000 hr of aging at a shelf 
(bake) stress of 100 °C. Because the leakage current 
changes of interest are small, high demands are placed 
upon measurement repeatability. A number of fac- 
tors are important in maintaining the necessary long- 
term repeatability. There is need for operating pro- 
cedures which provide a record of proof of stable 
performance through a history of frequent calibra- 
tion and through a history of measurements on a group 
of “control” transistors. Also, the measurement 
procedure should treat each transistor alike. A 
stabilization or “soaking” period of a fixed duration 
(24 hr+2 hr) under measurement ambient conditions 
prior to measurement should be incorporated in the 
test procedure so as to minimize the transient recovery 
effects associated with the removal of the transistor 
from a high-temperature aging stress environment. 
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Relatively close ambient temperature control during 
measurement is essential. The very close control 
of +0.01 °C of the NBS measurement facility, although 
proven to be a key feature in the aging studies, is 
unnecessarily stringent. Commercial equipment hav- 
ing control with +0.1 or 0.2 °C is adequate for the 
leakage-current bake-stress screening described in 
this paper. 

The initial choice of value of fractional change in 
junction leakage current as the screening limit is 
necessarily somewhat arbitrary. There is evidence 
that even transistors bearing identical type numbers 
but of different manufacture might require different 
screening limits. Once an initial choice is made, 
however, an a posteriori aging test of all rejected units 
and sample groups of accepted units is recommended 
as an evaluation procedure for monitoring the screen- 
ing effectiveness and for improving the choice of 
screening limits. These procedures are now being 
followed at NBS in a further study of the proposed 
screening procedure as applied to much larger groups 
of transistors and to germanium transistors of dif- 
ferent construction. 


10. Appendix 


A total of 516 germanium low-power alloy-junction 
switching-type transistors (type 2N396), obtained 
from a single manufacturer, were aged under 13 dif- 
ferent combinations of ambient temperature, collector 
voltage, and collector current. The thirteen aging 
conditions are summarized in table 3, from which it 
may be seen that these conditions form an incomplete 
factorial experiment. Each condition is identified by 
a number in the lower left corner of each box. Each 
condition contained 40 transistors mounted on four 
plugs. 

Immediately prior to being subjected to its aging 
environment, each transistor was measured twice, a 
morning group of 10 parameter measurements being 
replicated in the afternoon of the same day. During 
the course of the experiment, the transistors were 
periodically removed from their aging environments 
and the 10 parameters were measured four times 
before the aging treatments were again applied. 
The measurement period began with removal of the 
transistors from the aging environment, after which 
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they were permitted to stabilize for 24 hr at 25 °C. 
During the morning of the second day of the measure- 
ment period, the 10 parameters were measured and a 
replicate set of measurements made during the after- 
noon. The transistors were then permitted to stabilize 
for an additional seven days, the a.m. and p.m. sets of 
measurements being repeated during the eighth day 
of the measurement period. A computer was used 
to compare the a.m. and p.m. replicate measurements 
and print out all which changed by more than 2 percent. 
This allowed all such measurements to be reviewed 
and, if necessary, repeated before returning the 
transistors to the aging environment. 

The first three measurement periods, 340, 1,000, and 
2,000 hr, were selected to correspond to those com- 
monly used in aging schedules and specifications. 
The remaining periods were selected to provide ap- 
proximately equal log-time intervals between the 
measurement periods. The intervals so determined 


are a fixed fraction of the elapsed time, and for this 


experiment, were roughly 40 percent of the elapsed 
time. This procedure gave the following total elapsed 
times at which the measurement periods fell: 3.4, 
4.7, 6.0, 8.7, 11, 14, and 20 khr. 

The d-c parameters measured, using automatic 
equipment, were as follows: IEBO and ICBO at 2 and 
10 V; ICES, ICER (R=10 kQ), and ICEO at 10 V; 
normal and inverse 1+hFB at 10 mA and 5 V. The 
tenth parameter, 1+hfb, was measured using a 270- 
Hz small signal. From these measurements a num- 
ber of additional parameters were computed. The 
computed parameters provided estimates of: the 
voltage-sensitive components of the leakage currents, 
the intrinsic leakage currents, and (1+hFB)* (the 
asterisk indicating that the leakage current was sub- 
tracted from the measurement of 1 + hFB) and, finally, 


hFE*. 
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Effect of hydrostatic pressure on the refractive indices of 
some solids, R. M. Waxler and C. E. Weir. J. Res. NBS 694A (Phys. 
and Chem.), No. 4, 325-333 Yuly—4ug. 1965), 70 cents. 
Measurements were made on refractive index changes with hydro- 
static pressures between | bar and | kbar using the helium yellow 
line. The materials studied were: KBr, NaCl, LiF, diamond, MgO, 
quartz, AloOg. and three silicate glasses. All the materials increased 
in refractive index with pressure except diamond, MgO, and Al.Os 
which decreased, and LiF which did not change. The results were 
compared with photoelastic measurements, and Pockel’s geometric 
theory of photoelasticity was substantiated as well as Mueller’s 
physical theory. The data show that the ratio of change of polariz- 
ability with density is greater for solids having stronger interatomic 
repulsive forces. Volume and temperature coefficients of polar- 
izability were evaluated for the cubic crystals and glasses. The 
thermo-optic behavior of crystals and glasses is discussed. 


An absolute temperature scale from 4 °K to 20 °K deter- 
mined from measurements with an acoustical thermometer, 
H. H. Plumb and G. Cataland, J. Res, NBS 69A (Phys. and Chem.), 
Vo. 4, 375-377 (July—Aug. 1965), 70 cents. 

At the National Bureau of Standards an acoustical thermometer has 
been used to obtain values of temperature at every degree from 5 to 
20 °K as a basis for a temperature scale. This scale has been com- 
pared with four other temperature scales in the region from 10 to 
20 °K. Since the acoustical thermometer is an entirely new ap- 
proach to precision thermometry in this range the comparison of 
its resulting scale with scales based upon gas thermometry from 
12 to 20 °K is significant. Indications of inaccuracies in equi- 
librium hydrogen vapor pressure scales, and also in the He* vapor 
pressure scale, are presented. 


Current research on preservation of archival records on 
silver-gelatin type microfilm in roll form, C. S. McCamy and 
C. 1. Pope, J. Res NBS 694A (Phys. and Chem.), No. 5, 385-396 (Sept.— 
Oct. 1965). 70 cents. 

This is a progress report on a study of aging blemishes on micro- 
film and measures to prevent their formation. Recommendations, 
based on current information, include: avoiding excessive densities, 
avoiding physical damage to the film, the use of potassium iodide in 
the fixer, thorough washing. uniform drying, storage at low tem- 
perature and humidity in sealed containers, avoidance of air-borne 
reactants, increased archival use of positive copies, and periodic 
inspection. 

A wide-range (up to 10'°P) rotating cylinder viscometer, 
A. Napolitano, P. B. Macedo, E. G. Hawkins, J. Res. NBS 69A 
(Phys. and Chem.), No. 5, 449-455 (Sept.—Oct. 1965), 70 cents. 

The range of a high temperature rotating cylinder viscometer has 
been extended so that viscosity measurements can be made between 
10° to 10"° poises to within 2 percent accuracy. This involves three 
different techniques. After a calibration constant has been deter- 
mined for the apparatus in the standard oil range, the other two 
constants for the higher viscosities were self-determined. In making 
measurements at the low viscosities, 10° to 10° poises, the outer 
cylinder is rotated at constant speed and the torque measured on 
the inner cylinder; from 1045 to 107° poises the inner cylinder is 
rotated through an angle and timed as it returns to its zero position, 
and from 1055 to 10" poises the inner cylinder is driven through an 
angle at constant torque and timed as it traverses this angle. The 
Fulcher equation viscosity values obtained from previous measure- 
ments on Standard Glasses No. 710 and No. 711 by the fiber elonga- 
tion and rotating cylinder methods have been compared with the 
present results. 


Internal friction in ZrO: containing CaO, J. B. Wachtman, Jr., 
and W. C. Corwin, J. Res. NBS 69A (Phys. and Chem.), No. 5, 457- 
460 (Sept-Oct. 1965), 70 cents. 

Internal friction was measured on polycrystalline ZrO. containing 
4, 7, 10, 13, 16, and 20 mole percent CaO. In the range 10 to 20 
percent a symmetrical internal friction peak occurs with its maxi- 
mum at about 300 °C at 1 kHz; below 10 percent a nonsymmetrical 
peak occurs at somewhat higher temperature which may be asso- 
ciated with the presence of two phases. The symmetrical peak 
in as-sintered specimens has the same dependence on CaO content 
as that reported for electrical resistivity at 1000 °C and is probably 
caused by motion of oxygen vacancies. Annealing at 1000 °C for 
1000 hours causes reduction of the peak height but with a different 
composition dependence from that reported for annealing of elec- 
trical resistivity. The results suggest that oxygen vacancies in 
calcia stabilized zirconia exist in several states of binding but that 
detailed models are likely to be complex. 


On Kirchhoff’s law and its generalized application to absorp- 
tion and emission by cavities, F. J. Kelly. J. Res. VBS 69B (Math. 
and Math. Phys.), No. 3, 165-171 (July—Sept. 1965), 75 cents. 
Several authors have made the assumption that Kirchhoffs Law 
holds for the apparent local spectral emittance and apparent local 
spectral absorptance of any point on the interior surface of a cavity. 
The correctness of this assumption is demonstrated under certain 
general conditions, and its practical application to the calculation 
of the total flux absorbed by a cavity or spacecraft is discussed. 
A further application to the case of a nonisothermal cavity or space- 
craft is derived. By this derivation an easy method for determining 
the total flux emitted from such a nonisothermal cavity is found 
when the distribution of the apparent local spectral emittance of 
the isothermal cavity is known. The economy and versatility of 
this method is shown by the calculation of the total flux emitted 
from a nonisothermal cylindrical cavity for several arbitrary cases 
of temperature distribution on the interior surface of the cavity. 
Finally, the integral equation for a diffuse cavity whose wall emit- 
tance varies with position on the wall is transformed to an equation 
having a symmetric kernel. 


Ground-based passive probling using the microwave spec- 
trum of oxygen, E. R. Westwater, Radio Science J. Res. NBS/ 
USNC-URSI 69D, No. 9, 1201-1211 (Sept. 1965), $1.00. 

The determination of the kinetic temperature structure of the tropo- 
sphere from ground-based measurements of oxygen emission spectra 
in the microwave region is discussed. Molecular absorption prop- 
erties of oxygen and water vapor are reviewed. A new “inversion” 
technique is described. This technique uses a least squares- 
iteration solution which is applicable to Fredholm integral equations 
of the first kind. The inversion technique is used to reconstruct 
two types of tropospheric temperature profiles. The effects of cer- 
tain types of errors in the brightness temperature on the derived 
temperature distribution are computed. 


Response of NBS microwave refractometer cavities to at- 
mospheric variations, R. O. Gilmer, R. E. McGavin, and B. R. 
Bean, Radio Science J. Res. NBS/USNC-URSI 69D, No. 9, 1213- 
1217 (Sept. 1965), $1.00. 

The sampling cavity of the NBS refractometers has been evaluated in 
a wind tunnel, a waterflow tunnel, and in the free atmosphere. The 
cavity is slightly velocity-sensitive. It is also aspect-sensitive, an 
error being introduced when the wind is oblique to the axis of the 
cavity affecting the level of the spectrum. The characteristic 
flushing length of the cavity in normal operation is estimated as 
0.75 meters when the wind is into the cavity and considerably more 
when the wind is at an oblique angle. A modified cavity has been 
developed which appears to substantially improve the response of 
the refractometer. 
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Effects of rocket outgassing on rf experiments, W. Pfister 
and J. C. Ulwick. Radio Science J. Res. NBS/USNC-URSI 69D, No. 
9, 1219-1225 (Sept. 1965), $1.00. 

Peculiarities in the results of electron density measurements with 
rockets were notieed on several occasions. They were found to 
be effects of rocket outgassing. When the fuel valves of Aerobee 
rockets are not shut after burnout, residual fuel leaks out and dif- 
fuses into the atmosphere. Electrons attach to the fuel and de- 
crease in concentration while the positive ions are not affected. 
Outgassing effects were observed on probe experiments flown in 
rockets as well as on propagation experiments between rocket 
and ground. Quantitative estimates indicate that the leaking fuel 
is a very effective attaching agent. Some solid fuel rockets also 
show similar effects. 


Measurement of group velocity of 17.8-kHz VLF radio 
waves, W. D. Westfall. Radio Science J. Res. NBS/USNC-URSI 
69D, No. 9, 1235-1238 (Sept. 1965), $1.00. 

Group velocity of 17.8 kce/s VLF radio waves is determined by 
measuring the difference in time of arrival of the transmitted signals 
which travel over complementary great circle paths. The velocity 
is found to be 3.04 X 108 meters/second with an estimated error of 
approximately + 0.10 X 10° meters/second. Differences in measure- 
ment .echnique and general environmental conditions between this 
experiment and a previous measurement of 18 kc/s group velecity 
which was made in 1949 is discussed. 


Propagation of pulses in dispersive media, J. R. Wait, Radio 
Science J. Res. NBS/USNC-URSI 69D, No. 11, 1387-1401 (Nov. 
1965), $1.00. 

The paper deals with various approximate procedures for calculating 
the distortion of a pulse after it has propagated through a dispersive 
channel such as a waveguide. The methods used for evaluating the 
integrals utilize a stationary phase principle. Both impulsive-type 
and quasi-monochromatic sources are considered. It is shown that, 
in most cases, the transient response may be obtained straightfor- 
wardly from the shape of the phase versus frequency characteristics 
of the system. Some attention is given to the complication which 
arises when the group velocity has an extremum as a function of 
frequency 


Interference rejection capability of a switched radiometer, 
R. E. Clapp, Radio Science J. Res. NBS/USNC-URSI 69D, No. 11, 
1425-1429 (Nov. 1965), $1.00. 

The increasing congestion of the radio-frequency spectrum makes 
necessary a re-examination of current techniques used in radiometric 
receivers for radio astronomy. Square-wave switching, as currently 
used, modulates incoming CW signals and introduces sidebands 
which can intrude on the presently available windows in the spec- 
trum. This effect can be treated as an apparent deterioration of 
an ideal square-sided filter. Examples are given, for different filter 
bandwidths and chopping frequencies. Where radiometry is inter- 
ference-limited, sinusoidal chopping or a correlation technique 
avoiding chopping is strongly recommended. 


Atmospheric breakdown limitations to optical maser propa- 
gation, R. G. Tomlinson, Radio Science, J. Res. NBS/USNC-URSI 
69D, No. 11, 1431-1433 (Nov. 1965), $1.00. 

It is shown that air breakdown by focused optical maser pulses is 

dependent on gas pressure. It would appear that megawatt type 

og of 10 to 100 nanosecond duration could produce air break- 
own. 


Electrical engineering units and constants, VBS Misc. P1/ 
268 (June 1965). 10 cents a copy; $6.26 per 100 copies. 

This publication is a compilation of electrical engineering symbols 
and units and physical constants used in electrical engineering. 


Examination of liquefied petroleum gas liquid-measuring 
devices, M. W. Jensen, NBS Handb. 99 (Apr. 23, 1965), 35 cents. 

A manual for State and local weights and measures officials, de- 
scribing the devices, testing equipment and its calibration, inspecting 
and testing procedures, and a reporting system. 


Steel medicine cabinets, VBS CS267-65 (Mar. 1, 1965), 10 cents. 
This Commercial Standard covers the requirements for materials, 
construction, finishes, accessories, minimum sizes, methods of test 


and the labeling of sliding-door and swinging-door steel medicine 
cabinets. Cabinets having a combination of mirror and enclosed 
sliding door shelves, and known variously as mirror vanity, vanity 
mirror and vanity shelves, or otherwise designated, are not covered 
by this standard. 


Hide trim pattern for domestic cattlehides, VBS CS268-65 
(Apr. 4, 1965), 10 cents. 

This Commercial Standard covers the trimming of hides of the 
domestic bovine species generally described as “domestic packer 
cattlehide”, and other comparable hides. 


Phase and amplitude diversity in over-water transmissions 
at two microwave frequencies, H. B. Janes, A. W. Kirkpatrick, 
D. M. Waters, and D. Smith. NBS Tech. Note 307 (Apr. 12, 1965), 
10 cents. 

Propagation measurements at 9.4 and 9.2 Ge/s were made over a 
17 km, line-of-sight, over-water path at Eleuthera Island in the 
Bahamas. The purpose was to study the signal amplitude and phase 
variations at the two frequencies, and in particular, the variations 
in the phase difference and amplitude ratio (“differential amplitude”) 
of the two frequencies. The results include power speetra of phase 
and phase difference variations (from about 10~* to 4 c/s) and cumu- 
lative distributions of amplitude, differential amplitude, and phase 
difference. 


A bibliography of experimental saturation properties of the 
cryogenic fluids, N. A. Olien and L. A. Hall, NBS Tech. Note 309 
(Apr. 9, 1965), 60 cents. 

A bibliography of 507 references to experimental work is presented 
for the properties of the cryogenic fluids in the solid, liquid, and 
vapor phases at saturation. The cryogenic fluids included are 
helium, hydrogen. neon, nitrogen, oxygen, air, carbon monoxide, 
fluorine, argon, methane, and isotopes of helium, hydrogen, and 
methane. Each article has been reviewed and coded for proper- 
ties, method of presentation of the data, and temperature range. 
An index lists each fluid in five categories: solid-solid transition, 
solid-liquid transition, solid-vapor transition, liquid-vapor  transi- 
tion, and triple point. For helium the liquid-liquid transition is 
also included. Each category is indexed by the properties: pres- 
sure-temperature data, density and heat capacity of all phases at 
saturation, and latent heats. 


Watermarks for papermaking by electroforming, J. P. 
Young, Tappi 48, No. 1, 364—37A (Jan. 1965); Metal Finishing 63, 
No. 4, 79-82 (Apr. 1965). 

A process has been developed whereby a watermark is electro- 
formed on a papermaking screen. The pattern of the watermark is 
first outlined on the screen by use of a photosensitive resist or by 
the use of a suitable resist and the silk screen process. The water? 
mark is then formed directly on the screen with electrodeposited 
copper. This process allows transfer to the papermaking screen 
of any line or silhouette pattern, no matter how complicated, from 
a film transparency with a saving of time and skilled labor. 


Latensification in radiographic emulsions, W. L. McLaughlin, 
J. Photo, Sci. 13, 1-11 (1965). 
Latensification by means of post-exposures to selected wavelengths 
of light or near infrared is shown to provide as much as threefold 
enhancement of response of some commercial x-ray emulsions to 
ionizing photons. The sensitivity can be similarly enhanced by 
special development procedures. A combination of these tech- 
niques permits a five- or six-fold increase in sensitivity in the low- 
ensity region, for some emulsions. Spectral sensitivity changes in 
the x- and y-ray response are small; no rate dependence is observed; 
variations in total density, however, are found with differences in 
delay period between exposure and post-exposure and between post- 
exposure and processing. The effects may be explained by the 
change in the size and location of development centers resulting 
in an improvement of development activity, and therefore in a 
larger number of developed grains. In some instances, there may 
be an increase in developed grain size. j 


On the information in a microphotograph, C. S. McCamy, 
Appl. Opt. 4, No. 4, 405-411 (Apr. 1965). 
The information concentration of a legible microphotograph made 
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at a reduction ratio of 1200 is computed to be over 2 X 106 bits/mm?. 
Reasons are adduced for using a mosaic of squares having sides 
equal to the reciprocal of the resolving power, each being capable 
of assuming one of two density ranges, as a model of a high contrast 
material for the computation of information capacity. It follows 
that the information capacity is equal to the square of the resolving 
power. This relation gives values of information capacity in good 
accord with the values computed by the Jones theory. The infor- 
mation capacity of a Kodak High Resolution Plate is computed 
to be 3.25 10® bits/mm?. That the product of reduction ratio 
times quality index equals the product of resolving power times the 
height of the lower case e in type to be copied is verified at the high 
reduction ratio, and this relationship is considered useful for engi- 
neering computations involving the copying of commonly used styles 
of type. The use of literal charts designed on the basis of informa- 
tion theory is suggested for the evaluation of microcopying systems. 


Results of National Bureau of Standards corrosion investi- 
gation in disturbed and undisturbed soils, M. Romanoff, 
Proc. Appalachian Underground Short Course, Tech. Bull. 72, 
566-589 (West Virginia Univ., Morgantown, West Va., 1965). 
This paper gives (1) a summary of the extensive investigations which 
have been in progress at the National Bureau of Standards since 
1922 on the corrosion of metals in a wide variety of soils for exposures 
up to 17 years, and (2) a summary of the results of inspections made 
on steel pilings after exposures up to 40 years in various underground 
structures. 

The corrosion data show that there is an appreciable difference in 
the corrosion observed on steel pipe specimens buried at the test 
sites under disturbed soil conditions, and with that observed on 
steel pilings driven in undisturbed soils. 


Technique for determining spin-lattice relaxation times, 
G. A. Candela and R. E. Mundy, Rev. Sci. Instr. 36, No. 3, 338- 
342 (Mar. 1965). 

Experiments were performed so that the actual change in the static 
direct current susceptibility of a paramagnetic material cbuld be 
measured as a function of microwave power absorbed at electron 
spin resonance. The spin-lattice relaxation 7; can be expressed by 


— Xo — Xs) a 
1" Pm M Bae 


where P, is the microwave power necessary to change the field 
independent static susceptibility from Xo to xs. Xs is the steady state 
static susceptibility of the sample when it is absorbing P, amount 
of microwave power. Bac is the homogeneous dc field produced 
by the laboratory magnet, M, is the weight of the sample, and 
I'm is the magnetic constant. The relaxation time is determined 
by directly measuring Pa, Ms B, xo, and xs/xo. The apparatus, 
the experimental techniques, and the accuracy of the measurements 
are déscribed in some detail. 


Two cases of stress corrosion cracking in copper tubing, 
H. L. Logan and G. M. Uginsky, Materials Prot. 4, No. 3, 79-80 
(May 1965). 

A failure in annealed copper tubing used to transport natural gas 
from the main to an individual dwelling was investigated and found 
to be due to stress-corrosion cracking penetrating through the 
tubing wall from the exterior surface. A second failure also diag- 
nosed as stress-corrosion cracking is also reported. In each case 
the area of failure had been encased in a section of steel pipe. 


Nomographs for computing real, imaginary, and absolute 
values of vector ratios, H. S. Bowman, Acoust. Soc. Am. 37, 
No. 4, 751 (Apr. 1965). 

Two nomograms were designed and constructed to compute real, 
imaginary, and absolute values of two-dimensional vector ratios. 


A study of errors in the measurement of microscopic 
spheres, C. P. Saylor, Appl. Opt. 4, No. 4, 477-486 (Apr. 1965). 
Uncertainties in the microscopical measurement of small spheres 
as influenced by index of refraction of the mount, focus, aperture, 
and resolution have been variously reported and denied. The 
present study was conducted with large spheres in order to disclose 
the basic phenomena, but with such optical conditions as would 
ordinarily be used for much smaller objects. The reality of each 


reported error is verified. The principal uncertainties originate 
in the geometrical optics of thick bodies, but these effects tend to 
be confused by limitations in resolving power. 


Definition of “‘ampere” and “magnetic constant”, C. H. 
Page. Proc. IEEE 53, No. 1, 100-101 (Jan. 1965). 

Confusion about the definition of the MKSA system of electrical 
units is related to two meanings of the word “ampere.” This word 
refers to either the physical measurement unit of current, or the 
abstract unit in the mathematical representation of current. 

The magnetic constant is defined as the proportionality constant 
in an equation representing the interaction force between current 
elements: assignment of a numerical value to this constant implic- 
itly defines the physical ampere. 


Electromagnetic measurement, H. W. Lance, SPACE/AERO- 
VAUTICS Res. Develop. Handb. 42, No. 4, 199-202 (1964-1965). 
Important recent developments in electromagnetic measurements 
are helping to fill a serious “measurements gap.” Cesium beam 
frequency standards have been improved, and recently inaugurated 
LF and VLF standard frequency broadcasts permit improvements 
to be realized at receiving stations. In the low frequency region 
improved ratio devices (especially inductive voltage dividers) 
permit laboratories to perform a wide variety of measurements 
based upon relatively few reference standards. In the high fre- 
quency region improved coaxial connectors permit significant 
improvements in the accuracy of most types of measurements in 
coaxial systems. In the microwave region an accelerated radar 
standards program in the NBS Radio Standards Laboratory offers 
the probability of relief in this critical measurement area. The 
increasing maturity of the millimeter wave field and the extension of 
coherent radio measurement techniques to optical frequencies 
offer new challenges in electromagnetic measurements. 


A test of analytical expressions for the thermal emissivity 
of shallow cylindrical cavities, F. J. Kelly and D. G. Moore, 


Appl. Opt. 4, No. 1, 31-40 (Jan. 1965). 


Analytical expressions for the directional thermal emissivity of a 
shallow cylindrical cavity were tested by means of room tempera- 
ture reflectance measurements. The measurements were made 
by placing a paper-lined brass cavity of adjustable depth over the 
specimen opening of a recording spectrophotometer. Reflectance 
curves were obtained from 0.4 uw to 0.75 w with the plunger posi- 
tioned to give cavity depth-to-radius ratios of 0, 0.5, 1.0, 1.5, and 
2.0. The resulting reflectances were then converted to emissivities 
through use of the relation that, for an opaque material, the emissivity 
is equal to one minus the reflectance. Two theoretical expressions 
predicted the emissivity of a cavity with a diffusely reflecting wall to 
within 0.01. A closed-form expression of Gouffé, which forms the 
basis of a shallow-hole method for measuring the emissivity of 
nonmetals at high temperatures, was in as good agreement with the 
experimental measurements as a more rigorous expression derived 
by a different approach. Sizable departures of the cavity walls 
from good diffusers had a negligible effect on the emissivity of the 
cavity at a depth-to-radius ratio of 0.5, but the effect of specular 
components became progressively greater as the cavity depth was 
increased. 


A universal color language, K. L. Kelly, Color Eng. 3, No. 2, 
16-21 (Mar.—Apr. 1965). 

A method of color designation which provides a universal means of 
color communication for people of varying technical and occupa- 
tional backgrounds, using a newly available, basically different 
type of color chart. 


High silica glass, quartz, and vitreous silica, J. S. Laufer, 
J. Opt. Soc. Am. 55, No. 4, 458-460 (Apr. 1965). 

A discussion of confusion in the literature caused by twenty-eight 
terms used to describe high silica glass, natural quartz, cultured 
quartz, and the two classes of vitreous silica. Five clearly defined 
terms are recommended to eliminate error or confusion whenever 
these materials are mentioned. Appropriateness of the chosen 
terms is discussed on the basis of structure of the materials. 


Interference patterns in reverberant sound fields. II, R. V. 
Waterhouse and R. K. Cook, J. Acoust. Soc. Am. 37, No. 3, 424- 
428 (Mar. 1965). 
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Contour charts are given of the sound fields near one, two, or three 
reflecting surfaces at right angles, for random sound incidence. 
The charts show contours of the mean-squared pressure in decibels 
as a function of positien in the region within two wavelengths of the 
boundary surfaces. The boundary conditions at the reflecting 
planes are pressure-reflecting (normal velocity component equal to 
zero) or pressure-release (pressure equal to zero). By reciprocity, 
the charts also give the relative radiation resistance of a point source 
as a function of position near such reflectors. The formulas used 
are tabulated; they extend and in some places correct formulas 
given in an earlier paper [J. Acoust. Soc. Am. 27, 247-258 (1955)]. 
The features of the sound fields near the reflectors are discussed, 
and some applications to acoustical measurements in reverberant 
sound fields are considered. 


The use of terahertz photobeats for precise velocity-of-light 
measurement, Z. Bay and H. S. Boyne, Proc. Enrico Fermi Inter- 
national School of Physics: Course No. XXXI Quantum Electronics 
and Coherent Light 1963, pp. 352—371 (1964). 

A method for the determination of the velocity of light using high 
frequency beats of the order 10'® cps between laser lines is dis- 
cussed. The feasibility of accurately measuring the beat frequency 
by. means of a special phototube similar in design to a fast cathode 
ray tube is examined in detail. It is shown that one might expect 
to’ measure beat frequencies as high as 1.6 X 10" cps if there is no 
appreciable electron scattering in the photocathode. A determina- 
tion_of the wavelength separation AA corresponding to the beat 
frequency w can be made with a scanning Fabry-Perot interfer- 
ometer of a few meters length in order to obtain a precision of one 
in 10° for the velocity of light. 


Potential use of passive probing of atmospheric structure 
by thermal emissions at radio frequencies, B. R. Bean, E. R. 
Westwater, and R. L. Abbott, (Proc. 1963 Intern. Symp. Humidity 
and Moisture, Washington, D.C.), Book, Humidity and Moisture. 
Measurement and Control in Science and Industry, Vol. 2, Appli- 
cations, Sec. VII, Paper 71, pp. 595-608 (Reinhold Publ. Corp., 
New York, N.Y., 1965). 

Microwave resonance lines of water vapor (v=22.235 Ge/s) and 
oxygen (50 Ge/s < v < 65 Ge/s) give rise to strong thermal emissions 
by the atmosphere in the microwave region. The intensity of these 
emissions depends on the water vapor and temperature distributions 
within the atmosphere. This paper explores the possibility of using 
ground-based radiation measurements of these thermal emissions 
to determine the gross water vapor and temperature profiles of the 
lower atmosphere. 


Radio refractometry and its potential for humidity studies, 
R. E. MeGavin and M.J. Vetter, (Proc. 1963 Intern. Symp. Humidity 
and Moisture, Washington, D.C.), Book, Humidity and Moisture. 
Measurement and Control in Science and Industry, Vol. 2, Applica- 
tions, Sec. VI1, Paper 66, pp. 553-560 (Reinhold Publ. Corp., Neu 
York, N.Y., 1965). 

Radio refractometers have been in use for the past decade. Their 
primary use has been in the investigation of the refractive index 
structure as applied to radio wave propagation. There has been 
considerable improvement in these devices in recent years and the 
present status of radio refractometery indicates broader application 
of these techniques. Since the radio refractive index is a function of 
pressure, temperature, and humidity, these devices can be converted 
to hygrometers when suitable compensation for temperature and 
pressure is applied. This technique has been used successfully 
by a number of investigators. The results of some of these investiga- 
tions are reviewed. 


The use of the radio refractometer to measure water vapor 
turbulence, B. R. Bean and R. E. McGavin, (Proc. 1963 Intern. 
Symp. Humidity and Moisture, Washington, D.C.), Book, Humidity 
and Moisture. Measurement and Control in Science and Industry, 
Vol. 2, Applications, Sec. VII, Paper 67, pp. 561-568 (Reinhold 
Publ. Corp., New York, N.Y., 1965). 

The turbulent transport of the radio refractive index is studied theo- 
retically. This study indicates that, under a wide range of condi- 
tions, the turbulent fluxes of the radio refractive index and of the 
absolute humidity are linearly related. 

The initial results of an experiment specifically designed to check 
this conclusion are in agreement with the theoretical analysis. 


Emission spectroscopy, B. F. Scribner and M. Margoshes, 
Treatise on Analytical Chemistry, Ed. I. M. Kolthoff and P. J. 
Elving, 6, Pt. 1, Chapt. 64, pp 3347-3461 (UInter-Science Publ. Inc., 
New York, N.Y., 1965). 

The literature on optical emission spectrometry, including atomic 
absorption and flame photometry, is reviewed for the period since 
the last reviews which appeared in Analytical Chemistry in April, 
1962. Attention is given to cgntributions of basic importance or to 
analytical applications that are new, unique, or extensive. The 
review is divided for convenience into sections including Books 
and Reviews, Spectral Descriptions and Classifications, Spectro- 
scopic Instrumentation, Excitation, Standards and Calibration, 
Isotopic Analysis, Trace Analysis, Analysis of Gases, Analysis of 
Metals, Analysis of Non-Metals, and Literature Cited. The ref- 
erence list totals 565. 


Gas chromatographic determination of the moisture content 
of grain, E. L. Weise, R. W. Burke, and J. K. Taylor, (Proc. 1963 
Intern. Symp. Humidity and Moisture, Washington, D.C.), Book, 
Humidity and Moisture, Ed. A. Wexler, Vol. IV, Principles and 
Methods of Measuring Moisture in Liquids and Solids, Chapt. 1, 
pp. 3-6 (Reinhold Publ. Corp., New York, N.Y., 1965). 

Water is extracted from ground grain by anhydrous methanol as 
in the Karl Fischer titration method of Hari and Neustadt. Deter- 
mination of water in the methanol extract by gas chromatography 
yields a value for the moisture in the grain. Various gas chromato- 
graphic columns and operating conditions were explored to find 
combinations suitable for the adequate separation and accurate 
determination of small amounts of water in large excesses of meth- 
anol. Comparison is made with the Karl Fischer method. This 
procedure should be applicable to other materials if (1) methanol 
will extract all of the water classed as “moisture” and (2) methanol 
does not extract substances that would coincide with either the 
methanol or water peaks on the chromatogram. 


Oxidation of asphalt flux with oxides in nitrogen, P. G. 
Campbell and J. R. Wright. /nd. Engr. Chem. Prod. Res. Develop. 
4, No. 2, 121-128 (June 1965). 

An asphalt flux of Southeastern U.S.A. origin was treated with 
various oxides of nitrogen as the oxidizing gases. Changes in physi- 
cal and chemical properties were measured by infrared spectro- 
scopy, softening points, and microviscometry. Reaction tempera- 
tures did not exceed 166 °C. and gas flow rates were well below 
those used in the commercial air-blowing of asphalts. Nitrous 
oxide and nitric oxide, alone, gave very little reaction. Nitrogen 
dioxide, photosensitized nitrogen dioxide, and a nitric oxide-oxygen 
mixture reacted rapidly with the asphalt flux to produce coating- 
grade asphalts. The reaction time needed to produce a given 
asphalt softening point with photosensitized nitrogen dioxide, for 
example, was 1/20 that needed for air alone. Coating-grade as- 
phalts made with oxides of nitrogen were more stable to photo- 
oxidation, and to film rupture at any given rate of shear, than was 
the commercial asphalt from the same flux. A correlation was ob- 
served between the extent of photo-oxidation of exposed asphalts 
and the changes in viscoelastic properties. 


Radiation-induced polymerization and other reactions of 
n-perfluoropentadiene-1, 4 at high temperature and pres- 
sure, D. W. Brown, J. E. Fearn, and R. E. Lowry, J. Polymer Sci. 
3, Pt. 2, 1641-1660 (1965). 

The radiation-induced polymerization of n-perfluoropentadiene-1, 4 
was studied between temperatures of 100 and 170 °C and pressures 
of 8,000 and 15,000 atm. Kinetic evidence indicates that poly- 
merization occurs by a free radical reaction; the activation energy 
is between 14 and 17 kcal/mole and the activation entropy is —8 
+5 e.u./mole. Transfer with monomer limits the number-average 
degree of polymerization to values of 40 or less except in special 
circumstances. Dimerization and double bond migration occur to 
some extent; n-perfluoropentadiene-], 3 is formed in the latter 
process. It and the 1,4-diene copolymerize; the latter undergoes 
cyclic addition so that the polymers are soluble and have little per- 
fluorovinyl unsaturation. The polymers are’ brittle if the fraction 
of 1, 3-diene in the polymer is less than 0.1. They are rubbery and 
of considerably higher molecular weight if the fraction of 1, 3- 
diene is greater than 0.4. The thermal stability of the polymers 
decreases as the content of 1, 3-diene increases. 
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Strength of fingernails, R. W. Young, S. B. Newman, and R. J. 
Capott, J. Invest. Dermatol. 44, No. 5, 358-360 (Dec. 23, 1964). 
A method is described for determining strength values on finger- 
nails. Flexural stresses ranging up to 880 kilograms per square 
centimeter have been obtained by this technique. Data from finger- 
nails of a number of subjects are shown to be in agreement with 
strength values for other keratinous materials. The use of this 
method as a tool for evaluating the effects of diet on fingernails is 
proposed. 


A miniature helium turbo-expander for cryogenic refrigera- 
tion systems, M. T. Norton, (Proc. XI Intern. Congress Refrigera- 
tion, Munich, Germany, Aug. 1963), Prog. Refrig. Sci. Technol. 
1, 131-135 (1963). 

The Cryogenic Engineering Laboratory is investigating the use of 
small turbo-expanders, supported by gas bearings, which produce 
small amounts of refrigeration at 4.2 °K and 30 °K. The turbine 
and externally pressurized gas bearings used with the 4.2 °K sys- 
tem have previously been reported. 

This paper describes a helium turbine expander designed for use 
in a closed cycle system. The turbine was designed to extract 
about 200 watts at 30 °K from the helium flow. The turbine is 7.9 
mm in diameter, rotates at 9800 revolutions per second, and was 
designed for an inlet pressure of 4 atmospheres and a 4:1 pressure 
ratio. Tests covering pressure ratios from 4:1 to 6.4:1 indicate 
a peak efficiency of 65 percent. The turbine design philosophy 
and performance data of this turbine are given. 


A waveguide noise standard, C. K. S. Miller, W. C. Daywitt, 
and E. Campbell, (Proc. Intern. Measurements Conf.), Acta IMEKO 
3, 371-382 (1964). 

Radio astronomy and noise measuring systems have laid heavy 
stress on the need for noise standards covering an increasingly 
larger frequency range. In response to this need, the National 
Bureau of Standards recently has completed work in National 
Standards and comparison systems for calibrating microwave noise 
sources in WR90 (8.2 to 12.4 Gc) and WR62 (12.4 to 18.0 Ge) wave- 
guide. This availability has, in turn, placed more emphasis on 
the need for good interlaboratory noise standards. The gas- 
discharge noise tube makes a very stable, long-life source for use 
in an interlaboratory standard. However, the available noise 
source mounts examined did not have the required characteristics 
for use in this quality of standard. This paper discusses experi- 
mental data of gas-discharge noise sources gathered over a period 
of two years and indicates where improvements can be made. 
Using this knowledge, an experimental gas-discharge noise tube 
mount was designed and constructed. The data taken using this 
mount indicated that the new design was suitable for interlaboratory 
work and provided a possible means of using the improved mount 
with different gas tubes to form more than one standard. 


A standard for extremely low values of spectral irradiance, 
R. Stair, W. B. Fussell, and W. E. Schneider, Appl. Opt. 4, No. 1, 
85-89 (Jan. 1965). 

To serve the interests of workers in certain fields of extremely 
low irradiation such as is present in the air glow and in phosphores- 
cent, fluorescent, biochemical, and associated areas, a standard of 
spectral irradiance of extremely low intensity has been setup. This 
standard, in effect, combines the NBS standard of spectral radiance 
with a small aperture which approximates a point source of radiant 
energy of known spectral distribution. This combination, with one 
or more convex mirrors, will supply extremely low, spectral ir- 
radiances in wavelength range 0.25 to 2.6 microns with an uncer- 
tainty of about 10 percent. 


Commercial carbon composition resistors as pressure trans- 
ducers, C. E. Miller, J. W. Dean, and T. M. Flynn, Rev. Sci. Instr. 
36, 231-232 (Feb. 1965). 

This paper presents the results of a recent series of pressure tests 
conducted on a group of commercial carbon composition resistors 
and discusses the possible application of the unique piezoresistive 
behavior of these resistors to pressure measurement. 


Human engineering in the design of a console for the com- 
parison of volt boxes, P. H. Lowrie, Jr., 19th Annual Instr. Soc. 


Am. Conf. and Exhibit, Preprint No. 12.1-3-64, pp. 1-8 (Oct. 12- 
15, 1964). 

Too often the human element is not given adequate consideration 
in the design of precise measuring equipment. The dependence 
upon the operator for the attainment of high measurement accuracies 
cannot be denied. This paper describes a calibration console in 
which the human engineering factors were given the same considera- 
tion as the technical requirements. In this system the operator is 
considered to be a decision maker, and those functions not requiring 
judgment are handled automatically by the console. Most of the 
calculations are performed automatically by internal circuits, and 
the results are displayed digitally upon command. The console is 
an illustration in point that the principles of human engineering are 
as applicable to the precise measurements found in the standards 
laboratory as they are to the production line. The paper discusses 
the features of the console and points out some of the pitfalls that 
may be encountered in a design of this type. 


Maintenance of maring electronic equipment by personnel 
not trained in electronic maintenance, G. Shapiro and G. J. 
Rogers, Proc. Annual Tanker Conf. American Petroleum Institute, 
pp. A-1—A-18 (June 1965). 

The fault isolation technique developed at the National Bureau of 
Standards for the Navy could be used aboard the ships of the Mer- 
chant Marine to permit the electronic equipment used on the bridge 
and in the engine room to be maintained at sea by the ship's officers. 
A more sophisticated version of the test instrument can be used to 
provide a maintenance log. 

This maintenance technique is implemented by building the elec- 
tronic system of functional modules, by providing test sockets at 
which each module’s performance can be evaluated, and by using 
a test instrument which indicates that module performance is good 
if it is within prescribed limits, and bad if it is not. For the main- 
tenance log, a different test instrument would be used to present 
performance information in quantitative terms. 


Optical instrumentation for the biologist: microscopy, 
J. R. Meyer-Arendt, Appl. Opt. 4, No. 1, 1—9 (Jan. 1965). 
A great variety of optical instruments has been developed for use in 
biology. The microscope and its many special types are part of 
them. In this article, a survey is given on the history, development, 
current status, and future outlook of general microscopy, micro- 
spectroscopy, interference microscopy, and electron and x-ray 
microscopy. 


Precision detector for complex insertion ratio measuring 
systems, C. M. Allred and R. A. Lawton, JEEE Trans. Instr. Meas. 
IM-13, No. 2 & 3, 76-81 (June-Sept. 1964). 

A precision, low-noise, highly sensitive, synchronous rf detector is 
described that permits simultaneous phase and level nulling in 
complex insertion ratio measuring systems. This insertion ratio is 
a measure of both the phase shift and magnitude change produced 
by a network when inserted in a transmission path. 


Stress corrosion cracking of metals, H. L. Logan, Metals Engr. 
Quart. 5, No. 2, 32-38 (May 1965). 

Stress-corrosion cracking results from a combination of high, but 
usually residual, tensile stresses and the chemical action of a speci- 
fic but mild corrodent. The functions of the specific corrodcnt are 
to (a) form a protective film on the surface of the metzi and (b) 
reduce the energy to form new surfaces at the crack tip. Cracking 
is initiated and propagated in part by an electrochemical process. 
Small anodes are formed at grain boundaries or on grain faces as the 
result of plastic deformation which either ruptures a protective 
film and/or produces a heterogeneous condition at the grain boundary 
or on the grain face. The unstrained material forms a large cathode. 
Cracking proceeds by an electrochemical process until the stress 
concentration becomes sufficiently large to trigger a fracture. In the 
early stages of crack propagation, fracture will be stopped by the 
energy requirements to plastically deform the material at the crack 
tip. Eventually the structure may fail by ductile fracture. 

If upon deformation tangled arrays of dislocations are formed in the 
material, stress-corrosion cracking, if it develops, will be inter- 
crystalline. 

In face centered cubic material, if conditions of dislocation array, 
short range order and/or stacking fault energy prevent cross slip, 
stress corrosion cracking will most probably be transcrystalline. 
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The effect of a fatigue crack on the fatigue strength of an 
aluminum alloy, J. A. Bennett, Mater. Res. Std. 5, No. 5, 235-239 
(May 1965). 

Specimens of 2024-T4 aluminum alloy were subjected to reversed 
bending at a stress amplitude of 40 ksi, and the loss of fatigue 
strength due to this stressing was evaluated by subsequent tests at 
lower stress amplitudes. The phenomenon of gas evolution was 
used to determine whether or not a crack had developed during the 
damaging stressing: it was found that prior to the development of a 
crack the fatigue strength was reduced only 10-20 percent but that 
a small crack caused a further 50-60 percent reduction. These 
results demonstrate the necessity for considering the two stages 
of a fatigue failure separately in attempting to evaluate cumulative 
damage in this material. 


The role of cryogenics in the production of high and ultra- 
high vacuum, L. O. Mullen and M. J. Hiza, Cryogenics 4, No. 6, 
387-294 (Dec. 1964). 

Many of the variables affecting the equilibrium pressure in a vacuum 
system are temperature dependent. This paper, oriented toward 
unbakeable and elastomer sealed systems, reviews the effects of 
some of these variables. Some recognized limitations of cryopumps 
and diffusion pumps at pressures in the low 10-' torr range are 
pointed out. 


A comparison between the National Bureau of Standards 
two-pressure humidity generator and the National Bureau of 
Standards hygrometer, S. Hasegawa. R. W. Hyland, and S. W. 
Rhodes, (Proc. 1963 Intern. Symp. Humidity and Moisture, Wash- 
ington, D.C.), Book, Humidity and Moisture. Measurement and 
Control in Science and Industry, Ed. A. Wexler, Vol. 3, Funda- 
mentals and Standards, Chapt. 22, pp. 455-459 (Reinhold Publ. 
Corp., New York, N.Y., 1965) 

4 3X3 Graeco-latin square experiment was used for testing four 
parameters of the National Bureau of Standards two-pressure 
humidity generator to determine which parameters influence the 
accuracy of the generator. The external saturator temperature, the 
saturation temperature, the saturator pressure, and the air flow were 
varied and the calculated mixing ratio of the generator was compared 
against the mixing ratio as measured with the National Bureau of 
Standards’ standard gravimetric hygrometer. 

The results show that for the highest accuracy, the National Bureau 
of Standards’ pressure humidity apparatus has a maximum flow 
limitation of under five standard cubic feet per minute (scfm). 
Tentative values of the correction factor for the saturation pressure 
of pure water vapor when used to compute the humidity of water-air 
mixtures were obtained. The correction factor is a function of 
saturation temperature and pressure. It is estimated that the 
maximum uncertainty in the. mixing ratio produced by the generator 
is 1/2 percent over a range of final saturator temperature of 25 
to — 20 °C and final saturator pressure of 0 to 60 psia. 


A pneumatic bridge hygrometer for use as a working humid- 
ity standard, L. Greenspan, (Proc. 1963 Intern. Symp. Humidity 
and Moisture, Washington, D.C.), Book, Humidity and Moisture. 
Measurement and Control in Science and Industry, Ed. A Wexler, 
Vol. 3, Fundamentals and Standards, Chapt. 20, pp. 433-443 
(Reinhold Publ. Corp., New York, N.Y., 1965). 

A pneumatic bridge, utilizing two sets of critical flow nozzles in series 
> and arranged to jorm a network analogous to a Wheatstone bridge, 
is described. By inserting a desiccant between the two nozzles in 
one arm of the bridge and measuring the ensuing pressure difference 
across the bridge produced by the absorption of the water from the 
test gas in one arm of the bridge, a measure of the humidity of the test 
gas is obtained. The equation expressing the relationship between 
the bridge pressure ratio, 1, and the mixing ratio, r, defined as the 
mass of water per mass of dry gas, is given by 


— _1+Kr 
ait Ud +r)/2(1 +r/E)'? 

where E is the ratio of the molecular weight of water to that of the 
dry test gas and K is an empirically determined constant. The bridge 
pressure ratio is the ratio of a mecsured pressure difference to the 
sum of two measured absolute pressures. 


This instrument may be calibrated against a suitable reference 
standard. When calibrated against the NBS gravimetric hygrom- 
eter, a value of K=20.62 X10-° was obtained. The largest error 
measured, over the range of 0.15 to 19.1 mg of water vapor per 
gram of dry air, was 0.025 mg of water vapor per gram of dry air. 
The upper limit to the effect of random errors and possible systematic 
errors in terms of mixing ratio is estimated at 0.050 mg of water 
vapor per gram of dry air or 1/4 percent of full scale. The instru- 
ment was designed for use as a working or transfer standard for 
humidity measurements. 


The international system of units, A. G. McNish, Am. Soc. 
Quality Control Annual Tech. Conf. Trans. pp. 431-438 (1965). 
The International System of units is advantageous in engineering 
calculations because of its completeness and simplicity. No difh- 
culty is experienced in its application in the United States although 
it is inconsistent with our customary units for length, mass and 
temperature. 


The NBS standard hygrometer, A. Wexler and R. W. Hyland, 
(Proc. 1963 Intern. Symp. Humidity and Moisture, Washington, 
D.C.), Book, Humidity and Moisture. Measurement and Control 
in Science and Industry, Ed. A. Wexler, Vol. 3, Fundamentals and 
Standards, Chapt. 19, pp. 389-431 (Reinhold Publ. Corp., New 
York, N.Y., 1965). 

A gravimetric hygrometer is described that serves as the NBS 
standard instrument for the measurement of the moisture content 
of gases on an absolute basis in terms of mixing ratio (mass water 
vapor/unit mass of dry gas). The National Bureau of Standards 
and other laboratory reference and working instruments are inter- 
compared and calibrated with this instrument. The measuring 
operation involves the absorption of the water vapor from a water 
vapor-gas mixture by a solid desiccant and the determination of the 
mass of this water vapor by precision weighing: it also involves 
the determination of the volume of the associated gas of known 
density by counting the fillings of two calibrated stainless steel 
cylinders. An automatic system permits the sampling of the test 
gas at any desired flow rate up to 2 lpm (STP) and for any desired 
number of fillings. The instrument provides a value of the mixing 
ratio averaged over the time interval of a test. 

The construction and operation of the instrument is described. 
Discussions of the tests and calibrations of component parts, and of 
the sources of errors also are included. An analysis of the random 
and systematic errors effecting the overall accuracy in the deter- 
mination of mixing ratio shows that if 0.60 g of water vapor is col- 
lected from moist air, then the estimated maximum uncertainty 
expected for mixing ratios between 27 mg/g and 0.19 mg/g is 12.7 
parts in 104. 


Characteristics of a Raman laser excited by an ordinary 
ruby laser, H. Takuma and D. A. Jennings, Proc. IEEE 53, No. 2, 
146-149 (Feb. 1965). 


The stimulated Raman effect of benzene has been observed using 


an ordinary (non-giant) ruby laser. The buildup of oscillation at 
the ve, 2v2, 3v2, and 4v. Stokes lines and also at the first v; Stokes 
line and also at the first vy, Stokes line have been observed. The 
threshold exciting power for laser action in the v, Stokes line has 
been measured to be 9.5 kW. A rate equation for the Raman laser 
has been given, and the constant for Raman scattering of benzene 
has been determined as k= 9.5- 10-25 sec.-!. The estimation based 
on the results of this investigation indicates that it is possible to 
construct a Raman laser of benzene using a He-Ne gas laser. 


Drift-free Mossbauer spectrometer, F. C. Ruegg, J. J. Spijker- 
man, and J. R. DeVoe, Rev. Sci. Instr. 36, No. 3, 356-359 (Mar. 
1965). 

In order to eliminate drifts encountered in electro-mechanical 
Mossbauer spectrometers, a drive has been developed which is 
drift free with regard to velocity as well as position. The spec- 
trometer can be driven in two modes: constant acceleration and 
constant velocity. The constant acceleration input signal is the 
analog output of the address scalers of a multichannel analyzer 
running as a multi-scaler (time mode). The constant velocity is 
obtained from a photo-cell slit system which drives the input of the 
spectrometer with a rectangular wave. In both modes, the photo- 
cell system is used to eliminate all position drifts by starting each 
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new cycle from the same position. To ar ‘ ze the velocity of the 
spectrometer a voltage to frequency couverter and multichannel 
analyzer are used. 


Elastic constants of cubic lead fluoride at room temperature, 
J. H. Wasilik and M. L. Wheat, J. Appl. Phys. 36, No. 3, Pt. 1, 791- 
793 (Mar. 1965). 
Using an ultrasonic pulse echo technique operating at frequencies 
10 Me/sec, we have obtained the following elastic constants for 
cubic PbF.: C,,=8.88 x 10" dyn/em?, Cy. = 4.72 x 10" dyn/cem?, 
Cy, = 2.454 x 10" dyn/cm? at 27 °C. From these data a Debye 6 
of 221.2 °K was calculated. 


Electron microscopy and diffraction of synthetic corundum 
crystals. I. Pure aluminum oxide grown by the Verneuil 
process, D. G. Barber and N. J. Tighe, Phil. Mag. 11, No. 111, 
495-512 (May 1965). 

Single crystal aluminum oxide has been chemically thinned by 
means of a jet-polishing device and examined in the electron micro- 
scope. During thinning, the crystals could be etch-pitted, and the 
relationship between etch-pits and dislocations has been investi- 
gated; etch-pits on subgrain boundaries are often unrepresentative 
of the dislocation density. The characteristics of thermally etched 
surfaces have also been studied, and the results of evaporation ex- 
periments in the microscope are presented. 

The behaviour of dislocations in various plastically deformed crys- 
tals has been examined. Tangles of dislocations, commonly found 
in metals and in other ionic solids, are absent, and dislocation 
densities in glide bands are low. Grown-in dislocations apparently 
participate in the glide process: therefore, they are not strongly 
pinned by impurities at temperatures permitting plastic deformation. 
The occurrence of dislocation dipoles is taken as evidence for the 
existence of cross slipping of basal dislocations at temperatures 
below those required for slip on {1120} prism planes. No positive 
evidence for dislocation dissociations have been found and if split- 
ting occurs, it must be small. indicating that the stacking fault 
energy is high. 

Preliminary studies of dislocation motion in the microscope at 
temperatures ~ 1500 °C are reported and discussed. 


Ice nucleation by photolyzed silver iodide, G. Burley, Phil. 
Mag. 10, 527 Sept. 1964). 

\ progressive change of the shape of ice crystals nucleated by silver 
iodide particles at constant temperature has been found, with a 
linear dependence only on the cumulative time of exposure of the 
latter to light. The directional crystallographic properties of the 
photodeactivation process are applied to correlate this with the 
mechanism of epitaxial nucleation. 


The absolute photometry of the zodiacal light, L. L. Smith, 
F. E. Roach, and R. W. Owen, Planet. Space Sci. 13, 207-217 
(Pergamon Press Ltd., Northern Ireland, 1965). 

Interpretation is made of photometric records obtained at the 
Haleakala Observatory during six nights between August 1961 and 
June 1962. An interference filter centered on wavelength 5300A 
was used. After allowance for the airglow continuum and the inte- 
grated starlight the residual was interpreted as due to zodiacal light. 
The brightness of the zodiacal light is tabulated for that part observ- 
able in the night sky. Graphical representations of the data illus- 
trate the distribution of brightness in ecliptic coordinates. Evidence 
is presented supporting the existence of a terrestrial component 
superimposed photometrically on the general interplanetary 
component. 


X-ray spectrochemical analysis of material: cement and 
dental alloys, B. L. Bean and B. W. Mulligan, Am. Soc. Testing 
Vater. Spec. Tech. Publ. 373, Application of Advanced and Nuclear 
Physics to Testing Materials, pp. 25-45 (June 17, 1965). 

The application of X-ray spectrochemical analysis to various ma- 
terials has undergone marked expansion in recent years. Special 
examples from the authors’ experiences in the analysis of cements 
and dental materials are cited. Difficulties in presenting a repre- 
sentative sample to the primary X-ray beam are described, together 
with means for minimizing these difficulties. Interelement cor- 
rections for the X-ray analysis of cement are discussed. References 
are given to the work of numerous investigators in this field. 


Measurement and standardizatiun of LF to UHF electrical 
quantities. Sketch of recent world-wide developments, 
M. C. Selby, Prog. Radio Sci. 1960-1963, Vol. 1. Radio Standards 
and Measurements, pp. 45-55 (Elservier Publ. Co., Amsterdam, The 
Netherlands, 1965). 

An outline of world-wide activities and progress on measurement & 
standardization of electrical quantities over an approximate range 
of 30 Ke/s to 3 Ge/s is presented. The period covers approximately 
the years of 1960, 1961, and 1962. More prominent technical 
developments are noted and publication subjects in a number of 
countries are indicated to expose the apparent areas of more specific 
interest in these countries. Accelerating national and international 
organizational activities are also called attention to in order to show 
the degree of increasing attention to this subject. An attempt is 
made to present a digest which will be found more specific and more 
revealing than can be expected from a report of the International 
Chairman of URSI Commission I. On the other hand, it is far from 


the detailed account to be found in the national progress reports. 
The information on international organizational activities seems to 
be pointed out for the first time in these series of URSI reports. 


Radio meteorology, J. W. Herbstreit, Science 147, No. 3653. 
76-78 Yan. 1, 1965). 

A report of the highlights of the 1964 World Conference on Radio 
Meteorology held at the National Bureau of Standards, September 
14-18, 1964 


The use of electromagnetic signals emitted from nuclear 
explosions to study long-range VLF propagation, A. G. Jean 
and J. R. Wait, J. Geophys. Res. 70, No. 5, 1258-1251 (Mar. 1, 1965). 
Within recent years a number of reports [Newman, 1953: Kom- 
paneets, 1958; United Nations, 1958] have been published which 
describe the electromagnetic radiation from a nuclear detonation. 
The phenomenon has been discussed at the conference of “experts” 
at Geneva at various intervals since the summer of 1958 [United 
Nations, 1958}. 


Other NBS Publications 


Journal of Research 69A (Phys. and Chem.), No. 4 (July- 

Aug. 1965), 70 cents. 

Correlation of successive atomic steps in crystals by relaxation 
mode analysis. A. D. Franklin. 

Reduction of crystallographic point groups to subgroups by homoge- 

H.S. Peiser and J. B. Wachtman, Jr. 

Effect of hydrostatic pressure on the refractive indices of some solids. 
R. M. Waxler and C. E. Weir. (See above abstracts). 

Crystallization kinetics and polymorphic transformations in poly- 
butene-l. J. Powers, J. D. Hoffman, J. J. Weeks, and F. A. 
Quinn, Jr. 

Compliance-time-temperature — relationships from 
measurements on a pure-gum rubber vulcanizate. 
G. W. Bullman, and L. A. Wood. 

Distribution function of the end-to-end distances of linear polymers 
with excluded volume effects. J. Mazur. 

An additivity rule for the vapor pressure lowering of aqueous solu- 
tions. R. A. Robinson and V. E. Bower. 

Franck-Condon factors to high vibrational quantum numbers V: 
QO, band systems. R. W. Nicholls. 

An absolute temperature scale from 4 °K to 20 °K determined from 
measurements with an acoustical thermometer. H. H. Plumb 
and G. Cataland. (See above abstracts). 


neous stress. 


indentation 


F. L. Roth, 


Journal of Research 69A (Phys. and Chem.), No. 5 (Sept.- 

Oct. 1965), 70 cents. 

Current research on preservation of archival records on silver- 
gelatin type microfilm in roll form. C.S.McCamy and C. I. Pope. 
(See above abstracts). 

Franck-Condon factors to high quantum numbers VI: C. band 
systems. R. W. Nicholls. 

Matrices of spin-orbit interaction in the electron configurations 
p?d and p*d. J. L. Tech and R. H. Garstang. 

Heat of formation of calcium aluminate tricarbonate at 25 °C. 
Berman. 

Vapor pressure and heat of sublimation of tungsten. 
E. R. Plante, and J. J. Diamond. 


BAA. 


R. Szwarce, 
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Heat capacity and enthalpy measurements on aluminum carbide 
(Al,Cs) from 15 to 1173 °K. Thermodynamic properties from 0 
to 2000 °K. G. T. Furukawa, T. B. Douglas, W. G. Saba, and 
A. C. Victor. 

Properties of aqueous mixtures of pure salts: thermodynamics of 
the ternary system water-potassium chloride-barium chloride 
at 25 °C. R. A. Robinson and V. E. Bower. 

A wide-range (up to 10'° P) rotating cylinder viscometer. A. Napolli- 
tano, P. B. Macedo, and E. G. Hawkins. (See above abstracts). 

Internal friction in ZrO. containing CaO. J. B. Wachtman, Jr.. 
and W. C. Corwin. (See above abstracts). 

Splitting of equivalent points in noncentrosymmetric space groups 
into subsets under homogeneous stress. H. S. Peiser, J. B. 
Wachtman, Jr., F. A. Munley, and L. C. McCleary. 

Inclusion of perfluoromethyl groups in the lattice of copolymers 
of tetrafluoroethylene and hexafluoropropylene. L. H. Bolz and 


R. K. Eby. 


Journal of Research 69B (Math. and Math. Phys.), No. 3 

(July-Sept. 1965), 75 cents. 

A development of the theory of errors with reference to economy 
of time. M. D. Hersey. Prefatory note by Churchill Eisenhart. 

Transversals and matroid partition. J. Edmonds and D. R. Ful- 
kerson. 

Some L, Markoff inequalities. L. F. Shampine. 

Some theorems on the permanent. R. A. Brualdi and M. Newman. 

On Kirchhoff’s law and its generalized application to absorption and 
emission by cavities. F. J. Kelly. (See above abstracts). 

On a relation between two-dimensional Fourier integrals and series 
of Hankel transforms. J. V. Cornacchio and R. P. Soni. 

C. Witzgall. 

Some extensions of Banach’s contraction theorem. P. R. Meyers. 

A variant of the two-dimensional Riemann integral. A. J. Goldman. 

The use of finite polynomial rings in the factorization of the general 
polynomial. D. B. Lloyd. 

A primal (all-integer) integer programming algorithm. 


On convex metrics. 


R. D. Young. 


Radio Science J. Research NBS/USNC-URSI 69D, No. 9 
(Sept. 1965), $1.00. 
Ground-based passive probing using the microwave spectrum of 


oxygen. E.R. Westwater. (See above abstracts). 


Response of NBS microwave refractometer cavities to atmospheric 


variations. R.O. Gilmer. R 
above abstracts). 


E. McGavin, and B. R. Bean. 


(See 

Effects of rocket outgassing on ri experiments. W. Pfister and J. C. 
Ulwick. (See above abstracts). 

Further analysis of propagation of plasma waves in a “spoke-wheel” 
magnetic field. R. R. Gold 

Measurement of group velocity of 17.8-kHz VLF radio waves. W.D. 
Westfall. (See above abstracts). 

On the radio noise level at low and very low frequencies in polar 
regions. T.S. Jorgensen. 

Experimental study on the circular loop antenna immersed shallowly 
in a conducting medium. K. lizuka. 

Directivity of uniformly spaced optimum endfire arrays with equal 
sidelobes. M. T. Ma. 

Extension of Fock theory for currents in the penumbra region. 
V. H. Weston. 

Isometric-circle interpretation of bilinear transformation and its 
application to VSWR minimization. J. G. Rudolph and D. K. 
Cheng. 

Magneto-ionic propagation in inhomogeneous media. Part I. 
Transverse propagation. B. K. Banerjea. 

Magneto-ionic propagation in inhomogeneous 
Oblique propagation. B. K. Banerjea. 

Editorial comment on the scientific papers of Lord Rayleigh (John 
William Strutt). J. R. Wait. 


media. Part II. 


Radio Sci. J. Research NBS/USNC-URSI 69D, No. 10 (Oct. 

1965), $1.00. 

Irreversible power and radiation resistance of antennas in aniso- 
tropic ionized gases. K. S.H. Lee and C. H. Papas. 

Scattering resonances of a cylindrical plasma. W. M. Leavens. 

Radiation patterns from plasma enclosed cylindrical hypersonic 
vehicles. J. H. Harris, A. T. Villeneuve, and L. A. Broca. 


The Schumann resonances. R. K. Cole, Jr. 

Atmospheric radio noise bursts in the LF band at Bangalore. 
C. Aiya and K. N. Lakshminarayan. 

Influence of finite ground conductivity on the propagation of VLF 
radio waves. J. R. Wait and K. P. Spies. 

Model experiments on propagation of groundwaves across an abrupt 
boundary at perpendicular incidence. R. J. King and S. W. Maley. 


Bam. 


Radio Science J. Research NBS/USNC-URSI 69D, No. 11 

(Nov. 1965), $1.00. 

Propagation of pulses in dispersive media. 
abstracts). 

An anisotropic electron velocity distribution for the cyclotron ab- 
sorption of whistlers and VLF emissions. H. Guthart. 

Nose whistler dispersion as a measure of magnetosphere electron 
temperature. H. Guthart. 

Interference rejection capability of a switched radiometer. 
Clapp. (See above abstracts). 

Atmospheric breakdown limitations to optical maser propagation. 
R. G. Tomlinson. (See above abstracts). 

Phase steps and amplitude fading of VLF signals at dawn and dusk. 
D. Walker. 

Propagation in a model terrestrial waveguide of nonuniform height 
theory and experiment. E. Bahar and J. R. Wait. 

Comments on H. Volland’s “Remarks on Austin’s 
J. R. Wait. 

The path integrals of LF/VLF wave hop theory. 
M. E. Chrisman. 

Reactive loading of arbitrarily illuminated cylinders to minimize 
microwave backscatter. K. M. Chen. 

On the statistical theory of electromagnetic waves in a fluctuating 
medium (I1) Mathematical basis of the analogies to quantum field 
theory (a digest). K. Furutsu. 


J. R. Wait. 


(See above 


R. E. 


Formula.” 


L. A. Berry and 


lonospheric radio propagation, K. L. Davies, NBS Mono. 80 (Apr. 
1965), $2.75. 

Systems engineering in ceramics. Proceedings of a symposium 
April 19, 1964, NBS Misc. Publ. 267 (May 1, 1965), $2.00. 

Transmission predictions for tropospheric communication 
circuits. P. L. Rice, A. G. Longley, K. A. Norton, and A. P. 
Barsis, NBS Tech. Note 101, Vol. I (May 7, 1965), $1.00. 

Transmission loss predictions for tropospheric communication cir- 
cuits. P. L. Rice, A. G. Longley, K. A. Norton, and A. P. Barsis, 
NBS Tech. Note 101, Vol. Il (May 7, 1965), $1.00. 

Research on crystal growth and characterization at the National 
Bureau of Standards July to December 1964, ed. H. F. McMurdie, 
NBS Tech. Note 260 (May 8, 1965), 50 cents. 

Studies of solar flare effects and other ionospheric disturbances 
with a high frequency Doppler technique, V. Agy, D. M. Baker. 
and R. M. Jones, NBS Tech. Note 306, (Apr. 28, 1965), 75 cents. 

Data reduction for stable auroral red arcs observed at Rapid City, 
South Dakota, J. E. Cruz, R. Davies, L. K. Droppleman, E. Maro- 
vich, L. R. Megill, M. H. Rees, L. Reisbeck, and F. E. Roach, 
NBS Tech. Note 308 (May 3, 1965), 50 cents. 

Attenuation of the ground wave of a low frequency electromagnetic 
pulse, J. C. Morgenstern and J. R. Johler, NBS Tech. Note 310 
(Apr. 1965). 30 cents. 

On the effect of heavy ions on LF propagation, with special reference 
to a nuclear environment, J. R. Johler and L. A. Berry, NBS Tech. 
313 (June 7, 1965), 45 cents. 

Stability of two-phase annular flow in a vertical pipe, S. Jarvis, Jr., 
NBS Tech. Note 314 (June 7, 1965), 55 cents. 

Selected tables of atomic spectra. Atomic energy levels and 
multiplet tables Sill, Siti, Sitv, C. E. Moore, NSRDS-NBS3, 
Sec. 1 (June 25, 1965), 35 cents. 


loss 


A line-formula notation system for coordination compounds, P. M. 
McDonnell and R. F. Pasternack, J. Chem. Documentation 5, 
No. 1, 56-60 (Feb. 1965). 

Decyclization of cyclohexene, 4-methylcyclohexene, and 4-vinyl- 
cyclohexene in a single-pulse shock tube, W. Tsang, J. Chem. 
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